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Abstract 



We study selected aspects of quantum gravity phenomenology inspired by the grav- 
itational analogy in Bose-Einstein condensates (BECs). We first review the basic 
ideas and formalism of analogue gravity in BECs, with particular emphasis on the 
possibility of simulating black holes. We stress that the dispersion relation in a BEC 
has a relativistic form in the hydrodynamic limit, and acquires non-relativistic, 'su- 
perluminal' modifications at high frequencies. This makes it a particularly interest- 
ing model for many scenarios of quantum gravity phenomenology which consider 
a possible violation of local Lorentz invariance at high energies. In particular, it al- 
lows the study of kinematical corrections that such quantum gravity scenarios could 
impose on general relativity. 

We present a simple model for a (l+l)-dimensional acoustic black hole configu- 
ration in a BEC, and study its dynamical stability. We find that black hole horizons in 
a BEC are stable under quite general boundary conditions, inspired by similar stud- 
ies of the stability of black holes in general relativity. This allows us to study the 
stable dynamical modes or quasinormal modes of such configurations, as well as to 
confirm the long-expected possibility of (at least in principle) simulating Hawking 
radiation in a BEC. 

We find that there indeed exist short-lived quasinormal modes in our model, in 
spite of its (l+l)-dimensionality, and moreover that their spectrum covers a contin- 
uous region of the complex frequency plane. This is in sharp contrast to general 
relativistic black holes and to acoustic black holes in the hydrodynamic limit, where 
no quasinormal modes exist in 1+1 dimensions, whereas in higher dimensions the 
spectrum is discrete. The quasinormal modes that we find are characteristic of a 
relaxation at the microscopic level, and we trace their presence back to the superlu- 
minality of the dispersion relation, and the associated permeability of the horizon. 
We argue that a similar effect could show up for astrophysical black holes if quan- 
tum gravity turns out also to lead to superluminality at high frequencies. 

We then study the impact of superluminal dispersion relations on the Hawking 
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radiation for a collapsing geometry in which a black hole is created. A crucial aspect 
of modified dispersion relations is that the horizon becomes a frequency-dependent 
concept. In particular, with superluminal dispersion relations, at every moment of 
the collapse process, there is a critical frequency above which no horizon is experi- 
enced. If moreover there is an overall saturation level to the collapse (as is certainly 
the case in BECs due to the interatomic separation, and as is expected in the gravi- 
tational case due to quantum gravity resolving the general relativistic singularity), 
then there is also a global critical frequency, such that higher frequencies never see 
a horizon and hence do not participate in the thermal radiation process. This leads 
to several important differences with the standard Hawking radiation. Generally 
speaking, a finite critical frequency implies that the intensity of the late-time radia- 
tion is weaker than in the standard (Lorentz invariant) case. The Planckian form of the 
thermal spectrum is largely preserved (although with a decreased intensity) unless 
the Lorentz symmetry violation scale (the scale at which non-relativistic deviations 
in the dispersion relation become important) lies below the critical frequency. In the 
latter case, high frequencies (above the Lorentz symmetry violation scale) can ac- 
quire extremely large surface gravities, leading to important ultraviolet contributions 
in the radiation spectrum. Finally, the radiation originating from the collapse dies 
out on a relatively short time-scale. 

In the final part of this thesis, we discuss some questions related to the possibil- 
ity of constructing a serious toy model for Planck-scale gravity based on the con- 
densed matter analogy, namely a model of emergent gravity. The basic idea is that 
gravity (including its dynamical aspects, i.e., the Einstein equations) might emerge 
as an effective low-energy macroscopic description of the collective behaviour of 
the microscopic constituents of a system. A crucial consequence if this turns out to 
be a correct approach is that there would not be any direct quantisation procedure 
leading from the macroscopic to the microscopic degrees of freedom. We mention 
some general motivations for such an approach, in particular the problem of dark 
energy, and discuss in detail the question of diffeomorphism invariance in emergent 
gravity based on condensed matter models. We argue that this question should be 
split up into a kinematical and a dynamical part. We explain how diffeomorphism 
invariance at the kinematical level can actually be understood as a low-energy ef- 
fective symmetry in such models, just like Lorentz invariance. With respect to the 
dynamical aspects of diffeomorphism invariance, we argue that this is essentially 
identical to the problem of recovering the Einstein equations in a condensed-matter- 
like model, and make some remarks with regard to the possibility that Sakharov's 
induced gravity proposal could solve this issue. 



Chapter 1 
Introduction 



The fundamental observation underlying the study of models for gravity in Bose- 
Einstein condensates (BECs) is the following [117j. The equation of motion of acous- 
tic perturbations in a perfect (irrotational, inviscid and barotropic) fluid (such as a 
BEC — under certain approximations which we will detail extensively in due time) 
is described by a d'Alembertian equation in curved spacetime: 

Op = -^=d }i ^igV v d v cp = . (1.0.1) 
v Q 

These acoustic perturbations or phonons therefore travel along the null geodesies 
of the effective metric gf„ v , with g its determinant. This formula is well known from 
relativistic field theory as the equation of motion for a massless scalar field (p prop- 
agating in a curved spacetime. Moreover, up to a conformal factor, the effective 
metric has the following form: 




(1.0.2) 



where v is the velocity vector of the background fluid flow and c the speed of sound 
in the fluidQ Since g tt oc (v 1 — c 2 ), the sign of g tt depends on the fluid regime, 
indicating that black hole configurations are in principle possible, at least in radial 
flows. 

This basic observation is the essence of analogue gravity in condensed matter 
models, which offers the promising prospect of allowing to better understand and 
possibly simulate some phenomena related to high-energy physics, and to black 

1 Throughout this thesis, we will only explicitly use a different notation for the speed of sound 
and the speed of light when this seems useful to avoid confusion. We will however always write c 
explicitly rather than setting c = 1. 
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holes in particular, as we briefly discuss in the following section. Its interpreta- 
tion can also be extended in two directions, which we will call quantum gravity 
phenomenology and emergent gravity, in logical order, and which we will introduce 
thereafter. 



1.1 Analogue gravity in condensed matter 



The basic observation described above implies that one can study certain aspects of 
general relativity and quantum field theory, and in particular of black hole physics, 
by analogy with fluid systems such as a BEC B56l |57| . To put it in other words: 
perfect fluids mimic the kinematical aspects of general relativity, i.e., the propagation 
of fields in curved relativistic spacetimes, including the possibility of black hole 
configurations. Note that there are currently no known analogue systems which 
reproduce the Einstein equations, which is why the analogy does not immediately 
extend to the dynamical level. 

To take maximal advantage of the analogy, the microscopic physics of the fluid 
system should be well understood, theoretically and experimentally, preferably even 
in regimes where the relativistic description breaks down. Then, full calculations 
based on firmly verified and controlled physics are possible, even beyond the rel- 
ativistic regime. Additionally, laboratory experiments become feasible that could 
shed light on issues of high-energy physics. There exists a wide variety of phys- 
ical systems, both condensed matter and other (see lfT5ll for a review), in which 
curved effective relativistic spacetimes emerge. BECs stand out mainly for their 
conceptual clarity, both as a source of the gravitational analogy (excitations moving 
in a mean-field background provided by the collective behaviour of the condensed 
phase of the condensate) and as a generally well understood and experimentally 
well controlled system, thereby offering the exciting prospect of directly simulating 
the analogues of some phenomena predicted in high-energy physics. For instance, 
BECs are considered a good candidate for a possible future experimental detection 
of (phononic) Hawking radiation |[T3l l35t |43| or of particle creation in expanding 
spacetimes Il69l 11331 . 
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1.2 From analogue gravity to quantum gravity 
phenomenology 

However, the analogy is not limited to the propagation of relativistic fields in curved 
background spacetimes. A particularly interesting application is in terms of quan- 
tum gravity phenomenology Indeed, the geometric or relativistic regime described 
above is not universally valid in a BEC, but only in a certain approximation. Es- 
sentially, it is valid for the low-energy modes or phonons, i.e., in the hydrodynamic 
approximation. For higher energies, and also in the presence of a horizon (as we 
will emphasise repeatedly throughout this thesis), the relativistic or hydrodynamic 
approximation breaks down, giving way to the underlying microscopic description. 
This behaviour is of course similar to what is expected in the realm of quantum grav- 
ity, with the crucial difference that the underlying microscopic description is well 
known in the case of BECs. Therefore, deviations from the hydrodynamic or relativis- 
tic regime in BECs provide possible insights into the kinematical corrections that a 
hypothetical quantum theory of gravity might impose on general relativity^] for ex- 
ample on the local Lorentz invariance that characterises general relativity. They 
thus provide a way of studying quantum gravity phenomenology by analogy. Such 
extrapolations from condensed matter models to quantum gravity phenomenology 
are especially relevant if the deviations from classical general relativity imposed by 
quantum gravity turn out to be roughly similar to those that occur in BECs with re- 
spect to the relativistic regime. In BECs, the effective low-energy Lorentz symmetry 
is broken as follows. The full dispersion relation between the frequency to and the 
wave number k is of the 'superluminal' (strictly speaking: supersonic) type: 



to 



2 = c 2 k 2 + <x 2 c 2 k A , (1.2.1) 



with a. a parameter (of dimension [a] = L) sufficiently small such that a relativistic 
dispersion relation to 2 = c 2 k 2 is recovered at low energies. Note again, as in the met- 
ric ( |1.0.2| >, that c represents here the (low-frequency) speed of sound. This dispersion 



relation indicates that the high-frequency modes move at a speed higher than the 
speed of sound. Essentially, the gravitational analogy in a BEC is then of particular 
relevance for quantum gravity phenomenology if quantum gravity (independently 
of its fundamental structure) modifies the observable low-energy physics in a way 
that can be written as a modified dispersion relation with respect to some preferred 



2 We again stress the kinematical aspect: the condensed matter analogy does not reproduce the 
Einstein equations, and hence does not allow to study quantum gravitational corrections to the dy- 
namics of general relativity beyond the linear regime. 
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reference framej^ This is obviously an experimental issue which is far from being 
settled. Ultimately, it means that the condensed matter model is not a good phe- 
nomenological model if quantum gravity turns out to be strictly Lorentz invariant, 
although even then it could still teach us valuable lessons about the implications of 
this strict invariance. It is perhaps also of limited use for Doubly Special Relativity 
and related types of Lorentz modification which can also lead to modified disper- 
sion relations, but in general of a more complicated form and, more importantly, 
without breaking the equivalence of reference frames [2j. However, it is certainly a 
natural model for effective field theory approaches in general, which form the bulk 
of current efforts at quantum gravity phenomenology Il3ll85ll. 

Quantum gravity phenomenology from condensed matter models is therefore 
an illustration of what one could call a 'bottom-up' approach to quantum gravity: a 
study of the first deviations that could occur with respect to the known low-energy 
(Lorentz invariant) laws of physics of our universe, based on the kind of deviations 
that occur in laboratory systems which reproduce low-energy effective Lorentz in- 
variant spacetimes similar to the one of our universe. This approach is therefore 
complementary to the traditional 'top-down' approaches to quantum gravity such 
as string theory or loop quantum gravity, which postulate a basic structure for our 
universe or its geometry (or a basic set of mathematical guidelines for arriving at 
such a structure) and then attempt to connect these basic principles with our known 
physics. Recovering the (semi)classical limit for our universe (i.e., the "reconstruc- 
tion problem", as it is called in [32]) is precisely one of the main open problems in all 
such top-down approaches. Actually, at present, it is not clear whether the semiclas- 
sical limit is uniquely defined for any of those top-down approaches, and so none of 
them make any falsifiable prediction with respect to modified dispersion relations, 
or even with respect to Lorentz invariance in general. Therefore, to bridge the gap, 
it should be clear that there is a genuine interest for a bottom-up approach such as 
the one discussed here. 



1.3 From analogue gravity to emergent gravity 

Finally, the condensed matter analogy with gravity can be taken another step fur- 
ther. If relativistic spacetimes naturally emerge in a multitude of condensed matter 
systems, then perhaps gravity as such (now, including its dynamical aspects) is ac- 
tually a generic property of many physical systems, and not necessarily limited to 

3 Note that we use the term 'modified dispersion relation' in the sense traditional in quantum 
gravity phenomenology, i.e.: modified with respect to a relativistic dispersion relation. 
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the basic structure of our universe. In other words, gravity might be an effective 
description emerging at the macroscopic level as a consequence of the collective be- 
haviour of the microscopic constituents of a system. It could even emerge (in some 
limit) in a wide variety of microscopic systems. If this turns out to be a correct de- 
scription of gravity, then one important consequence would be that there is no direct 
quantisation procedure leading from the macroscopic to the microscopic degrees of 
freedom. Part of the motivation for such an approach is to be found in the previ- 
ous paragraph: contrarily to the case in top-down approaches for quantum gravity, 
the semiclassical limit is well defined in the condensed matter model. Although 
not central to this thesis, we will also briefly comment (in chapter |8| on a second 
important motivation, namely the dark energy problem. 

However, such an approach also raises a number of important questions. First, 
as observed repeatedly above, the gravitational analogy in condensed matter sys- 
tems is currently limited to the kinematical aspects of general relativity, i.e., to the 
propagation of (classical or quantum) fields in curved relativistic spacetimes. This 
seems to be an important limitation: any serious toy model for a quantum theory 
of gravity should lead to the Einstein equations in the adequate regime. There- 
fore, a crucial question in this respect is the following. Does the condensed matter 
analogy simply break down at the observation that the Einstein equations do not 
reproduce in a laboratory, or can we take this further and understand precisely why 
our universe is ruled by the Einstein equations while laboratory condensed mat- 
ter systems are not? A second question is related to diffeomorphism invariance and 
background independence. The effective geometries that emerge in condensed mat- 
ter models are obviously background dependent: they depend on the presence of 
the condensed matter atoms, which are localised in the geometry of the laboratory. 
When extending this to a toy model for the gravitation of our universe, one is natu- 
rally led to a structure reminiscent of the ether theories that were popular before the 
advent of special relativity ||2T1 |49i whereby the 'spacetime atoms' are localised in 
an absolute predefined geometry. In particular, this seems to suggest that emergent 
gravity models inspired by the condensed matter analogy are not diffeomorphism 
invariant in general, contrarily to other approaches such as loop quantum gravity, in 
which diffeomorphism invariance is regarded as a crucial requirement for quantum 
gravity MM- 

We will argue that the problem of diffeomorphism invariance can be separated 
in a kinematical part and a dynamical part. Since we have repeatedly stated that 
condensed matter models reproduce the kinematical aspects of general relativity, it 
perhaps comes as no surprise that the kinematical aspect of the diffeomorphism in- 
variance problem has a solution, which actually involves more conceptual subtlety 
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than technical machinery. The dynamical aspect of the diffeomorphism invariance 
question, however, is very much an open problem. We will limit ourselves to point 
out that it is intimately related to the problem of reproducing the Einstein equations, 
and discuss a possible way to tackle this problem, based on Sakharov's induced 
gravity [106], together with some of its problems. 



1.4 A very brief history of analogue gravity 

In this section, we will mention some historical landmarks in the development of 
analogue models of gravity with particular relevance to the present thesis. A com- 
prehensive review with an extensive history section can be found in IITBll . 

Although it had already been known for several years that relativistic metrics 
such as ( |1.0.2[ > could be written explicitly for the propagation of sound waves in 



inviscid fluids [135], the original observation that this allows (acoustic) black hole 
configurations, and thereby the possibility of studying high-energy (quantum grav- 
itational) effects on black hole physics is due to Unruh in a seminal paper published 
in 1981 fll 14] | • Perhaps curiously, this paper received little attention, until Jacobson 
took up Unruh's suggestion to study the transplanckian problem in Hawking ra- 
diation ten years later [65]. This led to a burst of interest and the appearance of 
several papers on the robustness of Hawking radiation (see chapter [7]), all based 
on more or less abstract models for the transplanckian physics. In a crucial paper 
from 1998 [ 1117] , Visser established the gravitational analogy on firm foot by deriv- 



ing (1.0.1) and ( 1.0.2| > as a rigorous theorem of mathematical physics and studying 



in detail the analogue equivalents of quantities associated with black holes such as 
the horizon, ergosphere and surface gravity. 

In parallel, several proposals were made for concrete physical systems in which 
black hole analogies could be simulated, including (see [15] and references therein) 
solid state models, slow light in dielectrics, gravity waves in shallow water, and 
two key proposals in the area of condensed matter: a series of works by Volovik 
on gravitational analogies in 3 He-A, culminating in the appearance in 2003 of the 
book |124[| . and the proposal by Garay et al. in 2000 to simulate black holes in Bose- 
Einstein condensates 11561 . Shortly thereafter, the first explicit discussion of how 
BEC black holes could be used to simulate Hawking radiation in a laboratory is 
due to Barcelo, Liberati and Visser |X3f . A workshop held in Rio de Janeiro in 2000, 
exclusively dedicated to analogue models of gravity, marked the maturity of the 
field and led to the publication of the reference book [95]. 
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Finally, to end this brief overview, the 2005 review IITBI firmly settled the field 
of analogue gravity as an active and well-established research programme within 
the gravity community, while it seems that the idea of possible laboratory simula- 
tions of Hawking radiation has matured sufficiently to be taken very seriously in 
the condensed matter community, as illustrated for example by Il35ll43ll . 

Note that most of the work appearing on gravitational analogies can indeed be 
divided into the three categories that we discussed above: analogue gravity as such 
(with a main focus on laboratory simulations), quantum gravity phenomenology, 
and emergent gravity: 

• An increasingly important number of papers, in the line of [13] cited above, 
discuss the possibility of laboratory simulations of phenomena such as Hawk- 
ing radiation or particle creation in expanding universes. 

• A second set of works, following the central message of Unruh's original pa- 
per |[114| , use the gravitational analogy in fluid and other systems mainly as an 
inspiration for quantum gravity phenomenology, the lion's share being ded- 
icated to the possible influence of transplanckian physics on Hawking radia- 
tion. 

• Finally, Volovik |124| , apart from studying the possibility of laboratory experi- 
ments on the gravitational analogy in liquid Helium, is also the main advocate 
of taking the condensed matter analogy seriously as a candidate toy model for 
emergent gravity. 

1.5 Synopsis 

Most of the work discussed in this thesis is relevant mainly for the second level of 
interpretation, namely that of quantum gravity phenomenology. It also has obvious 
consequences for analogue gravity and laboratory simulations as such, in the fol- 
lowing sense. We will see that departures from the hydrodynamic regime in a BEC 
(at high energies and in the presence of horizons) predict deviations from the stan- 
dard relativistic phenomena related to black hole physics, such as Hawking radia- 
tion. Therefore, if these deviations are sufficiently strong, then apart from providing 
a possible window to quantum gravity phenomenology, they could also show up 
spontaneously and even inevitably in laboratory experiments reproducing the BEC 
analogues of these relativistic phenomena. We will end this thesis by making some 
comments on the third and strongest level of interpretation, which takes the con- 
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densed matter model seriously as a possible toy model for a full theory of quantum 
gravity 

The structure of this thesis is then the following. In chapter [2j we discuss some 
general preliminaries with respect to effective spacetimes in condensed matter sys- 
tems and in BECs in particular. Since this thesis is mainly directed at the (quantum) 
gravitational community we start with a brief refresher of some basic BEC theory 
including the Gross-Pitaevskii equation which describes the evolution of the con- 
densed phase of the BEC in the mean-field approximation. We then describe the ob- 
tention of the effective acoustic metric, and stress that this allows the consideration 
of black hole configurations. Finally we derive the Bogoliubov dispersion relation 
and look at the hydrodynamic limit. We emphasise that in this hydrodynamic limit, 
the dispersion relation acquires a relativistic form, and examine the validity of this 
hydrodynamic or relativistic limit. 

In chapter [3j we present a simple model for effectively (l+l)-dimensional black 
hole configurations in a BEC. The model consists of two homogeneous regions with 
a step-like discontinuity in the density profile. Although such a model is highly ide- 
alised, we argue that it captures the essence of the physics that we are interested in, 
at least for the purpose of studying the dynamical behaviour of black hole horizons 
in BECs. We derive the matching conditions that connect the magnitudes describ- 
ing the condensate, such as its density and phase and their derivatives, across the 
discontinuity. Finally, we develop a numerical method to find the dynamical modes 
(instabilities and quasinormal modes) of such configurations. 

This setup is then applied to the dynamical modes of such (l+l)-dimensional 
black hole configurations in a BEC, namely an analysis of the dynamical instabilities 
in chapter |4j and a quasinormal mode analysis in chapter |5j respectively. 

In chapter |4j we first indicate the slightly confusing state of the literature with 
respect to the dynamical stability of black hole horizons in BECs, and highlight the 
importance of this issue with respect to the possibility of simulating Hawking radia- 
tion in a BEC. We then thoroughly discuss the boundary conditions that are needed 
to apply the algorithm developed in the previous chapter, and analyse the possible 
appearance of dynamical instabilities. We discuss the results both for black hole 
and white hole horizons, and show that black hole horizons are generally stable, 
although we stress that modifying the boundary conditions (e.g., by introducing a 
wall or sink into the condensate) could affect this result. For a white hole horizon, 
the stability or instability is shown to depend even more crucially on the type of 
boundary conditions applied. In any case, we conclude that stable black hole con- 
figurations are possible, at least in principle. This means that studying the stable 
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dynamical modes or quasinormal modes, as well as the Hawking radiation, of BEC 
black hole configurations indeed makes sense. 

The structure of chapter|5]is the following. After a general introduction on quasi- 
normal modes of black holes (essentially: relaxation modes of the black hole after 
perturbation), we examine the boundary conditions that are applied in general rel- 
ativity to calculate the quasinormal modes of a black hole, and discuss how these 
should be modified in the case of BECs due to the superluminality of the disper- 
sion relation. We emphasise that in general relativity, no quasinormal modes exist 
in (l+l)-dimensional configurations, whereas a discrete spectrum appears in 3+1 di- 
mensions. The same is true for acoustic black holes in the hydrodynamic limit. We 
then apply the procedure outlined in chapter [3] to verify the occurrence of quasi- 
normal modes in the simple (l+l)-dimensional model presented earlier when using 
the full dispersion relation, including non-relativistic corrections. Not only do we 
find that such quasinormal modes indeed exist, but quite remarkably, their spec- 
trum turns out to cover a continuous region of the complex frequency plane. A 
rough dimensional estimate shows that these quasinormal modes are short-lived, 
indicative of a relaxation at the microscopic scale. This leads us to speculate that the 
discrete spectrum resulting from the standard quasinormal mode analysis in (3+1)- 
dimensional general relativity will also develop continuous bands if the dispersion 
relations at high energies undergo superluminal modifications. Such additional 
modes, though perhaps too short-lived to be detectable directly, could in princi- 
ple leave a trace of transplanckian (quantum gravitational) physics in the spectrum 
of gravitational waves emitted by black holes. 

In a brief intermezzo, chapter |6j we discuss in slightly more detail the link be- 
tween modified dispersion relations of the type that occur in BECs, and the Lorentz 
symmetry violation scenarios that are expected in quantum gravity phenomenol- 
ogy. This link justifies why we speculate that some of the results obtained in the two 
previous chapters can be directly extrapolated to such scenarios of quantum gravity 
phenomenology. It also allows us, instead of developing the gravitational analogy 
starting from the microscopic theory of a BEC as we did in the previous chapters, 
to directly study the impact of modified dispersion relations. Such an approach is 
then valid both for quantum gravity phenomenology in a scenario with superlumi- 
nal dispersion relations, and for analogue gravity in a BEC. We adopt this strategy 
in the next chapter. 

Chapter [7] studies the impact of superluminal dispersion relations on the late- 
time radiation or Hawking radiation from the formation of a black hole through a 
dynamical collapse. We start with a detailed motivation of why this study is rel- 



18 



CHAPTER 1. INTRODUCTION 



evant in spite of the extensive literature already existing with regard to Hawking 
radiation and the transplanckian problem. In particular, we highlight that most ex- 
isting studies have focused on stationary black hole configurations and subluminal 
dispersion. We also stress, as already mentioned at the beginning of this section, that 
if strong deviations from the standard Hawking picture are predicted, then these 
could spontaneously show up in laboratory simulations of the phononic equivalent 
of Hawking radiation in a BEC. We then discuss the classical geometry of our col- 
lapse model, and propose two concrete profiles, which differ in the behaviour of 
their surface gravity. The first profile considers a black hole with a constant surface 
gravity, while in the second one — characteristic of a Schwarzschild black hole — the 
surface gravity increases with frequency. This frequency-dependence of the surface 
gravity can be traced back to the fact that the horizon itself becomes a frequency- 
dependent concept because of the modified dispersion relation, an issue which we 
discuss in detail. The main body of the chapter consists in a concise derivation 
of the standard Hawking radiation, i.e., with relativistic dispersion relations, and 
a demonstration of how this derivation can be adapted to account for superlumi- 
nal dispersion relations. We obtain an analytic expression for the late-time radia- 
tion from the collapse of a black hole with superluminal high-frequency dispersion. 
Three important differences with the standard case can be deduced from this ex- 
pression. First, we highlight the dependence on a critical frequency which appears 
during the collapse process due to the superluminal dispersion, and above which 
no horizon is experienced. Second, we stress the importance of the surface gravity, 
which is shown to become frequency-dependent as well, again due to the super- 
luminality of the dispersion relation. Finally, we point out that the radiation orig- 
inating from the collapse process dies out as time advances. We then numerically 
integrate the analytic expression for the late-time radiation, and graphically illus- 
trate the three issues just mentioned. In particular, it is clearly seen that the exis- 
tence of a critical frequency implies that the radiation is fainter than in the standard 
case, even in a profile with constant surface gravity, unless this critical frequency 
goes to infinity (i.e., unless a singularity is formed at the end of the collapse). We 
also illustrate that the dependence on the surface gravity could lead to important 
ultraviolet contributions, depending on the relation between the critical frequency 
and the Lorentz symmetry violation scale. Finally, we show the time evolution of 
the radiation originating from the collapse and indicate that, for a natural choice of 
magnitudes, this extinguishes extremely rapidly. We conclude with a comparison 
of our results with other recent results available in the literature, and speculate that 
some of the standard assumptions that are taken when arguing for the robustness of 
Hawking radiation, which are directly extrapolated from the standard case, might 
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be model-dependent and actually broken in a physically realistic collapse scenario 
with superluminal dispersion. 

Chapter[8]discusses some aspects of how the analogue gravity framework in con- 
densed matter models can be extended to a serious model for emergent gravity. We 
start with some general considerations on what quantum gravity could or should 
be, and mention that one of the main motivations for an emergent gravity approach 
is the dark energy problem. The bulk of this chapter is devoted to the question 
of diffeomorphism invariance in emergent gravity models based on the condensed 
matter analogy. We argue that this question should be split up into a kinematical 
aspect and a dynamical aspect. We discuss why the kinematical aspect of diffeomor- 
phism invariance poses no real problem, since diffeomorphism invariance emerges 
at the same level as Lorentz invariance, i.e., it is realised as an effective low-energy 
symmetry. This discussion involves some subtle concepts concerning the relation 
between internal observers and microscopic physics. We end with some remarks 
with respect to the dynamical aspect of diffeomorphism invariance. We argue that 
this is essentially identical to the problem of recovering the Einstein equations in 
analogue gravity models, and discuss the possibility of using Sakharov's induced 
gravity to such an end, an issue which is under current investigation. 

We round off in chapter [9] with an extensive general summary and outlook. 
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Chapter 2 
Preliminaries 



2.1 General considerations 

It is widely expected that underneath the general relativistic description of gravita- 
tional phenomena there is a deeper layer in which quantum physics plays an impor- 
tant role. However, at this stage we do not have sufficient intertwined theoretical 
and observational knowledge to formulate an appropriate description of what un- 
derlies gravity. Moreover, starting from structurally complete quantum theories of 
gravity, it could still be very difficult to extract the specific way in which the first 
'quantum' modifications to classical general relativity might show up. This hap- 
pens for example within the loop quantum gravity approach Il6l ll05ll , while in string 
theory a huge 'landscape' of quantum vacua seems to arise ||112 | 1 , although it is not 
clear yet which ones — if any — could correspond to our universe ||78| . 

Analogue models of general relativity in condensed matter systems provide spe- 
cific and clear examples in which effective spacetime structures ultimately emerge 
from (non-relativistic) quantum many-body systems. As mentioned in the intro- 
duction, the key point is the realisation that some collective properties of these con- 
densed matter systems satisfy equations of motion formally equivalent to those of 
a relativistic field in a curved spacetime. The simplest example is the equation de- 
scribing a massless scalar field (p, 



satisfied, for example, by acoustic perturbations in a moving perfect fluid. These 
acoustic perturbations therefore travel along the null geodesies of the acoustic met- 
ric g, lv . The finite speed of propagation of the perturbations in the effective fields 




(2.1.1) 
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(in our case, the speed of sound) plays the role of the speed of light in fundamental 
physics. 

This example illustrates how for certain (semiclassical) configurations and low 
levels of resolution, one can appropriately describe the physical behaviour of the 
system by means of a (classical or quantum) field theory in a curved Lorentzian 
background geometry However, when one probes the system with higher and 
higher resolution, the geometrical structure progressively dissolves into a purely 
microscopic regime [9j. Therefore, although analogue models cannot be considered 
at this stage complete models of quantum gravity (since they do not reproduce the 
Einstein equations), they provide specific and tractable models that reproduce many 
aspects of the overall scenario expected in the realm of real gravity. 

The main objective of studies such as the ones that we will present in this thesis 
is then to obtain specific indications about the type of deviations from the general 
relativistic behaviour to be expected when quantum gravitational effects become 
important. All, under the assumption that the underlying structure to general rela- 
tivity is somewhat similar to that in condensed matter systems. For reasons that we 
will elaborate in this chapter, we are interested in the case of Bose-Einstein conden- 
sates (BECs). 

There are two ways of proceeding. The first way is to start from a theoretical 
description of BECs, elaborate the gravitational analogy and see which type of de- 
viations arise with respect to the relativistic model. By extrapolation, one can then 
expect the same kind of deviations to occur in gravity if quantum corrections to 
gravity have a nature similar to the non-relativistic corrections that show up in BECs 
beyond the hydrodynamic or relativistic limit. 

The second way is to directly calculate some possible quantum corrections in 
gravity, by modifying the standard relativistic description of a (predicted) gravita- 
tional phenomenon in a way which immediately incorporates the kind of modifi- 
cations suggested by the condensed matter analogy. This essentially consists in re- 
placing the relativistic dispersion relations by superluminally modified dispersion 
relations in a consistent way, as we will discuss in detail further on. 

We will apply the first approach to an analysis of the dynamical stability (chap- 
ter [4]) and of the quasinormal mode spectrum (chapter |5j) of black hole configurations 
in Bose-Einstein condensates, based on a (l+l)-dimensional model which we intro- 
duce in chapter |3j Subsequently, the second approach will be applied to Hawking 
radiation for collapsing configurations (chapter [7|). First of all, though, in this chap- 
ter, we will review some basic elements of BEC theory and of how the gravitational 
analogy in a BEC is obtained. 
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2.2 Bose-Einstein condensates of dilute atomic gases 

In this section, a brief general refresher of the theory of Bose-Einstein condensation 
will be given fl36l |45l together with some relevant numbers related to their experi- 
mental achievement [ 17211111] . 



A Bose-Einstein condensate is an aggregation state of matter which arises in 
a many-body system of bosons, confined through an external potential, when a 
macroscopic amount of them are led to occupy the ground state determined by this 
external potential. Such bosonic many-body systems can be described in terms of 
quantum field operators Y(r) and Y + (r) that annihilate and create particles, together 
with the appropriate equal-time commutation rules 

[f (r),f CO] = , [Y + (r), Y + (r')] = , [f (r), Y + (r')] = S(t - r') , (2.2.1) 
and a Hamiltonian typically of the form 

H = J drY + (r)k (r)Y(r) + ^ J drdr'f + (r)f + (r') V(i' - r)Y(r)Y(r') . (2.2.2) 
Here, 

Mr) = -^V 2 + Vext(r) (2.2.3) 
Zm 

is the non-interacting single-particle Hamiltonian, m the atomic mass, V ext the exter- 
nal trapping potential, and it has been assumed that the many-body interactions are 
well modelled by a two-body interaction potential V(r' — r). 

When the system is cooled below a certain critical temperature T cr i t , it condenses]^] 
and Y can be separated into a macroscopic wave function xp plus quantum fluctu- 
ations. In atomic systems, the interaction between the atoms implies that some of 
them are excited even at zero temperature, in contrast to the ideal or pure Bose con- 



densate depicted in Fig. [Zl The proportion of excited atoms or quantum depletion 
increases from typically 0.1% at zero temperature to still less than 1% at the non- 
zero temperatures achieved in laboratory experiments. By contrast, in helium, the 
main other quantum liquid used for the study of gravitational analogies, the quan- 
tum depletion is of the order of 90%, making the theoretical treatment much more 
complicated. 

Note that the theory of Bose-Einstein condensation is in principle applicable to 
bosonic many-body systems in general. However, most systems would experience 



lr rhe condensation of a BEC takes place in momentum space: all atoms tend to the lowest energy 
state characterised by k = 0. For experimental implementations with atomic gases, the condensation 



in momentum space is accompanied by a condensation in position space, see Fig. 2.2 due to the use 
of harmonic traps. 
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a classical phase transition to liquid and solid state long before the quantum de- 
generacy occurs, thereby seriously complicating the simple theoretical treatment of 
the condensation described so far. In order to avoid this, the most successful ap- 
proach has been to use weakly interacting alkaline atomic gases with extremely low 
densities, i.e. dilute atomic gases, combined with rapid cooling techniques to prevent 
solidification. 

For such gases, the evolution equation for the condensed part xp can be obtained 
as follows. Binary collisions characterised by the s-wave scattering length a are the 
main relevant interaction in a cold dilute gas. One can then replace the interaction 
potential V(r' — r) in the Hamiltonian ( 2.2.2| > by the effective prescription 



V(r' - r) = gS(t' - r) , (2.2.4) 
with g = 4nh a/m the effective coupling constant. The Hamiltonian becomes 

H = J d 3 r [Y%Y + |Y + Y f YY] , (2.2.5) 

and the evolution equation for Y is 

ihd t Y = [Y, H] = (h Q + gY + Y)Y . (2.2.6) 

Since the quantum depletion is small, it makes sense to use a mean-field approach 
and expand the theory as a perturbation series around the macroscopic wave func- 
tion xp. The zeroth-order term can be obtained simply by substituting xp for Y in 
eq. ( 2.2.6| >, which gives the non-linear Schrodinger equation or Gross-Pitaevskii equa- 



tion 

ihdtf = (h + g\xp\ 2 )xp . (2.2.7) 

This classical mean-field equation is known to provide an excellent approximation 
for the evolution of the condensed phase xp in the case of weakly interacting systems 
(s-wave scattering length a much smaller than the interatomic distance, or in other 
words: na 3 <C 1, with n the particle density) containing a large number of atoms. 

Bose-Einstein condensation was first experimentally observed in 1995 at JILA 



in rubidium ( 87 Rb), and at MIT in sodium ( 23 Na), see Fig. 2.2 Since then, about 
a dozen other atomic gases including 7 Li, 1 H and 85 Rb have followed, 87 Rb and 
23 Na still being the favourites from an experimental point of view Apart from 

atomic gases in weakly interacting regimes, Bose-Einstein condensation has also 
been achieved in strongly interacting regimes [44J by tuning the scattering length 
through the use of a magnetic Feshbach resonance 11521 ISilig This technique was 

2 Note that condensation of bosons in a strongly interacting regime is also responsible for the 
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Wh at is Bose-Einstein condensation (BEC)? 





High 
Temperature T: 

thermal velocity v 
density d" 3 
"Billiard balls" 




Low 
Temperature T: 

De Broglie wavelength 
^B=h/mv « T ia 
"Wave packets" 




T=T cr jt: 
Bose-Einstein 
Condensation 

>WJB»d 

"Matter wave overlap" 

T=0: 
Pure Bose 
condensate 

"Giant matter wave" 





Figure 2.1: Criterion for Bose-Einstein condensation, taken from |72|. d is the inter- 
atomic distance and A^g the thermal or de Broglie wavelength, which is a measure 
for the extent of the wavepackets associated with the atoms (A^g oc T~ 1/2 ). 



also used to suddenly switch from a weak repulsive to a strong attractive inter- 
action, producing a violent collapse followed by a partial explosion process called 
"Bose nova" because of its similitude with astrophysical supernova explosions [47|. 
As another striking example, in magnons, a collective excitation of the atomic spin 
structure of a crystal in the presence of an external magnetic field, Bose-Einstein 
condensation has recently been achieved at room temperature [46 J. 

A few illustrative numbers associated with BECs in atomic gases are the follow- 
ing |5l The critical temperatures involved are typically from a few hundred 
nK to a few ]iK, at densities between 10 14 and 10 15 atoms/ cm 3 . The total amount of 
atoms in a BEC ranges from only a few hundred to roughly 100 million in Na and 
Rb and about 1 billion in atomic hydrogen (note that larger amounts are desirable 
mainly because they make the experiment more robust with respect to impurities 
and misalignments), with typical condensate sizes varying from 10 }im to a few hun- 
dred ^m. Finally, two parameters that we will regularly need in this thesis are the 
speed of sound c and the healing length £, which, roughly speaking, is the distance 
needed for the condensate to smoothen out a sharp inhomogeneity in the atomic 
density. Typical values for the velocity of sound c in BECs range between 0.5mm/ s 
and 10mm/ s, while the healing length £ usually lies between lOOnm and 2^m. 



superfluidity of 4 He, a phenomenon which has been observed as early as in 1938 [71J. However, 
4 He is not well described by the standard mean-field theory of Bose-Einstein condensation since less 
than 10% of the atoms are in the condensed state, due to the strong interaction between the atoms. 
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Figure 2.2: Experimental observation of Bose-Einstein Condensation in 87 Rb at JILA 
(above) and in 23 Na at MIT (below). In both pictures, the left frame represents the gas 
at a temperature just above the critical condensation temperature, the centre frame 
the same gas just after the condensation and the right frame a nearly pure conden- 
sate. The upper picture is a computer image of the momentum distribution, the 
lower one of the spatial distribution. 

Apart from being a very active field of research in themselves, Bose-Einstein 
condensates have also become an "ultralow-temperature laboratory for atom op- 
tics, collisional physics and many-body physics" [5J. And one might add: "for high- 
energy physics", as illustrated not only by the many works on analogue gravity in 
BECs, but also for example on braneworlds in warped spacetimes fHf and a sugges- 
tion to probe certain aspects of supersymmetry in a Bose-Einstein condensate 111101 . 

Here, we will focus on the gravitational analogy in BECs in regimes that can be 
described by mean-field theory. 



2.3 Gross-Pitaevskii equation and quantum potential 

In the previous section, we arrived at the Gross-Pitaevskii equation 




(2.3.1) 
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which describes the evolution of the condensate's macroscopic wave function or 
order parameter ip, and where we have now explicitly indicated the spatial and 
temporal dependencies]^] We remind the reader that m is the boson mass, V ex t the 
external potential and g a coupling constant which is related to the corresponding 
s-wave scattering length a through g = 4nh a/m. 

The Gross-Pitaevskii equation can be expressed in terms of hydrodynamic quan- 
tities such as the local speed of sound c and the velocity of the fluid flow v. This 
proceeds by first introducing the hydrodynamic or Madelung representation for the 
order parameter in the Gross-Pitaevskii equation: 

ip = ^ne i9/h e-^ t/n . (2.3.2) 

Here n is the condensate density ji the chemical potential and 8 a phase factor. Sub- 
stituting in eq. ( 2.3.1| > we arrive at 



d t n = _lv-(nV0), (2.3.3) 
m 



d t d = -^Z( VQ ) 2 ~S n ~ Vext-ji- Vquantum, (2.3.4) 



1 

2m 

where the so-called 'quantum potential' is defined as 



h 2 V 2 ^/n 

^quantum — ^ p= • (2.3.5) 

2m yjn 



In many situations the quantum potential in eq. ( |2.3.4[ ) can be neglected, a point 



to which we will return in section 2.6 The resulting equations ( 2.3.3| > and ( 2.3.4| > 



are then equivalent to the hydrodynamic equations, i.e., a continuity equation and 
an Euler equation, for a perfect classical fluid. In this case, it is well known that 
the propagation of acoustic waves in the system can be described by means of an 
effective metric, thus providing the analogy with the propagation of fields in curved 
spacetimes ||114t I117| , as we will illustrate in the next section. Given a background 
configuration (hq and 6q), this metric can be written as 




c 2 



(2.3.6) 



where 1 is the 3x3 unit matrix, c 2 = gno/m and v = 'VOq/m. These magnitudes, 
c and v, represent the local speed of sound and the local velocity of the fluid flow 
respectively. 



3 Note that the external potential V ex t can also be time-dependent, and can even be used to tune 
the coupling constant g through a Feshbach resonance, as mentioned in the previous section. For 
simplicity, we will assume both V^t and g time-independent. 
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2.4 Effective acoustic metric 



Acoustic waves or phonons are linear perturbations of a background solution to the 



Gross-Pitaevskii equation. In order to linearise the eqs. (2.3 1, we write 

n(r,f) = n Q {x)+g- l n 1 {x,t) , (2.4.1) 

0(r,f) = o (r)+0 1 (r,f), (2-4.2) 

where n\ and 6\ are small perturbations of the density and phase of the BEC (the 
coupling constant g is included in the linearisation of n for notational convenience). 



The eqs. ( |2.3| > then separate into two time-independent equations for the background: 








-V • (c 2 v) 



1 



-rav 



mc 



Vc 



n • ext r 1 r, 

2 2m c 

plus two time-dependent equations for the perturbations: 

dfii = -V • (niv+ c 2 V0i) , 
d t di = -v • V0i 



]i + 



ft 2 v 2 c 



ni + ^ 2 V 



c 2 V 



ni 

r2 



(2.4.3a) 
(2.4.3b) 

(2.4.4a) 
(2.4.4b) 



ft/ (mc) is the healing length of the condensate, which, as we mentioned 
2.21 can be understood as the distance over which an otherwise homogeneous 



where £ 
in 



condensate recovers or 'heals' from a sharp inhomogeneity. The effective acoustic 
metric ( 2.3.6| > can then be derived as follows. 

In the hydrodynamic or acoustic limit, the last term on the right-hand side of 



eq. ( |2.4.4b[ ), which is a direct consequence of the presence of the quantum poten- 
tial (2.3.5) in eq. ( 2.3.4| >, can be neglected. Then, extracting n\ from the remaining 



terms in eq. (2.4.4b I and plugging into eq. (2.4.4a), the result can be written in four- 
dimensional notatiorEJas 

where 

n ( -1 



0, 



(2.4.5) 




Define y/—gg }ll/ = p v , where g = det(gL v ), and hence also g = det(/ 
ing gl lv , we obtain the acoustic metric 



flV) 




(2.4.6) 
By invert- 

(2.4.7) 



4 We use the convention x ?< = (f, r). 
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Note that we have derived this acoustic metric from the Gross-Pitaevskii equation. 
However, the same result can be obtained directly by combining and manipulating 
the Euler and continuity equations for an irrotational flow of any non-viscid and 
barotropic fluid ||117| . Some indications of why we are nevertheless specifically in- 
terested in BECs have already been given in the introduction. We will postpone a 



more detailed motivation until section 2.7 



From eq. ( 12.4.71 ), the wave equation for the propagation of the fluctuation field Q\ 



of the velocity potential (the sound wave), is 

ne 1 = (TV^Vvfli = — La^v^^A) = o . (2.4.8) 

V ~~ 9 

In other words, acoustic perturbations propagate according to a d' Alembertian equa- 
tion of movement, which is formally identical to the wave equation of a massless 
scalar field propagating in a curved spacetime. 

Before moving back to the main discussion, a few remarks might be in place. 
First, note that the effective or acoustic metric has a Lorentzian signature (—,+,+,+) 
So the collective behaviour of the condensate really gives rise to the emergence of 
an effective curved background spacetime for the propagation of the phonons. 

Second, it is important to stress that we are dealing with a bi-metric system. Of 
course, each part of the system (background and perturbations) can only couple to 
one metric. The first one is the usual Minkowski metric for flat spacetime 

r\w = (diagh4 ght , 1, 1, 1] ) V v • (2.4.9) 

This metric governs the background fluid, in which relativistic effects are totally 
negligible, since for the fluid flow velocities involved: v <C Cii g ht- The background 
fluid couples only to this metric. The phonons or acoustic perturbations, on the 
other hand, do not see this flat-space metric but are instead governed only by the 
effective acoustic metric g„ v . It is thus seen that the classical Newtonian dynamics of 
the background fluid leads to an effective low-energy relativistic curved Lorentzian 
spacetime metric for the sound waves. 

Third, the sign of gtt in ( 2A.7) depends on the values of c and v, i.e., on the fluid 



regime. A change of sign in gtt characterises an ergoregion. For a purely radial 
flow, ergoregion and horizon coincide. Therefore, if we can manipulate our system 
in such a way that for some r^, a subsonic outside region and a supersonic inside 
region are created, i.e.: 

c 2 (r) > v 2 (r) for r > (subsonic outside region) , 
c 2 (r) < v 2 (r) for r < (supersonic inside region) , 



30 



CHAPTER 2. PRELIMINARIES 



then there will be a horizon at r^: 

c(r h ) = \v(r h )\ . (2.4.10) 

Indeed, for c 2 — v 2 > 0, gtt is negative, indicating the timelike character of the time 
coordinate, whereas for c 2 — v 2 < 0, gtt becomes positive and so this coordinate 
becomes spacelike. As time marches on, any particle in the supersonic inside region 
will irrevocably move towards the centre of the system (r = 0), where the singularity 
would lie in the case of a gravitational black hole. 

Finally, in general we will be mainly interested in fluctuations on top of a sta- 
tionary background, with the notable exception of chapter [7} where we examine the 
Hawking radiation in a collapsing geometry. In stationary cases, we can write 

d t n = d t c = d t v = 0. (2.4.11) 

As a further simplification, we will focus on effectively (l+l)-dimensional config- 
urations and write c(t, x) and v(f, x). This means that we consider perturbations 
propagating in a condensate in such a way that the transverse degrees of freedom 
are effectively frozen. In other words, the only allowed motions of both the pertur- 
bations and the condensate itself are along the x-axis. 

We now move back to the main line of discussion. 



2.5 Dispersion relation 



A crucial element in the quantum gravitational analogy is the phononic dispersion 
relation in a BEC and how it compares to its relativistic counterpart. 

To derive the dispersion relation, consider a region in which the condensate is 



homogeneous (with c and v constant), and seek for solutions of eqs. (2.4.4) in the 
form of plane waves: 

ni(f,x) = Ae^-^, (2.5.1a) 
e x (t,x) = Bjfa-o*) f (2.5.1b) 

where A and B are constant amplitude factors. We stress that the frequency co and 
the wave number k are not restricted to real values, but are allowed to be complex. 
The imaginary parts of co and k will encode the dynamical behaviour of the system. 
Substituting into eqs. ( |2.4.4j ) we find 

(i{co-vk) c 2 k 2 \ /A N 



\l + lZ 2 k 2 -i(co-vk)J 



0. 



(2.5.2) 



B 
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For a non-trivial solution to exist, the determinant of the above matrix must vanish. 
This immediately leads to the Bogoliubov dispersion relation. 

(cv - vkf = c 2 k 2 + ^c 2 £ 2 fc 4 . (2.5.3) 

This equation is globally valid in homogeneous condensates. In non-homogeneous 
configurations, it is still valid as a strictly local dispersion relation, provided that the 
quantities involved (c, v, £ and k) all vary slowly with x. 



The full dispersion relation (2.5.3) is quartic in k. It should be compared to the 



hydrodynamic dispersion relation obtained by neglecting the quantum potential: 

(cv - vk) 2 = c 2 k 2 . (2.5.4) 

Such a hydrodynamic dispersion relation is of the usual quadratic relativistic form. 
This fundamental observation allows us to refer to the hydrodynamic or the relativistic 
limit interchangeably. 

However, we are not so much interested in this hydrodynamic or relativistic 
limit, but precisely in deviations from it. Therefore we will always work with the full 
dispersion relation, and only use the hydrodynamic limit as a means of comparison. 
The key point is that, because the full dispersion relation is quartic, when writing 
the mode u w corresponding to a certain frequency co: 

u lv = e-^ t Y J A j e lk i x , (2.5.5) 

;' 

there are now four contributions to this mode, stemming from the four values of k 
associated with each value of a;. In the hydrodynamic limit there are obviously only 
two contributions to each u^-mode: 

fc 12 = cv/(v±c), (2.5.6) 

In the configurations that we will analyse, the additional modes will turn out to play 
a crucial role at all frequencies because of the presence of black hole horizons, as we 
will discuss in the next section. 

Another crucial aspect of the non-hydrodynamic corrections at high frequencies 
is that these make the dispersion relation 'superluminal' (i.e., superphononic or su- 



personic). Indeed, the group velocity for the modes which behave like eq. (2.5.6) in 
the hydrodynamic limit can in general be written as 

Vg = -^ = v±c k , (2.5.7) 
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Figure 2.3: Behaviour of the effective speed of sound C\ as a function of the wave 
number A; in a BEC. Due to high-frequency non-hydrodynamic (non-relativistic) 
corrections, the effective speed of sound starts to increase significantly around the 
Lorentz symmetry breaking scale ki = 2/£, and becomes effectively superphononic 

(c k > c). 



with an effective propagation speed 



1 + \^ 2 



c k = c ^ , (2.5.8) 

; i + k 2 k 2 



(see Fig. |2.3[ >. For sufficiently large wave numbers, becomes significantly larger 
than c, i.e., it becomes effectively superphononic. 

The connection with quantum gravity phenomenology is the following. In a 
wide range of phenomenological scenarios for quantum gravity, violations of the 
local Lorentz symmetry are expected at high frequencies. When these Lorentz viola- 
tions are associated with the existence of privileged reference frames, they typically 
lead to modified dispersion relations of the kind 

rA 

oo 2 = c 2 k 2 ±c 2 -j , (2.5.9) 

with ki the local Lorentz symmetry breaking scale (and obviously c represents the 
speed of light in this last equation). BECs therefore provide a concrete example of 
such a dispersion relation with 'superluminal' Lorentz-breaking corrections, charac- 
terised by a Lorentz violation scale ki = 2/£. We will come back to the connection 
with quantum gravity phenomenology in chapter [6j 



2.6 Validity of the hydrodynamic approximation 

The functions c(t,x) and \(t,x) completely characterise the acoustic metric. In gen- 
eral relativity, any metric has to be obtained by solving the Einstein equations. Here, 
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however, the magnitudes c(t, x) and v(f, x), and so the acoustic metric, are those sat- 



isfying the continuity and Euler equations of hydrodynamics, eqs. (2.3 1 without the 
quantum potential. Thus, these equations play a role analogous to the vacuum Ein- 
stein equations in general relativity. Of course, at the global non-linear level these 
hydrodynamic equations are completely different from the real Einstein equations. 
But we are interested in their behaviour when linearised around a background so- 
lution. Although this does not lead to the full complexity that is expected in general 
relativity (if only because the fluctuations of the metric here are characterised by a 
scalar or spin-0 field, versus a spin-2 tensor field in general relativity), it still cap- 
tures much of the essence of a proper linearised general relativistic behaviour. 



There exist, however, situations in which the quantum potential in eq. (2.3) can- 
not be neglected, and the hydrodynamic approximation is not valid. To analyse 
this, we start from the full dispersion relation. One can obtain the hydrodynamic 
approximation by writing a; as a series expansion in terms of gk, and truncating the 
higher order terms, i.e.: 



co = \v±c^l + ^k 2 jk 

~ (v±c)kT \ctfk 3 + 0{^k 5 ) ~ (v±c)k+0(£ 2 k 3 ). (2.6.1) 
8 

For the series expansion to make sense, one must have £k <C 1, i.e., the frequencies 
must be small compared to the characteristic scale of the system, or in other words 
the wavelengths must be much larger than the healing length £ = %/ (mc). The 
truncation moreover requires 

? 2 k 2 

* <1. (2.6.2) 



8\l±v/c\ 



This second condition is never satisfied around a horizon (c = \v\). Therefore, near a 
horizon, the additional modes that are not present in the hydrodynamical limit are 
always important, even at low frequencies. Whenever probing the configuration 
either with high frequencies, or near the horizon (with any frequency), the geomet- 
rical picture will dissolve and the underlying microscopic structure will appear to 
some degree. 



Without forgetting this subtlety, we will continue to call eq. ( |2.3.6| > the 'effective' 
or 'acoustic' metric in the system, even when analysing the full Gross-Pitaevskii 
equation. Then, we can consider the equations ( |2.3| > to play a role similar to the vac- 



uum Einstein equations. Their treatment is classical, but they incorporate correc- 
tions containing h. Therefore, the standard treatment of BECs based on the Gross- 
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Pitaevskii equation provides an example of a way of incorporating quantum cor- 
rections to the dynamics of a system without recurring to the standard procedures 
of back-reaction. Again, at the global non-linear level these equations bear no rela- 
tionship whatsoever with any sort of 'semiclassical' Einstein equations, nor do they 
present sufficient complexity to lead to spin-2 fields that could be the analogue of 
gravitational waves. However, at the linear level, that is, in terms of linear tenden- 
cies of departing from a given configuration, equations (2.3) encode much of the 
essence of the linearised general relativity behaviour (a Lorentzian wave equation 
in a curved background), semiclassically modified to incorporate a superluminal 
dispersion term, as we have already discussed. 



2.7 Why Bose-Einstein condensates? 



We mentioned right after eq. (2.4.7) that the effective metric could have been derived 
for any irrotational, non-viscid and barotropic fluid, real or idealised. There are 
nevertheless several good reasons why we are interested specifically in BECs. 

First of all, real fluids have the obvious advantage over idealised fluids of per- 
mitting real experiments. It is not too difficult to simulate a black hole horizon in 
a real fluid system. For example, engineers routinely produce acoustic horizons in 
wind tunnels. Actually, at first sight, any moving fluid, including plain water, might 
serve as the background fluid in which to conduct the gravitational analogy, at least 
at the level of creating a supersonic flow and associated acoustic horizon. However, 
water and other 'normal' fluids show a high degree of viscosity which introduces 
terms that cannot be moulded in the relativistic form described earlier, and hence 
break the gravitational analogy even at low energy. Furthermore, many quantum ef- 
fects of interest, such as Hawking radiation, would become indistinguishable from 
the background noise and impurities in normal fluids such as water. Finally, es- 
timates for the Hawking temperature in water indicate an order of magnitude of 
1 fiK, which is clearly hopeless to detect at temperatures for which water behaves as 
a fluid. Because of these problems, the quest for a useful background fluid leads in 
the direction of superfluids. Superfluids have a vanishing viscosity, and some super- 
fluids can be made extremely pure. Two main candidates stand out: Bose-Einstein 
condensates in dilute atomic gases and superfluid 3 He. 

BECs have the comparative advantage of being conceptually well understood, 
and relatively simple to describe theoretically and manipulate experimentally. This 
is partly due to the extremely low densities usually involved, which means that the 
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interatomic distance is typically a lot greater than the scattering length]^] leading to 
rather weak interactions. In superfluid 3 He systems, the density is several orders 
of magnitude larger than in the BEC, leading to an interatomic separation of the 
same order of magnitude as the scattering length. As a consequence, only a frac- 
tion of all the particles are inside the condensate, while the others form a second, 
normally fluid phase. Furthermore, the fundamental degrees of freedom in 3 He are 
fermionic. The superfluidity arises from Cooper-pair formation between weakly 
bound fermions. These pairs behave as bosons, but some of the properties of the 
individual fermions persist. All of this stimulates the existence of many secondary 
effects, which make 3 He a complex and rich model, in many ways much more ap- 
propriate for thinking in analogue terms about our universe than the simple model 
provided by BECs [124]. But in terms of conceptual simplicity, Bose-Einstein con- 
densates clearly stand out as the best option. As a simple example, in Bose-Einstein 
condensates under normal energetic circumstances, atomic interactions only take 
place in the s-wave. A much higher temperature than the condensation temper- 
ature, or a far higher atom density, would be required to excite higher multipole 
wavesj^] In superfluid 3 He on the other hand, the size of the atom, the thermal de 
Broglie wavelength and the interatomic separation are all of the same order of mag- 
nitude. This creates a complex and rich situation, but also blurs and mixes effects 
which might be interesting to study separately. 

Finally, the third reason why we are interested in BECs, is that they provide a 
real system in which modified dispersion relations arise. As we mentioned in the 
introduction, and will further detail in chapter [6j they therefore offer an interesting 
model for a prominent line of research within quantum gravity phenomenology. 
The mere occurrence of Lorentz-breaking modifications of the dispersion relations 
shows that there exist concrete examples in nature in which Lorentz invariance is 
realised as a low-energy effective symmetry, thereby strengthening the assumption 
that this might be the case for the local Lorentz invariance of general relativity as 
well. Furthermore, such modified dispersion relations also model most types of 
possible Lorentz symmetry violations at high energies that are expected from a the- 
oretical quantum gravity phenomenological point of view (at least the ones associ- 
ated with a preferred reference frame). Also, in contrast to sub luminal dispersion, 
which gradually dampens the influence of high frequencies, superluminal disper- 



5 Typical values for 23 Na at a density of approximately 10 14 atoms / cm 3 are 3 nm for the scattering 
length, versus 200 nm for the interatomic distance f72\ . 

6 For the same example of 23 Na at a density of approximately 10 14 atoms/cm 3 , the limiting tem- 
perature below which only s-wave scattering occurs is ~ lmK, while the condensation temperature 
is ~ 2}iK 121. 
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sion relations such as those that arise in BECs amplifies the effects of high-frequency 
('transplanckian') physics. This makes BECs an interesting probe for the possible in- 
fluence of microscopic corrections to classical gravity, as well as for the robustness 
of phenomena related to black holes, such as Hawking radiation, see chapter U\ 



Chapter 3 



Black hole configurations in 
Bose-Einstein condensates 



3.1 A simple (l+l)-dimensional model 



We will now apply the theory developed in the previous chapter to analyse the dy- 
namical behaviour of effectively (l+l)-dimensional BEC configurations with density 
and velocity profiles containing a single acoustic black hole horizon, see Figs. 3.1 and 



3.2 This analysis will be separated into two parts: the dynamical stability will be 



examined in chapter |4| and a study of the quasinormal modes is contained in chap- 
ter |5j First, in the present chapter, we will discuss the background configurations 
and numerical method used in these analyses. 

The type of (l+l)-dimensional background profiles that we will examine con- 
sists of two regions each with a uniform density and velocity, connected through 
a steplike discontinuity, see Fig. 3.1 This is obviously an idealised case, obtained 
from Fig. 3.2 in the limit when the intermediate transition region L — > 0. We always 
consider flows moving from right to left, hence a negative value of v. 

The parameters of the configuration are chosen such that there is an acoustic 
horizon at the x = discontinuity, which is the formal equivalent of a general rel- 
ativistic black hole horizon, as described in section 2.4 The left-hand asymptotic 
region (x — ► — oo) mimics a general relativistic black hole singularity in the sense 
that nothing can escape from it towards the horizon (we will always impose an ade- 
quate boundary condition in order to ensure this). Note, however, that in our model, 
this region is perfectly regular, so the analogy with a general relativistic singularity 
is limited to the previous aspect. 
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Figure 3.1: Flow and sound velocity profile consisting of two homogeneous regions 
with a steplike discontinuity at x = simulating a black hole configuration in a BEC. 
The solid blue line represents the speed of sound c, the dashed red line the fluid 
velocity v. In the upper part of the picture, the negative value of v indicates that the 
fluid is moving leftwards. For x > 0, the fluid is subsonic since c > \v\. At x < 
it has become supersonic. At x = 0, there is a sonic horizon. This picture can be 
considered the idealised limiting case of Fig. 3.2 when L — > 0. 



At each discontinuity, matching conditions apply that connect the magnitudes 
describing the condensate at both sides of the discontinuity. Furthermore, we need 
a set of boundary conditions, which determine what happens at the far ends of 
the condensate. All these elements determine the characteristics of the system, and 
hence its eigenfrequencies. 

In the remainder of this short chapter, we will describe this particular model of 
(l+l)-dimensional configurations with abrupt discontinuities. 



3.2 Background configuration and numerical method 

The profiles that we will study are piecewise homogeneous with a single discon- 
tinuity. In order to find their dynamical modes, we will first of all need a way to 
describe the solutions that characterise the profiles under study. The idealised form 
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Figure 3.2: Flow and sound velocity profile consisting of two homogeneous regions 
connected through a transition region of width L around the sonic horizon at x = 0. 
This profile simulates a black-hole-configuration in a BEC. As in Fig. 3.1 the solid 
blue line represents the speed of sound c and the dashed red line the fluid velocity v. 
Actual calculations presented in this work will be based on the limiting case L — > 0, 
see Fig. |3T 



of these profiles allows a simplified description in each homogeneous region. Then, 
the various sections have to be linked through matching conditions at the discon- 
tinuity. After finding all the mode solutions, the dynamical modes (instabilities or 
quasinormal modes, respectively) will be those that moreover satisfy the adequate 
boundary conditions. We will describe the matching conditions in the current sec- 
tion, together with a straightforward algorithm for numerical implementation. The 
boundary conditions depend on the concrete problem, the calculation of instabili- 
ties or of quasinormal modes, and will therefore be described in chapter |4] and |5j 
respectively. 

In order to calculate the matching conditions that link the description of the con- 
densate in two homogeneous regions around a discontinuity, we integrate the differ- 
ential equations ( |2.4.4 > that describe the evolution of the perturbations around this 
discontinuity (say x = 0). This results in the following four independent, generally 
valid, matching conditions: 

[0l] =0, [wi + c 2 d x 0i 



0, 



c 2 d x 



ill 



= 0, 
0. 



(3.2.1a) 
(3.2.1b) 



The square brackets in these expressions denote, for instance, [0j] 



0l| x =O+ ~~ #lL=0- 
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We can simplify the second condition in eq. ( |3.2.1a| > to [c 2 9 x #i] = by noting that 
[vn\] = because of the background continuity equation ( 12.4.3a! ), while for our 
choice of a homogeneous background the last condition becomes simply = 0. 

For a given frequency co, the general solution of eqs. ( 2.4.4[ > can be written as 



"l 



£ Af^i*-^) (x < 0), 
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[x < 0), 

(3.2.2) 

[x > 0), 



where {kj} are the roots of the corresponding dispersion equations (four roots for 
each homogeneous region), and the subscripts L and R indicate the values of c and 
v in the left-hand-side and the right-hand-side region respectively. The constants 
Aj have to be such that the matching conditions ( 3.2.1| > are satisfied. We can write 

0, where 



down these conditions in matrix form as Ay-Ay 
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Since there are eight free parameters Aj for each discontinuity (four in each re- 
gion), and only four conditions, matching is in principle always possible. However, 
there are additional conditions: the boundary conditions which we will discuss in 
the next chapters. These boundary conditions can be added as additional rows into 
the matrix Aj,-, to obtain a (4 + N) x 8 matrix Ay, with N = N(co) the number of 
constraints on the coefficients Aj (i.e., the number of modes that do not fulfil the 
boundary conditions and hence should be forbidden)^] Depending on the value of 
N(co), we can have the following cases. 



If N(co) < 4, then there will always exist a non- trivial solution Aj, since there 
are more degrees of freedom than constraints. Therefore, a region of the com- 



1 To give a simple example, assume that a boundary condition consists in the requirement: v g > 
at the right asymptotic infinity. One can then simply check, for the four right-hand-side modes, 
whether they fulfil this condition. If not, the mode in question is forbidden and N(co) is incremented 
by one. 
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plex a;-plane where N < 4 will represent a continuous region of eigenfrequen- 
cies. 

• If N(co) = 4, then we have an 8 x 8 system of equations and there is a non- 
trivial solution Aj only if det(A) = 0. Generally speaking, this means that 
in a region where N = 4, we might at the very most expect isolated discrete 
eigenfrequencies, but no continuous zones of solutions. 

• If N(to) > 4 then the system can be split into two (or more) 8x8 subsystems, 
each with a sub determinant A;. For a non-trivial solution to exist, each of these 
subsystems should fulfil the previous condition det(Aj) = 0. 

To summarise, we can define a non-negative function F(co) such that 

• F(co) = if N(co) < 4, 

• F{a>) = |det(A)| if N(co) =4, 

• F{cv) = £|det(A i )| if N(co)>4. 

The eigenfrequencies of the system are those values of co for which F(co) = 0. In- 
stabilities correspond to eigenfrequencies with Im(o;) > and are therefore located 
in the upper half of the complex co plane. Quasinormal modes on the other hand 
are relaxation modes, i.e., they have Im(o;) < and will therefore be located in the 
lower half complex co plane. 

The essential ingredient that is still missing to calculate the instabilities and the 
quasinormal modes in these BEC black hole configurations is the set of boundary 
conditions. Since these are different for both problems, we will postpone them to 
the relevant chapters. 
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Chapter 4 

Analysis of dynamical instabilities 



Now that we have set down the general framework, it is time to examine the linear 
dynamical stability of various (l+l)-dimensional configurations in Bose-Einstein 
condensates with steplike sonic horizons, as illustrated in Fig. 3.1 We will exam- 
ine both black hole and white hole horizons. We will start with a motivation for 
the problem. Since most of the framework has been discussed in the previous sec- 
tion, we will then immediately discuss the boundary conditions, with particular 
attention to their meaning in gravitational terms. We highlight that the stability 
of a given configuration depends not only on its specific geometry, but especially 
on these boundary conditions. Under boundary conditions directly extrapolated 
from those in standard general relativity, both black hole and white hole config- 
urations are shown to be stable. However, we show that under other (less strin- 
gent) boundary conditions, configurations with a black hole horizon remain stable, 
whereas white hole configurations develop instabilities associated with the presence 
of sonic horizons. 



4.1 Motivation 

The stability analysis presented in B57| for acoustic black hole configurations with 
fluid sinks in their interior concluded that these configurations were intrinsically 
unstable. However, the WKB analysis of the stability of horizons in [|80| suggested 
that black hole horizons might be stable, while configurations with white hole hori- 
zons seem to inevitably possess unstable modes. Regarding configurations in which 
a black hole horizon is connected with a white hole horizon in a straight line, the 
analysis presented in [ j42j ] concluded, also within a WKB approximation, that these 
configurations were intrinsically unstable, producing a so-called "black hole laser". 
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However, when the white hole horizon is connected back to the black hole horizon 
to produce a ring, it was found that these configurations can be stable or unstable 
depending on their specific parameters ||56tl57| . This suggests that periodic bound- 
ary conditions can eliminate some of the instabilities associated with the black hole 
laser. So there seems to exist a variety of different results, and hence a bit of confu- 
sion in the related literature. 

The importance of this issue is easy to see when one bears in mind that BECs 
have long been speculated to be good candidates for an experimental detection of 
Hawking radiation fl9tll3|], see chapter|7j In this context, dynamical stability is an im- 
portant condition. Indeed, Hawking radiation is (by definition 116210 the 'late-time' 
radiation of a black hole, once all transitory phenomena have vanished. However, 
if dynamical instabilities unavoidably exist, then these will necessarily start domi- 
nating from a certain time on. Depending on the relation between the time scales 
characteristic for the onset of Hawking radiation and the one for the dynamical in- 
stabilities, it might still be possible to achieve a 'metastable' experiment, i.e., one 
with a sufficient lifetime to detect the black hole radiation in spite of the existence 
of dynamical instabilities. However, this might require additional fine-tuning, and 
would in any case add an important source of background 'noise', thereby seriously 
complicating the detection of the black hole radiation. Therefore, if one could show 
that there exist black-hole configurations in BECs which are in principle devoid of 
dynamical instabilities, this would remove an important obstacle in the experimen- 
tal realisation of Hawking radiation. 

In this chapter, we will therefore try to shed some light on these issues and clear 
up some of the apparent contradictions. To simplify matters, we will consider (1+1)- 
dimensional profiles that are piecewise uniform with a single steplike discontinuity, 
as described earlier (see Fig. 3.1 1. In terms of dynamical (in)stability, there seems 
to be no crucial qualitative difference between the present case and a profile with 
smooth but narrow transitions between regions with an (asymptotically) uniform 
density distribution, at least not from a theoretical point of view [57J. Therefore, 
here we will not take into account the effects of the finite size of the transition re- 
gion and consider only idealised steplike cases. The specific way to examine the 
kind of instability we are interested in, consists basically in seeking whether, under 
appropriate boundary conditions, there are complex eigenfrequencies of the system 
which lead to an exponential increase with time of the associated perturbations, i.e., 
a dynamical instability. Throughout the section we will use a language and nota- 
tion as close to general relativity as possible. In particular, we will use boundary 
conditions similar to those imposed in the standard stability and quasinormal mode 
analysis of black holes in general relativity B76l [ 
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4.2 Boundary conditions 

In order to extract the possible intrinsic instabilities of a BEC configuration, we have 
to analyse whether there are linear mode solutions with positive Im(o;) that satisfy 
outgoing boundary conditions^] By 'outgoing' boundary conditions we mean that 
the group velocity is directed outwards (toward the boundaries of the system). The 
group velocity for a particular k-mode is defined as 

fdco\ fc 2 k+^ 2 c 2 k 3 
v g = Re ( ^77 ) = Re ( 77 ^ T _ — +v | , (4.2.1) 



dk ) \ (xi — vk 

where we have used the full Bogoliubov dispersion relation ( 2.5.3[ >j^] The physical 
idea behind this outgoing boundary condition is that we are only examining the 
intrinsic instabilities originating inside the system, since these are equivalent to dy- 
namical instabilities in general relativity. 

To illustrate this assertion, let us look at the classical linear stability analysis of a 
Schwarzschild black hole in general relativity 117611931 . When considering outgoing 
boundary conditions both at the horizon and in the asymptotic region at infinity, 
only negative Im(o;) modes (the quasinormal modes) are found, and thus the black 
hole configuration is said to be stable. If the presence of ingoing waves at infinity 
were allowed, there would also exist positive Im(o;) solutions. In other words, the 
Schwarzschild solution in general relativity is stable when considering only inter- 
nal rearrangements of the configuration. If instead the black hole was allowed to 
absorb more and more energy from outside the system (i.e., ultimately coming from 
infinity), its configuration would continuously change and appear to be unstable. 

The introduction of modified dispersion relations adds an important difference 
with respect to the traditional boundary conditions used in the linearised stability 
analysis of a general relativistic black hole. In a BEC black hole, one boundary is the 
standard asymptotic region, just like in general relativity. To determine the other 
boundary, the following must be taken into consideration. For a relativistic or hy- 
drodynamic (quadratic) dispersion relation, nothing can escape from the interior of 
a sonic black hole: the acoustic behaviour is analogous to linearised general relativ- 
ity. But due to the superluminal corrections, the high-frequency effective co-moving 
speed of sound has a higher value than the zero-frequency speed of sound c, 
see 2.5.8 and Fig. 2.3 These superluminal corrections are incorporated in the full 



Recall that we follow the convention u w = E ; Aje'( k i x ~ wt \ see eq. ( [2.5.5) or j3.2.2\ . 



2 Note that this definition is equivalent to the simpler definition [ 2.5.7) only for the modes that 
behave like eq. 1 2.5.6) in the hydrodynamic limit. 
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Bogoliubov dispersion relation but not in the relativistic approximation. Due to this 
superluminal character, information from the interior of the acoustic black hole can 
escape through the horizon and affect its exterior. Therefore, since we are precisely 
interested in deviations from the relativistic regime, this permeability of the horizon 
must be taken into consideration. The second boundary then is not the black hole 
horizon itself, but the internal 'singularity', which corresponds to x — ► — oo in our 
model. We repeat that, in our model, all magnitudes are regular in this region, so 
there is no singularity in the strict sense. The outgoing boundary condition there 
reflects the fact that no information can escape from the singularity. 

There is another complication that deserves some attention. In the case of a hy- 
drodynamic dispersion relation, the signs of Vg = v ± c and Im(fc) = (v ± c)lm(co) 
coincide for Im(o;) > 0. For example, in an asymptotic x — > +oo subsonic region, an 
outgoing fc-mode has v g > 0, so that Im(fc) > 0. Such a mode is therefore damped 
towards this infinity, giving a finite contribution to its norm. In the linear stability 
analysis of black hole configurations, stable dynamical modes therefore correspond 
to non-normalisable perturbations, while unstable modes correspond to normalis- 
able ones. When considering modified dispersion relations, however, this associa- 
tion no longer holds. In particular, with the BEC dispersion relation, in an asymp- 
totic x — ► +oo region, among the unstable (lm(co) > 0) outgoing (v g > 0) fc-modes, 
there are modes with Im(fc) > as well as modes with Im(fc) < 0. An appropri- 
ate interpretation of these two possibilities is the following. The unstable outgoing 
modes that are convergent at infinity (those with Im(fc) > 0) are associated with per- 
turbations of the system that originate in an internal compact region of the system. 
Unstable outgoing modes that are divergent at infinity are associated with initial 
perturbations acting also at the boundary at infinity itself. 



Take for example a black hole configuration of the form described in Fig. 3.1 



The right asymptotic region, which can be interpreted as containing a 'source' of 
BEC gas in our analogue model, simulates the asymptotic infinity outside the black 
hole in general relativity. The convergence condition at this right-hand side then 
implies that the perturbations are not allowed to affect this asymptotic infinity ini- 
tially. However, for the left asymptotic region, this condition is less obvious. In our 
BEC configuration, this left asymptotic region can be seen as representing a 'sink'. It 
corresponds somehow to the general relativistic singularity of a gravitational black 
hole. The fact that in general relativity this singularity is situated at a finite distance 
(strictly speaking, at a finite amount of proper time) from the horizon, indicates that 
it might be sensible to allow the perturbations to affect this left asymptotic region 
from the start. We will therefore consider two possibilities for the boundary con- 
dition at x — > —oo. (a) Either we impose convergence in both asymptotic regions, 
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thereby eliminating the possibility that perturbations have an immediate initial ef- 
fect on the sink, or (b) we allow the perturbations to affect the sink right from the 
start, i.e., we do not impose convergence at the left asymptotic region. The option of 
imposing convergence at the left asymptotic region could be interpreted as exclud- 
ing the influence of the singularity on the stability of the system. In other words, 
condition (b) would then be equivalent to examining the stability due to the com- 
bined influence of the horizon and the singularity, while under condition (a) only 
the stability of the horizon is taken into account. 

As a final note to this discussion, since we are interested in the analogy with 
gravity, we have assumed an infinite system at the right-hand side. In a realistic 
condensate other boundary conditions could apply, for example taking into account 
the reflection at the ends of the condensate (see e.g. ||56ll57| V 



4.3 Case by case analysis and results 

Now that we have discussed the boundary conditions, we have all the ingredients 
necessary to perform the numerical analysis described in the previous section. Since 
we are looking for instabilities (Im(o;) > 0), we plot the function F(to), defined in 



Sec. 3.2 in the upper half of the complex plane, and look for its zeros. Each of these 
zeros indicates an unstable eigenfrequency, and so the presence (or absence) of these 
zeros will indicate the instability (or stability) of the system. 

All figures are given in dimensionless units. Typical values for the velocity of 
sound c in BECs range from 0.5 to 10mm/ s, while the healing length £ typically lies 
between lOOnm and 2fim. Taking, for concreteness, c=lmm/s and £=l^m, our nu- 
merical results can be translated to realistic physical numbers by using milliseconds 
as natural units. For example, the typical timescale for the development of an insta- 
bility with Im(o;) — 0.1 would be about 10ms. We have checked within a reasonable 
range of values that our results do not depend on the particular quantities chosen 
for the velocities of the system. 

We will now discuss the specific configurations that have been analysed case by 
case. 



4.3.1 Black hole configurations 



Consider a flow accelerating from a subsonic regime on the right-hand side to a 
supersonic regime on the left, see Fig. 3.1 For observers at the right-hand side, 
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this configuration possesses an acoustic black hole horizon. For such configurations 
with a single black hole horizon, when requiring convergence in both asymptotic re- 
gions [case (a)], there are no zeros (see Fig. 4.1 ), except for two isolated points on the 
imaginary axis (see Fig. 4.3 which is a zoom of the relevant area in Fig. 4.1 1. Note 
that we always check the existence of a zero by zooming in on the area around its 
location up to the numerical resolution of our program. These two isolated points 
seem to be of a rather special nature. They are located at the boundary between 
regions with different number N of forbidden modes in the asymptotic regions. The 
zeros that we will find for other configurations are of a totally different nature: they 
are either sharp vanishing local minima of F(co) situated well inside an area with a 
constant value of N (N = 4 to be precise), or continuous regions where the algebraic 
system represented by the matrix Au is underdetermined (and hence N < 4, see 
Sec. 3.2). A technical analysis of these special points can be found in [16]. Here, let it 
suffice to mention that points of this kind are always present in any flow, indepen- 
dently of whether it reaches supersonic regimes or not. Hence it is clear that they do 
not correspond to real physical instabilities, since otherwise any type of flow would 
appear to be unstable. Accordingly, in the following, we will not take these points 
into consideration. When we assert that a figure is devoid of instabilities, we will 
mean that the function F(co) has no zeros except for the special ones just mentioned. 



Fig. 4.1 also shows that the system remains stable even when eliminating the con- 
dition of convergence at the left-hand side, i.e., at the sink simulating the singularity 
[case (b)]. 

To sum up, configurations possessing a (single) black hole horizon are stable 
under the general boundary conditions that we have described, i.e., outgoing in 
both asymptotic regions and convergent in the upstream asymptotic region, inde- 
pendently of whether convergence is also fulfilled in the downstream asymptotic 
region or not. 



4.3.2 White hole configurations 

Let us now consider flows decelerating from a supersonic regime at the right-hand 
side to a subsonic one on the left. From the point of view of observers at the left- 
hand side, the geometric configuration possesses a white hole horizon. In general 
relativity, a white hole corresponds to the time reversal of a black hole. Therefore, in 
general relativistic terms, unstable modes of a white hole configuration would cor- 
respond to stable (quasinormal) modes of the black hole. From this point of view, 
any quasinormal mode found in the analysis of a black hole configuration (see next 
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Black hole 

Csuper = 0.7, C su b = 1.8. 

iV(oi) log[FM] 



White hole 

*-sub = 1-8, <-super = 0.7. 





log[F(<a)] 




Re ((d) 



2 4 
Re(m) 



Figure 4.1: Plots of N(co) and logF(o;) for black hole (left) and white hole (right) 
configurations (white points or regions, where F(co) —> 0, represent instabilities). 
The speed of sound c is indicated for each region and the velocity v is then obtained 
from the constraint vc 2 = 1. Note that c > 1 corresponds to a subsonic region, c < 1 
to a supersonic region. In the upper pictures [case (a)], convergence has been im- 
posed in both asymptotic regions. In the lower pictures [case (b)], convergence has 
been imposed only in the upstream asymptotic region, but not at the sink or 'singu- 
larity' downstream. It is seen that black hole configurations are stable in both case 
(a) and (b), as are accelerating subsonic flows in general (see Fig. |4.2| below). White 
hole configurations are stable in case (a), but develop a huge continuous region of 
instabilities in case (b). Only a small strip of instabilities subsists in the decelerating 



subsonic flow (see Fig. 4.2 below), indicating that the major part of this unstable re- 
gion is a genuine consequence of the existence of the white hole horizon. Note that 
continuous regions of instability correspond to N(co) < 4. 
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Accelerating subsonic flow 

Csubl = 1-8, c sub2 = 1.9. 



tag[F(m)] 



Decelerating subsonic flow 

Csubl = 1-9, C su b2 = 1-8. 

IV (to) log[F(io)] 






n 



Re(co) 
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Figure 4.2: Plots of N(to) and log F(co) for accelerating (left) and decelerating (right) 
subsonic flows (i.e., without acoustic horizons). In the upper pictures [case (a)], 
convergence has been imposed in both asymptotic regions. In the lower pictures 
[case (b)], convergence has been imposed only in the upstream asymptotic region. 
The decelerating subsonic flow exhibits a small strip of instabilities in case (b), which 
was already found in the white hole case (note that these continuous regions of 
instability correspond to N(co) < 4). 



log[F(a>)] 




-0.1 0.0 
Re(co) 



Figure 4.3: Two special zeros of the function F(co) appear in the stability analy- 
sis of a black hole configuration (this plot is a zoom of the corresponding region in 



Fig. 4.1 1. They are located at the boundary between regions with different number 



N of prohibited modes. These points appear in the stability analysis of all configu- 
rations studied, including purely subsonic ones, and hence do not seem to represent 
real physical instabilities of the system. 
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chapter) could immediately be interpreted as an instability of the corresponding 
white hole. Note that such a point of view implies that the boundary conditions ap- 
propriate for the analysis of white hole configurations correspond to only allowing 
ingoing waves at the boundaries, due to time reversal. 

But from the point of view of acoustic models in a laboratory, the analysis of the 
intrinsic stability of the flow (under the outgoing boundary conditions described 
above) is also interesting. This analysis also has particular relevance with regard to 
configurations with two horizons (see |fl6|V 



In Fig. 4.1 we see that under outgoing boundary conditions the flow is stable 
when convergence is required in both asymptotic regions [case (a)], but exhibits a 
continuous region of instabilities at low frequencies when convergence is fulfilled 
only at the right-hand side [case (b)]. Indeed, in this continuous region, N = 3, in 
other words the algebraic system AuAu = is underdetermined and any frequency 
is automatically an eigenfrequency. 

When looking at the case of a completely subsonic flow undergoing a deceler- 
ation (see Fig. 4.2| , we find something similar. The system is clearly stable when 
convergence is imposed at the left-hand side, i.e. in case (a). Without convergence 
at the left-hand side, case (b), there is a small continuous strip of instabilities which 
corresponds, as in the white hole case, to a region where N = 3. However, this 
region is located at relatively high frequencies, disconnected from co ~ 0. It seems 
that part of the continuous region of instabilities found in the white hole configu- 
ration has its origin merely in the deceleration of the flow (giving rise to this high- 
frequency strip). But there is still a complete region of instabilities that is a genuine 
consequence of the existence of a white hole horizon. In fact, by decreasing the heal- 
ing length parameter £, the strip associated with the mere deceleration moves up to 
higher and higher frequencies, becoming less and less important as one approaches 
the hydrodynamic or acoustic limit. The main continuous region of instabilities that 
can be associated with the horizon, i.e., the white region expanding from co ~ 0, on 
the other hand, does not change its character in this process. 

We can therefore conclude the following with regard to decelerating configura- 
tions. When convergence is fulfilled downstream, the configuration is stable, re- 
gardless of whether it contains a white hole horizon or not. When this convergence 
condition is dropped, there is a tendency to destabilisation. In the presence of a 
white hole horizon, the configuration actually becomes dramatically unstable, since 
there is a huge continuous region of instabilities, and even perturbations with ar- 
bitrarily small frequencies destabilise the configuration. In the absence of such a 
horizon, only a small high-frequency part of this unstable region subsists. 
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Re(eo) 



Figure 4.4: A discrete set of instabilities appears in a black hole configuration when 
the left asymptotic region simulating the singularity is replaced by a wall or sink 
at a finite distance inside the black hole. In this plot we have used c su b = 1.9 and 
Csuper = 0.7 (with the corresponding v = 1/c 2 ). In addition we have taken L = 6 as 
the size of the finite internal region. 



4.3.3 Black hole configurations with modified boundary conditions 

We have seen earlier that configurations with a black hole horizon do not possess 
instabilities under either of the sets of boundary conditions that we considered. 
However, instabilities can show up when instead of extending the black hole con- 
figuration towards — oo, a wall (or sink) is introduced at a finite distance inside the 
supersonic region, leading to a possible bouncing and amplifying effect of some 
modes inside the black hole. Such a case would be described by boundary condi- 
tions different from the ones we have considered so far. For example, by replacing 
the boundary conditions at the left-hand side by 8\ x= _i = 0, it can be seen from 



Fig. WA that a discrete set of unstable modes appears. 



4.4 Discussion 



Let us start by discussing the stability of configurations with an acoustic black hole 
horizon in analogue systems that incorporate superluminal dispersion relations. We 
have seen that by requiring purely outgoing and convergent boundary conditions in 
both asymptotic regions, these configurations do not show any signs of instability. 
The same applies when dropping the convergence condition downstream (i.e., on 
the left-hand side simulating the singularity). This seems to contradict the results 
in ||57ll . There, the existence of a future (spacelike) singularity inside the black hole, 
from which no information is allowed to escape, was implemented by introducing 
a sink in the supersonic region at a finite distance from the horizon. Then, it was 
found that there were discrete instabilities in the system. However, these instabil- 
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ities correspond to the following particular set of boundary conditions: i) At the 
asymptotic region, only convergent boundary conditions were imposed, without 
any condition about the direction of propagation (in- or outgoing) of the perturba- 
tions; ii) At the sink, two types of boundary conditions were required, specifically 
designed for dealing with symmetric and anti-symmetric configurations. In our 
language these boundary conditions correspond to {0'lsmk = 0,n'\ S mk = 0} an d 
I sink — 0/ n lsink = 0}, respectively Comparison of the results obtained with these 
boundary conditions ||57| with the results that we obtained here, leads to two impor- 
tant conclusions. On the one hand, the unstable modes found in B57l have ingoing 
contributions at the asymptotic region. On the other hand, the boundary conditions 
at the sink are such that they combine outgoing and ingoing contributions. The sink 
implementation causes waves reaching the sink to bounce back towards the hori- 
zon. These two facts are responsible for the unstable behaviour of such black hole 
configurations. If no energy is introduced into the system from the asymptotic re- 
gion (in other words, if only outgoing perturbations are allowed) and moreover any 
bouncing at the sink is eliminated, then these configurations are stable. This is in 
agreement with the result found in [80 [. 

In the case of configurations with a white hole horizon, we have seen that with 
outgoing and convergent boundary conditions in both asymptotic regions, there are 
no instabilities in the system. However, when eliminating the convergence con- 
dition in the downstream asymptotic region (at the sink), one finds a continuous 
region of instabilities surrounding co = 0. Thus, we see that these white hole config- 
urations are stable only when the boundary conditions are sufficiently restrictive. 

When analysing configurations connecting two different subsonic regions, i.e., in 
the absence of acoustic horizons, we have also seen that, again, when convergence is 
required at the left-hand side, they are stable. But when this convergence condition 
is relaxed, globally decelerating configurations tend to become unstable, whereas 
globally accelerating ones remain stable. The instabilities of these decelerating con- 
figurations without horizons (i.e., purely subsonic ones) show up, however, in a 
small strip at high rfequencies. In contrast, white hole configurations present insta- 
bilities for a wide range of frequencies, starting from arbitrarily small values. This 
shows that the presence of a white hole horizon drastically augments the instability 
of the configuration. 
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4.5 Summary and conclusions 

To sum up, we have shown the high sensitivity of the stability not only to the type of 
configuration (the presence of a horizon, the accelerating or decelerating character 
of the fluid), but particularly to the boundary conditions. In general, the bound- 
ary conditions that we have considered are directly inspired by similar analyses in 
general relativity, and so we have always imposed outgoing boundary conditions, 
i.e., only the possible existence of instabilities arising within the system should be 
taken into account. For similar reasons, we argued that convergence must also be 
imposed at the upstream (right-hand side) asymptotic region, i.e., at the source of 
the condensate. 

When moreover requiring convergence at the left-hand or downstream asymp- 
totic region (which can be interpreted as simulating a singularity, since nothing can 
escape from it towards the horizon), both black hole and white hole configurations 
are stable (and also the combination of both into a black hole-white hole configu- 
ration, see [16]). When relaxing this convergence condition at the left-hand side, 
configurations with a single black hole horizon remain stable, whereas white hole 
configurations (and black hole-white hole configurations, see ||T6l ) develop instabil- 
ities not present in flows without horizons. 

Some final comments with relevance for the remainder of this thesis might be 
useful. We have seen that black hole configurations are stable under rather general 
boundary conditions. This implies the following. 

First of all, we have explicitly worked with the Bogoliubov dispersion relation 
valid for a BEC. However, since we have checked a wide range of parameters, the 
results that we have obtained can be extended to superluminal dispersion relations 
of the Bogoliubov type in general. One can then tentatively conclude that adding 
such superluminal corrections does not radically alter the dynamical stability of a 
relativistic black hole, whether it is an acoustic black hole or a gravitational one. 

Note, however, that this result is 'fragile' in the following sense: configurations 
with a black hole horizon and a bounce inside the system can become unstable. 
Also, the stability is guaranteed if only ingoing modes are allowed. Although this 
is a sensible condition from a theoretical analogue-gravity point of view, it is not so 
clear that this condition makes sense in a BEC experiment. Nevertheless, even in 
a realistic BEC experiment, it should be possible to make a condensate sufficiently 
long in order for reflection at its boundaries to be negligible for the duration of the 
experiment |43tl57| . Let us also remark that dynamical instabilities are not the only 
instabilities that can appear in an experiment. Surface instabilities, for example, 
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could also appear. Such surface instabilities are sometimes a desired effect because 
they can lead to vortex formation. However, they are also generally easy to avoid, 
in particular in effectively one-dimensional condensates [102J. 

Second, since we have shown that there can exist dynamically stable acoustic 
black hole configurations, it indeed makes sense to analyse the stable dynamical 
modes of such configurations, i.e., the quasinormal modes, as we will do in the next 
chapter. 

Finally, it also means that such configurations can indeed be expected to repro- 
duce analogue (i.e., phononic) Hawking radiation. Although this obviously does 
not solve all experimental issues, at least it removes the need to recur to metastable 
configurations, fine-tuned in the sense that the dynamical instabilities would not 
disrupt the configuration nor obfuscate the black hole radiation in the time neces- 
sary to achieve the experiment. 
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Chapter 5 

Quasinormal modes 



In this chapter, we perform a quasinormal mode analysis of black hole configura- 
tions in Bose-Einstein condensates (BEC). In this analysis we use the full Bogoliubov 
dispersion relation (see Sec. 2.5 I, not just the hydrodynamic or relativistic approxi- 
mation. We restrict our attention to (l+l)-dimensional black hole configurations in 
BEC with steplike discontinuities of the type illustrated in Fig. 3.1 For this case we 
show that in the hydrodynamic or relativistic approximation, quasinormal modes 
do not exist. The full dispersion relation, however, allows the existence of short- 
lived quasinormal modes. Remarkably, the spectrum of these modes is not discrete 
but continuous. 



5.1 Introduction 

In the previous chapter, the stability of several (l+l)-dimensional configurations in 
BECs was studied. These configurations were chosen in analogy with the geome- 
tries associated with gravitational black holes. Here we will elaborate further on 
black hole configurations, which were shown to be devoid of instabilities, and exam- 
ine their quasinormal or relaxation modes. The analysis we carry out is analogous 
to the standard quasinormal mode (QNM) analysis in gravity Il76ll93l . According to 
general relativity, when black holes are perturbed, they emit gravitational waves in 
the following way (see Fig. 5.1[ >. The initial regime depends mainly on the concrete 



form of the perturbation. Then, an oscillatory decay phase is attained in which the 
precise form of the emitted wave depends only on the properties of the black hole 
itself. Finally, a polynomial tail characterises the return to equilibrium. In the in- 
termediate oscillatory decay phase, a discrete set of complex frequencies is excited: 
the quasinormal modes of our concern. Since quasinormal modes are modes of de- 
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Figure 5.1: Typical wave form of gravitational radiation emitted by a perturbed 
black hole as a function of time t. The initial pulse depends on the concrete form 
of the perturbation. Then, a damped oscillation consisting of a discrete set of fre- 
quencies, the quasinormal modes, sets in. (Often, as in this figure, a single frequency 
clearly dominates). Finally, a power-law tail takes over at late time. 



cay, i.e., of energy dissipation, they are outgoing: they move towards the exterior 
of the spacetime region connected to an asymptotic observer. In general relativity 
this means that their group velocities are directed outwards both in the asymptotic 
region and at the horizon. Moreover, as we will illustrate, in general relativity this 
also automatically implies that they are non-normalisable or divergent. 

As we will briefly point out, in standard general relativity in 1+1 dimensions, 
quasinormal modes do not exist. In higher dimensions, a discrete spectrum appears 



|93[ , see Fig. 5.2 A general investigation of quasinormal modes in analogue 
black holes was carried out in [25] . Results qualitatively identical to the general 
relativistic case were found. Here the essential element that we add is the influence 
on these QNMs of a non-relativistic contribution to the dispersion relation. 



5.2 Boundary conditions 

The general framework to calculate the quasinormal modes in these BEC black hole 
configurations is the same as for the calculation of the instabilities, and has been 
described in Chapter [3] However, we still need to complement this by the adequate 
boundary conditions. We will first take a look at the boundary conditions that are 
used to calculate QNMs in general relativity. This can be immediately extrapolated 
to BECs in the hydrodynamic limit. As we will show, in this hydrodynamic limit 
in 1+1 dimensions, there exist no QNMs. Finally, we will come to the boundary 
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Figure 5.2: Discrete set of quasinormal modes in the complex frequency-plane for 
a Schwarzschild black hole and different values of the angular momentum index L. 
Figure taken from l93ll . (Note that we use the opposite convention for the sign of 
Im(cu).) 



conditions with the full dispersion relation. 



5.2.1 Quasinormal modes of a gravitational black hole 

Quasinormal modes of a gravitational black hole are the relaxation modes or modes 
of energy dissipation that characterise the pulsations of the black hole after pertur- 
bations initiated internally in the spacetime surrounding the black hole. 

Since the QNMs are relaxation modes, they must decay in time and moreover be 
outgoing. The first condition means that, for the sign convention followed here (see 
e.g. eq. ( |2.5.5| >), they can be represented by complex frequencies to with Im(o;) < 0. 



The outgoing character must be imposed in the asymptotic region. At the horizon, 
QNMs are also required to be outgoing because in general relativity nothing can es- 
cape through the horizon^] So the boundary conditions that are imposed in general 
relativity to find the QNMs of a black hole are simply that these modes should be 



*As we mentioned in the introduction to this chapter, in QNM terminology, 'outgoing' means 
directed towards the exterior of the spacetime region connected to an asymptotic observer. In other 
words, at the horizon, 'outgoing' means directed towards the singularity. 



60 



CHAPTER 5. QUASINORMAL MODES 



outgoing both in the asymptotic region and at the horizon. 

Because QNMs of a gravitational black hole are outgoing, they are also diver- 
gent. Indeed, due to the quadratic form of the relativistic dispersion relation, co 2 oc 
k 2 , the group velocity v g = dco/ dk and the imaginary part of k are related by 

Im(fc) cx1m(cv)/v g , (5.2.1) 

and therefore, since QNMs have Im(o;) < 0: 

sign[Im(fc)] = -signfz^] . (5.2.2) 



5.2.2 Hydrodynamic limit in 1+1 dimensions: absence of QNMs 

In BECs in the hydrodynamic limit, the dispersion relation is also quadratic: (co — 
vk) 2 = c 2 k 2 , see eq. ( |2.5.4| >, and therefore decaying outgoing modes are also auto- 



matically divergent. Indeed, 

v g = v±c, Im(fc) = Jm(co)/Vg , (5.2.3) 

and therefore sign[Im(fc)] = — signfy^] as before. Let us discuss whether QNMs can 
exist in this case. 

There are two ways to deal with this question. Either one can define a connection 



matrix based on the hydrodynamic limit of the matching conditions (3.2.1 1, and see 
that the appropriate outgoing boundary conditions can never be fulfilled, see |[2~T|. 
Alternatively, one can reason as follows. Because of the presence of the horizon, we 
know that, in the hydrodynamic limit, none of the modes at the left-hand side can 
move upstream across the horizon, i.e., they cannot connect to outgoing modes at 
the right-hand side. In other words, when the (x < 0) region is supersonic, it be- 
comes causally disconnected from the subsonic region (x > 0). But then, just like in 
general relativity, the requirement that QNMs must be outgoing should be imposed 
at the horizon, and not at the left asymptotic infinity. Because of the homogeneity 
of the profile, the mode solutions for each particular frequency co can be written as 
piecewise sums of plane waves: 

ni = Aie , 'fa*- & *> + A 2 i^ x -^ , (5.2.4) 
9 1 = Biffa*-"*) + B 2 e i(klX -^ . (5.2.5) 

Because of this plane-wave character, none of these modes can be outgoing both at 
the right asymptotic region and at the horizon. QNMs can therefore not exist in the 
hydrodynamic limit in 1+1 dimensions. 
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Although the argument presented here seems to depend on the homogeneity of 
the profiles, it is actually more general. Indeed, in the hydrodynamic limit, the ge- 
ometric description developed in section 2.4 is valid. It is well known from general 
relativity that all (l+l)-dimensional metrics are conformally flat. In addition, the 
d'Alembertian equation in two dimensions is conformally invariant ||26l . Therefore, 
all the solutions of the system are conformally equivalent to plane waves, and so 
they do not satisfy the requirement of being outgoing both in the asymptotic re- 
gion and at the horizon. Hence QNMs cannot exist in (l+l)-dimensional relativistic 
configurations, and this argument is valid both for general relativity and for acous- 
tic black holes in the hydrodynamic limit. The situation is completely different in 
3+1 dimensions, in which both general relativistic and hydrodynamic acoustic black 
holes exhibit a discrete QNM spectrum, as we mentioned earlier. 



5.2.3 Boundary conditions with the full dispersion relation 

A crucial feature of the superluminal character of the dispersion relation is that it 
makes the horizon permeable, so that outgoing boundary conditions should be im- 
posed at both asymptotic infinities (rather than at the right asymptotic infinity and 
at the horizon). Moreover, since we have two homogeneous regions, the ingoing or 
outgoing character of a mode is conserved along each side of the horizon (x = 0). 
An outgoing mode is then one that has a group velocity v g < for x < (and in 
particular for x — > — oo) and Vg > for x > (and in particular for x — > +oo). Vg is 



obtained from the full Bogoliubov dispersion relation, see eq. 4.2.1 



,'dco\ fc 2 k+^ 2 c 2 k 3 
''W^lg = Re 0,-vk (5 ' 2h) 



Note that the immediate connection between the sign of Vg and that of Im(fc) which 
existed in the hydrodynamic limit clearly does not subsist, so the automatic diver- 
gent character of outgoing modes is lost. 

In any case, quasinormal modes should be outgoing, and this is the essential 
boundary condition. Additionally, one might wonder whether the QNMs will be 
convergent, or divergent after all (in spite of the modification of the dispersion re- 
lation). An easy way to check this is by applying convergence as an additional 
boundary condition and comparing the resulting spectrum with the one obtained 
without this additional condition. 
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5.2.4 Summary 

To sum up, quasinormal modes should be outgoing. In general relativity and in 
acoustic black holes in the hydrodynamic approximation, this outgoing condition 
must be imposed at the external (right-hand side) asymptotic infinity and at the 
horizon. Moreover, the outgoing character of QNMs automatically implies a diver- 
gent behaviour in the asymptotic regions. However, in both these cases, there are no 
QNMs in 1+1 dimensions, whereas in 3+1-dimensional black holes, a discrete QNM 
spectrum appears. 

In BECs with the full modified dispersion relation, outgoing modes are no longer 
automatically divergent. The boundary condition to determine QNMs in the (1+1)- 
dimensional BEC black hole configurations under study is simply that they should 
be outgoing in both asymptotic regions. 

Application of the method developed in chapter [3] then leads to the results that 
we will now discuss. 



5.3 Results and discussion 

The quasinormal mode spectrum of a (l+l)-dimensional acoustic black hole in a 



BEC, obtained with the full dispersion relation, is illustrated in Fig. 5.3 Remember 
from the previous chapter that such a black hole did not possess unstable eigenfre- 
quencies, so it indeed makes sense to look for stable modes. Surprisingly, as can 



be seen from Fig. 5.3 not only does such a QNM spectrum indeed show up, but 
this spectrum is actually continuous. Indeed, in the whole continuous region where 
N(to) =3, the algebraic system composed of the matching and boundary conditions 
is underdetermined. Therefore, according to the procedure outlined in chapter[3j ev- 
ery frequency in this region is automatically a quasinormal mode. We have checked 
that no further (discrete) quasinormal modes appear in the N = 4 region. 

We have also checked whether these quasinormal modes are convergent or not. 
We have found that the quasinormal modes here contain both convergent and diver- 
gent contributions, and actually disappear as soon as either type of contribution is 
prohibited (i.e., as soon as a strict convergence or divergence condition is imposed 
at either side of the configuration). Overall, this makes these quasinormal modes 
divergent. However, this divergent character is not obtained in the same automatic 
sense, i.e., it is not due to a strict relation between the signs of v g and Jm(k), as in the 
case of the QNMs that are found in the (3+l)-dimensional relativistic and hydrody- 
namic cases discussed earlier. 
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Figure 5.3: Plot illustrating the quasinormal mode spectrum for a (l+l)-dimensional 
black hole configuration in a BEC. We represent the number N of constraints (i.e., 
the number of forbidden modes) following from the boundary condition that QNMs 
should be outgoing. The QNM spectrum consists of the continuous region where 
N(co) <4, marked with a red arrow. The numerical values used for this plot, in units 
such that the healing length £ = 1 and the speed of sound c = 1 are v = 0.7 in the 
subsonic region and v = 1.8 in the supersonic region. 



We have furthermore calculated the QNM spectrum for other black-hole-like 
configurations that were shown to be devoid of instabilities in lrT6l , and we have 
checked that a similar continuous spectrum of QNMs appeared in all of them. This 
indicates that the continuous character of the QNM spectrum is an essential conse- 
quence of the modified dispersion relation, and not just of the particular characteris- 
tics of a precise configuration. In fact, the presence of these QNMs can essentially be 
traced back to the permeability of the horizon due to the superluminal character of 
the dispersion relation. In our particular model, this is easy to see. Indeed, impos- 
ing an outgoing boundary condition at the horizon would prohibit the existence of 
QNMs due to the piecewise homogeneity of the profile, regardless of whether one 
considers the full dispersion relation or only its acoustic limit. 

Finally, a straightforward dimensional estimate shows that these quasinormal 
modes are short-lived. Indeed, as can be seen from Fig. 5.3 the obtained quasinor- 
mal modes typically have Im(o;) ~ 1 (in dimensionless units). This means that these 
quasinormal modes are characteristic of a relaxation at the microscopic level, i.e., at 
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the scale of the healing length (which, as we have seen in section 2.5 is also char- 
acteristic of the Lorentz symmetry breaking scale). Adding dimensions, the lifetime 
of such quasinormal modes is of the order of £/c, the healing length of the conden- 
sate divided by its speed of sound, typically around Ijim and 1mm/ s, respectively, 
resulting in lifetimes £qnm ~lms. 



5.4 Summary and conclusions 

We have examined the quasinormal modes of (l+l)-dimensional black hole configu- 
rations in BECs, and laid particular emphasis on the importance of the quartic form 
of the full dispersion relation, compared to the quadratic form of a hydrodynamic 
or relativistic one. In the latter case, outgoing relaxation modes are automatically 
divergent. The modification of the dispersion relation means that this relation is no 
longer automatic. Nevertheless, the quasinormal modes in the systems discussed 
here are also divergent. 

More importantly, the QNM spectrum that was obtained in these systems with 
a modified dispersion relation consists of a continuous region in the complex fre- 
quency plane. The lifetime of these QNMs was estimated to be of the order of £/c, 
the healing length of the condensate divided by its speed of sound. These quasinor- 
mal modes are therefore characteristic of a relaxation at the microscopic scale, with 
typical values of the lifetime in the millisecond range. 

The importance of this result lies in the following observation. Both in general 
relativity and in BECs in the hydrodynamic limit, QNMs simply do not exist in 1+1 
dimensions, while in higher dimensions they form a discrete spectrum. So due to 
the modification of the dispersion relation, the QNM spectrum in (l+l)-dimensional 
BEC black holes changes from non-existent to a continuous region of frequencies. 
It should be noted that this result does not depend on fine-tuned choices of the 
parameters of the configuration, but is a generic consequence of the permeability 
of the horizon in the presence of superluminally modified dispersion relations. It 
is therefore a straightforward speculation that the discrete spectrum which is ob- 
tained in the usual QNM analysis in standard 3+1-dimensional general relativity 
will also develop continuous regions if it turns out that there exist superluminal 
modifications of the dispersion relations at high energies. The lifetimes of these ad- 
ditional, continuous bands of quasinormal modes will probably be extremely short, 
since they would characterise a relaxation at the microscopic or '(trans)planckian' 
level, after which a relaxation at the classical level of general relativity takes over. 
Their existence could therefore in principle be a signal of transplanckian (quantum- 
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gravitational) physics encoded in the spectrum of gravitational waves emitted by 
black holes. 
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Intermezzo: Lorentz violation and 
modified dispersion relations 



In the previous chapters, we have studied the stability and quasinormal modes of 
BEC black holes, and speculated that the results could be extrapolated to quantum 
gravity. In this brief intermezzo, we will make a more precise argument for this 
extrapolation to quantum gravity in scenarios with Lorentz symmetry breaking — 
see ||68t |85| for general introductions on Lorentz violations at high energy and their 
possible detection. 

From a phenomenological point of view, one of the first modifications that quan- 
tum gravity could impose on classical physics is a modification of the dispersion 
relations for particles. Assuming that this effect is perturbative at low energy, this 
can be written in the following way (in a given observer's frame). The usual locally 
Lorentz invariant energy-momentum dispersion relation 

E 2 = m 2 c A + c 2 p 2 (6.0.1) 

is replaced by a modified dispersion relation 

E 2 = m 2 c 4 + c 2 p 2 + £/( ! V|p| ; Ej~ l , (6.0.2) 

i 

where Ei indicates the energy scale of local Lorentz symmetry breaking and 
the strength of z th order Lorentz symmetry breaking at the energy E^. Here, we are 
in the first place interested in examining whether there exist concrete physical phe- 
nomena that break Lorentz invariance around a certain well-defined fundamental 
energy scale. In order to extract concrete boundaries on the possible existence of 



such phenomena, eq. ( |6.0.2[ > is then the adequate starting point. Actually, from a 
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'naive' quantum gravity perspective, this Ei — if it exists — is expected to be around 
the Planck scale, and so the previous relation is usually written as 

E 2 = m 2 c* + c 2 p 2 + J^f^c'lpl'E 2 - 1 , (6.0.3) 

i 

where El has been replaced directly by the Planck energy Ep/. 

The fact that we immediately wrote, for example, f^p 2 in the previous expres- 
sions rather than the more general p x pl means that we are only looking for 
Lorentz symmetry breaking in the boost sector. The motivation is the following. 
Breaking the rotational subgroup of the Lorentz group would either entail observ- 
able effects at everyday macroscopic scales or imply that the boost subgroup is also 
broken. The opposite is not true: the boost subgroup can be broken without the 
rotation subgroup being broken too, and without leading to effects that should have 
been detected already at energies experimentally available. Therefore, phenomeno- 
logically speaking, it is sufficient to focus on Lorentz symmetry breaking in the boost 
sector. Also note that/W can be particle- and helicity-dependent. 

The scheme just summarised is most adequate for effective field theory approaches 
(EFTs) |68l |85| . There are two reasons [0 why we are interested here in studying 
Lorentz symmetry breaking from an EFT perspective rather than for instance Dou- 
bly Special Relativity (DSR)j^] First, the EFT approach is the simplest one from a 
phenomenological point of view. It is also an approach for which the framework is 
sufficiently well understood and precisely defined to allow the extraction of sensi- 
ble constraints on concrete sets of physically meaningful parameters. The second 
reason originates in analogue gravity: the form of Lorentz symmetry breaking dis- 
cussed here is the one which corresponds to how the effective low-energy Lorentz 
symmetry in condensed matter and other systems is broken at high energies. In 
other words, there are plenty of examples in nature where an approximate Lorentz 
symmetry is obtained in the low-energy limit of an equation of the form ( |6.0.2| > (and 
in particular, ( 6.0.4| > below). It therefore seems worth examining whether the same 



scheme could be applicable to gravity as well. 

The standard approach within the EFT framework is the following. One con- 
structs a Lagrangian using the standard model of particle physics in combination 
with Lorentz violating operators of physical interest. These additional operators are 



1 Doubly Special Relativity [2: 84 J consists in the proposal that the Lorentz group might be de- 
formed at high energies rather than broken, by the introduction of a second observer-independent 
quantity (apart from the speed of light), namely a minimum length. DSR leads to modified dis- 
persion relations, though generally more complicated than the form we discuss here, without the 
equivalence of reference frames being broken. 
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classified by their mass dimension. For example, all renormalisable terms that can 
be added to the standard model Lagrangian (forming the so-called minimal Stan- 
dard Model Extension or mSME) have been studied more than a decade ago |38|. 
They are all of mass dimension 3 or 4 and are further classified by their behaviour 
under CPT transformations. The full set of dimension 5 operators has also been 
recently classified ||28| . Observational constraints on the various for different 
particle species can be obtained, for example, from ultra-high-energy cosmic rays, 
in particular by examining effects such as reaction threshold energies, and disper- 
sion or birefringence effects in long distance propagation. These constraints then 
allow to calculate the corresponding limits on the various Lorentz violating oper- 
ators. The lowest-dimensional operators which are not essentially ruled out at the 
Planck level by current observations, are CPT even operators of mass dimension five 
and six (see |83l[86i| and references therein). It can be shown that these correspond 
to a lowest-order allowed term of order p 4 in the above framework. 

An intuitive way to understand this is the following, (z) = 1,3 terms are es- 
sentially ruled out because they would lead to strong parity violations which are 
not observed. An (z) = 2 term would lead to Lorentz violation at all energies, 
which is also not observed. The first reasonable Lorentz symmetry breaking term 
is therefore (z) = 4. As a matter of fact, note that (z) = 1 terms would be ampli- 
fied by the prefactor Ei in the simple scheme of eq. ( |6.0.2 |, which should be ruled 
out even if it did not lead to parity violation. In a more refined scenario, (z) = 1 
would be suppressed, for example by writing fi 2 f^c 3 \p\/Mpi, with ]i the particle 
mass and Mpj the Planck mass. Even so, non-zero (z) = 1 terms are still so tightly 
constrained from observations as to be essentially ruled out ||85| . A similar remark 
applies to (z) = 2 terms. These are not strictly ruled out by observations when writ- 
ten in the form fi 2 f( 2 ^c 2 p 2 / Mpi, but still so tightly constrained that, for our present 
purposes, we can essentially forget about them. There exist several additional com- 
plications ||85ll86l . For example, dimension 3 and 4 Lorentz-violating operators are 
not strictly ruled out, because of the comment just given. They must be extremely 
small, but it is not excluded that there exists some theoretical scheme which 'natu- 
rally' accomodates this smallness. Another complication is that non-renormalisable 
operators of dimension 5 and higher induce lower-order terms (such as p 2 contribu- 
tions) via radiative corrections. CPT odd operators of dimension 5 are very tightly 
constrained, but again not completely ruled out, at least if one assumes the existence 
of some additional suppression mechanism (for example, a residual symmetry such 
as supersymmetry). Note that CPT even dimension 5 and 6 operators would still 
require some protection mechanism between the TeV scale and the Planck scale to 
suppress the lower-order radiative corrections. 
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For our purposes, we can simply conclude that all modifications of the disper- 
sion relation with (z) < 4 are essentially ruled out at the Planck level, or at least 
require some extraordinary fine-tuning, while even the (z) =4 terms can already be 
quite tightly constrained at the Planck level for the most common Lorentz-violating 
scenarios B53l |82~1 |83l . In any case, we can concentrate on (z) = 4 as the lowest- 
order allowed modification. Since we will focus on massless particles (photons and 
phonons), we can write a modified dispersion relation of the type 

E 2 = c 2 p 2 + f^cY/E 2 L/ (6.0.4) 

or, in terms of frequencies to and wave numbers k: 

cv 2 = c 2 k 2 + f^c 2 k 4 /k 2 L . (6.0.5) 



Several remarks might be useful. 

First, this type of Lorentz symmetry breaking implies a violation of the equiva- 



lence of reference frames. Indeed, as mentioned just above eq. ( |6.0.1[ ), a particular 
modified dispersion relation is valid in a particular observer's reference frame (or 
rather: a set of concordant frames). In this context, it should be mentioned that in the 
condensed matter models, it is natural to assume that the preferred reference frame 
is the laboratory frame, since this is the frame in which the privileged observer (the 



experimenter) measures things. In the laboratory frame, eq. (6.0.5) becomes 



(to - vk) 2 = c 2 k 2 + f^c 2 k A /k 2 L , (6.0.6) 



which immediately leads to the Bogoliubov dispersion relation for a BEC ( 2.5.3| > by 



setting ki = 2/£, with £ the healing length of the condensate, and where v repre- 
sents the flow velocity of the backgrounf fluid and c the co-moving speed of sound. 
In a relativistic context, the question of preferred reference frames is more subtle. 
For example, in the case of a black hole configuration, assuming that Lorentz in- 
variance is broken, there are at least two natural candidates for preferred reference 
frames: the rest frame of an observer free-falling into the black hole, and the rest 
frame of the black hole itself, which in the fluid analogy correspond essentially to 
an observer co-moving with the fluid and to the laboratory frame, respectively. We 
will see in the next chapter that this subtlety can play an important role. Note also 
that most of the work in the EFT framework uses the rest frame determined by the 
cosmic microwave background radiation to constrain the possible Lorentz-violating 
parameters. 



A second remark is that the previous equation (6.0.61 and the observation that 



ki = 2/£ illustrate an important lesson from condensed matter models, namely that 
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naive dimensional estimates in the domain of quantum gravity should be handled 
with care. Indeed, as we mentioned above, in quantum gravity phenomenology, it 
is usually assumed that Lorentz violation, if it occurs, should occur approximately 
at the Planck scale. However, condensed matter models teach us that the Lorentz 
violation scale Ei and the Planck scale Epi are not necessarily related. Indeed, in a 
BEC, the Lorentz violation scale Ei is determined by the healing length £, which de- 
pends ultimately on the s-wave scattering length. The analogue of the Planck scale 
Epi, on the other hand, can be thought of as the scale at which the granularity of 
the background becomes apparent. In this sense, then, the analogue Planck scale is 
determined by the interatomic distance flj, which depends essentially on the den- 
sity: cii ~ n -1 / 3 . It can then be seen intuitively from the weakness of the interatomic 
interaction in dilute gases that Ei <C E Pi for all BECs, but the important point is 
that Ei and Epi are in principle unrelated. This point might well extend to quantum 
gravity: different types of quantum gravity phenomenology could be characterised 
by different, mutually independent energy scales, which are not necessarily accessi- 
ble to an observer who is limited to the effective low-energy physics. Our previous 
mention of observational constraints on Planck-scale Lorentz violating operators in 
the EFT framework just before eq. ( |6.0.4[ ) could then be rephrased by saying that 
observations strongly support Ei ^> E P \. 

A third remark is that the sign of in eq. (6.0.4 1-(6.0.6 1 determines the prop- 
agation at high energies: < stands for subluminal propagation, i.e., high- 
energy modes propagate slower than low-energy ones, whereas the BEC model that 
we are interested in gives > and hence superluminal propagation: the high- 
energy modes propagate faster than the low-energy ones. 

A fourth and final remark is the following. Apart from the intrinsic interest of 
knowing whether Lorentz invariance is a fundamental symmetry of our universe or 
not, it is of course also to be hoped that phenomenological approaches to quantum 
gravity will allow us to discriminate between various theoretical frameworks for 
quantum gravity. Unfortunately, in these theoretical frameworks or 'top-down' ap- 
proaches for quantum gravity, it is currently not really clear whether Lorentz invari- 
ance is fundamental or effective, and in the latter case, how its breaking scale is re- 
lated to the Planck scale. For example, while many string theory scenarios axiomat- 
ically incorporate Lorentz invariance, it has been argued that in certain situations, 
violations of Lorentz invariance may occur in a way consistent with world-sheet 
conformal invariance [|87|, thus leading to acceptable string theory backgrounds. In 
the context of loop quantum gravity, in [54] it has been argued that quantum ef- 
fects should modify the relativistic dispersion relations, although the issue seems 
far from settled (see for some general remarks). It is arguably only in scenar- 
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ios of emergent gravity based on condensed matter analogies that the situation is 
clear: Lorentz invariance is a low-energy effective symmetry, and so it is certainly 
expected to break at some scale, although not necessarily related to (and therefore 
possibly much higher than) the Planck scale 11731 . 

To sum up, the possibility of Lorentz symmetry violations at high energy due to 
quantum gravity effects is an active branch of research, both theoretically and exper- 
imentally There are several frameworks to define such violations. The one that we 
discussed here is based on effective field theory. There are numerous examples of 
real physical systems, in particular condensed matter systems such as BECs, where 
this approach is the natural one to describe the effective low-energy Lorentz sym- 
metry and its high-energy breaking. This allows us to link studies of gravitational 
analogies in BECs, such as the ones presented so far in this thesis, directly to quan- 
tum gravity phenomenology with superluminal dispersion relations. It also means 
that we can take the complementary approach and study the effect of Lorentz vi- 
olation such as it occurs in BECs (and such as it might occur in quantum gravity) 
directly based on modified dispersion relations. We will use this second approach 
in the following chapter. 



Chapter 7 



Hawking radiation in a collapse 
scenario with superluminal dispersion 
relations 

In this chapter, we analyse the Hawking radiation process due to collapsing con- 
figurations in the presence of superluminal modifications of the dispersion rela- 
tion. With such superluminal dispersion relations, the horizon effectively becomes a 
frequency-dependent concept. In particular, at every moment of the collapse, there 
is a critical frequency above which no horizon is experienced. We show that, as 
a consequence, the late-time radiation originating from the collapse suffers strong 
modifications, both quantitative and qualitative, compared to the standard Hawk- 
ing picture. Concretely, we show that the radiation spectrum becomes dependent 
on the critical frequency, on the surface gravities associated with different frequen- 
cies, and on the measuring time. Even if the critical frequency is extremely high, 
important modifications can still show up. 

7.1 Introduction 

Hawking radiation is the process whereby a gravitational collapse which leads to 
the formation of a black hole is predicted to produce late-time thermal radiation due 
to the quantum creation of particles with an approximately Planckian (black-body) 
spectrum ||6"Tll62l . The calculation consists essentially in writing down Bogoliubov 
transformations to estimate the particle creation between inequivalent vacua, as we 
will detail further on. In 1976, Unruh developed a technique for replacing the col- 
lapse by a stationary spacetime representing an eternal black hole, by imposing ad- 
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equate boundary conditions on the past horizon (the so-called 'Unruh state')- This 
leads to exactly the same result as in the collapse scenario while alleviating some of 
its technical difficulties 11131 - 

Several alternative ways have since been developed to derive Hawking radia- 
tion (see [34] for a recent selective overview), for example through the stress-energy 
tensor trace anomaly ||37| or via a tunneling method II101H123H . Most of these deriva- 
tions are based on stationary black holes, and it is not always clear how they are 
physically related to Hawking's original collapse scenario. Nevertheless, it is re- 
markable that they all seem to lead to essentially the same result. However, one 
may wonder whether there is no hidden assumption in all these derivations, for 
example the implicit requirement of Lorentz invariance up to arbitrary frequencies. 
Such an implicit assumption about the physics at arbitrarily high frequencies was 
certainly the case in Hawking's original derivation, as was pointed out almost im- 
mediately after its publication, for example in [58j, and which has become known as 
the "transplanckian problem" with respect to Hawking radiationj^] It was also the 
main motivation for Unruh's 1981 paper ||114| , in which he pointed out that some 
of the high-frequency issues related to black holes can be studied by analogy with 
the propagation of sound waves in a fluid flow. In spite of the transplanckian prob- 
lem, Hawking radiation was and still is considered a crucial landmark in quantum 
field theory in curved spacetimes, and its implications for the connection (predicted 
earlier by Bekenstein fl24j) between black hole physics and thermodynamics are an 
essential touchstone for the development of quantum gravityn 




Here, we will stay quite close to Hawking's original derivation, and modify it 
carefully in a way that allows us to calculate the late-time black hole radiation for a 
collapsing scenario, but now in the case of superluminal dispersion relations at high 
frequencies. 

The motivation for such an analysis is threefold. First, as we indicated in the 

1 A related transplanckian problem exists in cosmology. 

2 It is perhaps useful to point out that Hawking radiation does not depend on the dynamical as- 
pects of gravitation, i.e., it is a purely kinematic effect, or in other words: a prediction of quantum 
field theory in curved spacetimes in the presence of black hole horizons. The thermodynamical as- 
pects of black hole physics, on the other hand, do crucially depend on the dynamical aspects of grav- 
ity, i.e., ultimately, on the Einstein equations. This is essentially because the first thermodynamic 
law of black holes, which establishes the relation between changes in its mass M and changes in 
the surface A of its horizon: dM = (k/8tc) dA (for a Schwarzschild black hole with surface gravity 
k), depends on the relation between the mass and the geometry of the black hole, and thus on the 
Einstein equations. A similar problem applies to the zeroth law, which establishes the constancy of 
the surface gravity over the horizon. Thus analogue black holes are expected to allow simulation of 
Hawking radiation, but not of their thermodynamic aspects IU17I I118I. 
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foregoing intermezzo |6j from the point of view of quantum gravity phenomenol- 
ogy, increasing attention has focused on the consideration that maybe Lorentz in- 
variance is not a fundamental law, but an effective low-energy symmetry which is 
broken at high energies. From extrapolation of current experiments, we know that 
there exist stringent bounds on the most commonly expected types of Lorentz viola- 
tions at the Planck scale. Nevertheless, even violations at much higher energy scales 
might still significantly affect black hole physics. On the other hand, it could also be 
that Hawking radiation is recovered precisely because of this relation between the 
(hypothetical) Lorentz violation scale and the Planck scale. In order to better under- 
stand this issue, a simple way of modelling a wide range of Lorentz violating effects 
(and a quite natural one from the point of view of condensed matter analogies) con- 
sists in modifying the dispersion relations at high energy. This modification can be 
subluminal or superluminal, depending on whether high-frequency modes move 
slower or faster than low-energy ones. 

This brings us to the (closely related) second motivation, namely the intrinsic 
question of the robustness of Hawking's prediction, regardless of whether the high- 
frequency theory describing our universe will indeed turn out to incorporate super- 
luminal corrections. Even if our universe turns out to be perfectly Lorentz invariant, 
this question is still interesting as a matter of principle, namely with regard to how 
inevitable or fundamental Hawking radiation really is, rather than an 'accidental' 
consequence of Lorentz invariance. The question of the robustness of Hawking's 
prediction with respect to modifications of the transplanckian physics has been tack- 
led principally within the same framework that we are using, namely by analysing 
effective field theories such that at high energies a modification of the dispersion re- 
lation is incorporated (see however [1J for a different take on the problem). Histori- 
cally, however, attention has mainly been given to subluminal dispersion relations. 
Important contributions such as EH EH E] E3 El HI5I . and more recently 
l|108t 11161 seem to have settled the robustness of Hawking radiation with respect to 
subluminal modifications of the frequency spectrum, at least in the case of station- 
ary black-hole scenarios, although this conclusion still rests on certain assumptions, 
usually related to the behaviour of the fields near the horizon. In any case, it is im- 
portant to remember that this only solves part of the transplanckian problem with 
regard to Hawking radiation. Indeed, subluminal modifications gradually dampen 
the influence of ultra-high frequencies, and so they do not explore arbitrarily large 
frequencies. So even assuming that it has been demonstrated that Hawking's result 
can be recovered in the stationary case in the presence of subluminal modifications, 
the question remains whether it is also robust with respect to a (non-damped) mod- 
ification of the transplanckian physics in a collapse scenario. Superluminal modi- 
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fications differ conceptually from subluminal ones, in that they gradually magnify 
the influence of ultra-high energies, and thereby offer an interesting test-case for the 
transplanckian robustness of Hawking radiation. 

A third motivation for our study, as already mentioned several times throughout 
this thesis, comes from a possible connection with experiment. BECs are expected to 
be good candidates for a possible experimental verification of Hawking radiation [8; 
m |35t |H |56j |57l ESS S^Zl IHa- Furthermore, we have shown that dynamically 
stable black hole configurations are possible in a BEC, thereby removing a major 
potential complication at the moment of carrying out such an experiment. But since 
the physics of BECs automatically leads to a superluminal dispersion relation at 
high energies, the question is again which kind of modifications are to be expected in 
the laboratory realisation of a BEC black hole with respect to the standard Hawking 
picture. 

Of the above-mentioned works on the robustness of Hawking radiation, a few 
have also tried to address superluminal modifications ||4Tt [108 : 116]. Various prob- 
lems make this case quite different from the subluminal one. These problems can be 
related to the fact that the horizon becomes frequency-dependent when modifying 
the dispersion relations. This is also true in the subluminal case, though qualita- 
tively in a very different way. With superluminal modifications, the horizon lies 
ever closer to the singularity for increasing frequencies, and asymptotically coin- 
cides with it. This causes the 'apparent' interior of the black hole (the interior of 
the zero-frequency horizon) to be exposed to the outside world. Since it seems un- 
reasonable to impose a condition arbitrarily close to the singularity as long as we 
do not have a solidly confirmed quantum theory of gravity, most of the approaches 
used for subluminal dispersion are invalid, or at least questionable, for superlumi- 
nal dispersion. Moreover, since it seems reasonable to expect that quantum effects 
will remove the general relativistic singularity, an overall critical frequency might 
appear above which no horizon would be experienced at all. 

In the analysis that follows we will try to avoid making any further assumptions 
about the physics near the singularity. For instance we will analyse the characteris- 
tics of the radiation at retarded times at which the singularity has not yet formed, 
and consider both final options, either with or without a singularity. Our approach 
will be based on a derivation of the radiation very much in the spirit of Hawking's 
original calculation through the relation between the asymptotic past ray trajecto- 
ries and the asymptotic future ones in the case of a collapsing configuration. The 
only assumptions about these asymptotic extremities are the standard ones, namely 
a Minkowski geometry in the asymptotic past, and flatness at spatial infinity also in 
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the asymptotic future. The language we will use is related to the fluid analogy for 
black holes, which provides the intuitive picture of the spacetime vacuum flowing 
into the black hole and getting shredded in its centre. 

Our main results can be summarised as follows. Three crucial elements distin- 
guish the late-time radiation with superluminal dispersion relations from standard 
relativistic ones. First, at any instant there will be a critical frequency above which 
no horizon has yet been experienced. This critical frequency will induce a finite limit 
in the modes contributing to the radiation, which will therefore have a lower inten- 
sity than in the standard case, even if the critical frequency is well above the Planck 
scale. Second, due to the effective frequency-dependence of the horizon, the sur- 
face gravity will also become frequency-dependent and the radiation will depend 
on the physics inside the black hole. Unless special conditions are imposed on the 
profile to ensure that the surface gravity is nonetheless the same for all frequencies 
below the critical one, the radiation spectrum can also undergo a strong qualita- 
tive modification. Depending on the relation between the critical frequency and the 
Lorentz violation scale, the radiation from high frequencies is no longer negligible 
compared to the low-frequency thermal part, but can even become dominant. This 
effect becomes more important with increasing critical frequency. Finally, a third 
effect is that the radiation originating from the collapse process will extinguish as 
time advances. 



The remainder of this chapter has the following structure. In section 7.2 we will 
describe and motivate the classical geometry of our model and the concrete kinds of 
profiles that we are interested in. Note that these configurations are related, but not 



identical to the ones discussed in the chapters [3]- [5] In section 7.3 we will briefly re- 
view how Hawking's standard result can be obtained in the formalism of our choice 
for the case of standard relativistic dispersion relations. This calculation will be 



adapted in section 7.4 to superluminal dispersion relations, and we will analyti- 



cally obtain a formula for the late-time radiation for this case. Then, in section 7.5 
we will present graphics obtained from this analytic formula by numerical integra- 
tion. These graphics will illustrate the results mentioned in the previous paragraph, 
which we will discuss in more depth and compare with other results in the recent 



literature in section 7.6 



7.2 Classical geometry 

We will study radiation effects in simple (l+l)-dimensional collapsing configura- 
tions. It is well known that in the Hawking process, most of the radiation is pro- 
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duced in the s-wave sector. Hence, a spherically symmetric treatment, effectively 
(l+l)-dimensional, suffices to capture the most relevant aspects of the process. We 
will work in Painleve-Gullstrand coordinates, where 

ds 2 = -[c 2 -v 2 (t,x)]dt 2 -2v(t,x)dtdx + dx 2 , (7.2.1) 

and make the further simplification of taking a constant c (we will always write 
c explicitly to differentiate it from the frequency /wave-number dependent that 
will show up in the presence of a superluminal dispersion relation). In this man- 
ner all the information about the configuration is encoded in the single function 
v(t,x). The Painleve-Gullstrand coordinates have the advantage, compared to the 
Schwarzschild form, of being regular at the horizon. Moreover, they suggest a nat- 
ural interpretation in terms of the language of acoustic models. In such a model, as 
before and throughout this thesis, c represents the speed of sound and v is the ve- 
locity of the fluid (which corresponds to the velocity of an observer free-falling into 
the black hole). Also as before, we will consider the fluid to be left-moving, v < 0, 
so that the outgoing particles of light/ sound move towards the right. 

The most relevant aspects of the analysis presented hereafter depend only on 
the qualitative features of the fluid profile. However, to justify some specific calcu- 



lations later on it is helpful to use the following concrete profiles, see Figs. 7.1 and 



7.2 Consider a velocity profile v{x) such that v(x — > +oo) = 0, v(x = 0) = — c and 
further decreasing monotonically as x — > — oo. In the fluid image, the fluid nearly 
stands still at large distances and accelerates inwards, with a sonic point or horizon 
at x = 0. The further decrease in v can either be linear until a constant limiting value 



is almost achieved, as in Fig. 7.1 or of the form 



, - 2M/c 2 

°« = -^/^T2m7?' {722) 



again up to a constant limiting value, as in Fig. 7.2 Let us recall that this second 
velocity profile corresponds to the Schwarzschild line element with M the central 
mass, as can be seen by reparameterising the time coordinate and using r = x + 2M/ c 2 
as radial coordinate ||117| . 

Up to here, we have described a stationary profile. Now, to incorporate the dy- 
namics of the collapse, let us introduce a monotonically decreasing function £(f), 
£(f — ► — oo) — > +oo, and define v(t,x) as 



(7.2.3) 
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Figure 7.1: Velocity profile of a black hole with a linear slope of the velocity v(x), 
and hence a constant surface gravity k w /, from the classical (zero-frequency) hori- 
zon at x = down to some predefined limiting value. The auxiliary function 
£(f) separates each instantaneous velocity profile v{t,x) into a dynamical region I: 
v(t,x) = v(t) = v an d a stationary region II: z>(t,x) = z;(x). 




Figure 7.2: Velocity profile characteristic of a Schwarzschild black hole. The slope of 
v(x) increases from the classical horizon at x = leftwards up to a constant limiting 
value, leading to a surface gravity k w i which increases with the frequency. 




Figure 7.3: [From top to bottom and left to right] Various snapshots of the dynamical 



evolution of a collapsing configuration with variable surface gravity, as in Fig. 7.2 
The classical horizon is formed at £4. 



Imagining the collapse of a homogeneous star, the function £(£) represents the dis- 
tance from the star surface to its gravitational (Schwarzschild) radius. The dynami- 
cal configuration that we obtain consists of a series of snapshots, see Fig. 7.3 In each 
snapshot, |p(f,x)| increases (i.e., the fluid accelerates) from v = for x — > +00, up 
to a point xq = £(£), and then remains constant as x further decreases towards —00. 
For consecutive snapshots, the point Xo = £(0 moves leftwards, so v(t, x) covers an 
ever larger part of v{x). 

From the point of view of an outgoing (right-moving) particle of light/ sound, 
the configuration is nicely split up into, first, a dynamical or £-dependent region, 
and then a stationary x-dependent region, as defined in eq. ( |7.2.3| >j^] 

The difference between the two types of profiles, Figs. 7.1 and |7.2[ can best be 

3 In a physically realistic model, the apparent kink in the profile where the transition between 
both regions takes place will of course be smoothed out. 
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explained in terms of the surface gravity. When modifying the dispersion relation, 
the horizon becomes a frequency-dependent concept: each frequency experiences 
a different horizon, as we will see explicitly in section 7.4[ In particular, for super- 
luminal modifications, the horizon forms later (i.e., at higher values of \v\, or more 
negative values of x) for increasing frequencies. The surface gravity will then also 
become frequency-dependent: 



HO' 



dv 
dx 



(7.2.4) 



where % w i is the frequency-dependent location of the horizon. This surface gravity 



will be seen in section |7l4| and 7.5 to play a crucial role. This explains our choice of 



the two types of profiles: In the first profile, Fig. 7.1 we have considered a linear 
velocity profile \v(x)\ from the classical or zero-frequency horizon (c = \v\) down to 
some given limiting value, |u| max , from where it stays constant. This maximum ve- 



locity defines a critical frequency co' c as we will discuss in more detail in section 7.4 



Then the surface gravity k^i will be frequency-independent up to this critical fre- 
quency co' c after which it will rapidly vanish. In the second profile, see Fig. 7.2 we 



have taken a velocity profile typical for a Schwarzschild black hole, and therefore 
k w i increases with a/, again up to a predefined limiting value corresponding to the 
horizon for a critical frequency to' c . 

In the context of the fluid analogy, it seems obvious that some mechanism will 
avoid the formation of a singularity, while in gravity it is also usually expected that 
quantum effects will ultimately resolve the general relativistic singularity at the cen- 
tre of a black hole. But in any case, we will also take into account the limit for an 
infinite critical frequency, which corresponds to the eventual formation of a singu- 
larity, and is characterised by a velocity profile with |u| max — > oo. 



7.3 Standard Dispersion Relations 

In this section we will briefly review a way of deriving Hawking's formula for the 
radiation of a black hole with standard (relativistic) dispersion relations of the form 
co 2 = c 2 k 2 . For the sake of simplicity, we will only consider a massless scalar field. 
We will summarise the main steps: we start with a sketch of the general idea and 
then go from the Klein-Gordon inner product and the definition of the Bogoliubov /3 
coefficients, over the relation between past and future null coordinates, to the black- 
body radiation in the wave-packet formulation. Our aim is to present the key points 
of the procedure in such a way that they can easily be adapted to superluminally 
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modified dispersion relations — the subject of the next sections. 
7.3.1 General idea 

Hawking radiation can be understood as an example of mode mixing between pos- 
itive and negative-energy modes in a dynamical curved spacetime l|26ll5Ttl62| . 

In flat spacetime, due to the high amount of symmetry present, all inertial ob- 
servers agree on the definition of the vacuum and of what constitutes a particle. 
But this is no longer the case in curved spacetime. In particular, the definition of 
positive-energy modes with respect to some choice of time coordinate t is that 
they should be eigenfunctions of the operator d t such that 

d tt p(+) = -icvtp( +) for co > . (7.3.1) 

This definition assumes the existence of a timelike Killing vector field associated 
with t. In Minkowski spacetime, the Poincare group guarantees that all inertial 
observers agree on the definition of positive-energy modes. However, in curved 
spacetime the Poincare group is no longer a symmetry of the spacetime, and so the 
previous conclusion does not hold in general. We will therefore compare the vacua 
defined by two different observers. 

Assume two complete orthonormal sets of modes {i/?-} and {ipj}, and express 
the scalar field (p(x) in terms of each of them: 

<p{x) = Lihfiix) + ajtp'^x)} = Lj[bjfj(x) + b]fj(x)) . (7.3.2) 

Each decomposition defines a corresponding vacuum, 1 0/) and |0jj) respectively, 
through 

OflOj) =0 (Vz) (7.3.3) 

bj\0 n ) =0 (V/) (7.3.4) 

Since both sets are complete and orthonormal, one can write: 

fi = EM! + htf) > b i = E(4- fl i + ^4) > (7-3.5) 

i i 

where = {ipj,ip[) and jSyj = —(ipi,ip'*) are the so-called Bogoliubov coefficients. 

If 7^ then there will be mode mixing: a mode expressed in terms of tpj will 
necessarily contain contributions from modes in terms of ip'*. This implies 

&y|0f> = E^Kl°Jr>^°- (73 - 6) 
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The vacuum |0j) associated with {i/^} is then not a vacuum state for an observer 
associated with {tpj}. One can calculate the amount of particles that such an ob- 
server will see in the Vacuum' state \0j) as follows. From the number operator 
Nj = b~jbj, and assuming that all one-particle states are properly normalised, one 
obtains (0j|N ; -|0i) = Ei \Pji\ 2 °- The 

vacuum |0j) associated with {ip'^j therefore 
contains a number yj, \fiji\ 2 of /-particles in {ipj}. 

When indexing the modes with a continuous index such as the frequency, as we 
will do in the remainder of this chapter, instead of the discrete indices i and /, this 
sum becomes an integral. But similar technical complications left aside, the brief 
discussion given here captures the essence of particle creation in curved spacetimes. 

The particular case of Hawking radiation is obtained by associating {i/^} to an 
asymptotic past observer, before the onset of the collapse, and {xpj} to an asymptotic 
future observer. For both observers, 9f is (asymptotically) a timelike Killing vector. 
The observer in the past can then define the ip'i as the positive-energy modes and ip'* 
as the negative-energy ones with respect to t, and likewise for the future observer. 
Even though, for both observers, the spacetime is asymptotically stationary, the dy- 
namical evolution of the spacetime inbetween means that they will not agree on the 
definition of the vacuum. In particular, there will be particle creation as observed 
by the second (asymptotic future) observer originating in the quantum state which 
the first (asymptotic past) observer would call the vacuum state. 

Note that the definition of the Bogoliubov coefficients assumes that there is a 
well-defined inner product, which we will discuss now. 



7.3.2 Inner product 

The d'Alembertian or wave equation for a massless scalar field in 1+1 dimensions 



for the metric (7.2.1 ) with constant c can be written as 

(3 f + d x v) (df + vd x )cf> = c 2 d 2 x cp . (7.3.7) 

This conformally invariant theory is equivalent to the dimensionally-reduced 3+1 
spherically-symmetric theory if one neglects the backscattering due to the angular- 
momentum potential barrier (responsible for the so-called grey-body factors). 

In the space of solutions of this wave equation, we can define the Klein-Gordon 
pseudo-scalar product 



(7.3.8) 
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which is independent of the choice of the spatial slice E. For a t = constant slice, 
and in particular for t — > +00, this becomes 

(<PV<P2) = ~i J dx [<Pl@t + vd x )cp2 - (pl(d t + vd x )cp 1 ] . (7.3.9) 
We can define future null coordinates u(t,x) and w(t,x), such that, when t,x — > +00, 

u(t,x) — ► t — x/c , w(t,x) — ► t + x/c . (7.3.10) 
Writing the inner product in terms of these null coordinates, gives, for the limit 

ic ( f +0 ° 

(<Pl><P2)= ~2 [./-co ^ ^ ldu<P ^ ~ <P2 du( Pl\w=+oo 
r+00 ~\ 
+ J co ^ ][CpldwCp * 1 ~ ^2 a i^l] M =+oo f ■ (7.3.11) 

Note that, in the derivation of this formula, we have only made use of coordinate 
transformations, without relying on their null character or on geometrical tools such 
as conformal diagrams or the deformation of Cauchy surfaces. This is important 
because in the case of modified dispersion relations, such geometrical concepts be- 
come problematic and actually can only be maintained in the context of a rainbow 
geometry — if at all. 

Similarly we could have used past null coordinates (U, W), which obey, when 

t — > —00, 

U(t,x) -> t-x/c, W(t,x) -> * + x/c , (7.3.12) 
to calculate the inner product in the asymptotic past. 



7.3.3 Bogoliubov /3 coefficient 

The right-moving positive-energy solutions associated with the asymptotic past and 
with the asymptotic future, normalised in the Dirac-delta sense, can be expressed as 

= ^=L=e-™'u , ^ = ~J=e- i(VU , (7-3.13) 
y2nc co 1 yznc co 

respectively, where we use primes to indicate asymptotic past values. 

The mode mixing relevant for the Hawking process occurs in the right-moving 



sector. Therefore, we only need to calculate the first term in the scalar product ( 7.3.11 > 
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Then, plugging eq. ( |7.3.13| > into the definition of /3 given just below eq. ( |7.3.5| > and 



integrating by parts we obtain the simple expression 

/W = ^{^i j du e- lJu ^e~ wu , (7.3.14) 

so that all the information about the produced radiation is contained in the relation 

U = U(u) = U(u,w — > +oo). 

7.3.4 Relation U(u) 

For a standard relativistic dispersion relation, it is well known (see, e.g., [26]) that 
the relation between U and u for a configuration that forms a horizon (in our case 
at x = 0) can be expressed at late times as U = Uh — Ae~ KU ^ c (note that we use 
a subscript H for all quantities associated with the horizon), where Uh, A and the 
surface gravity 



dv 



dx 



(7.3.15) 

x=0 



are constants. 



We can define a threshold time uj at which an asymptotic observer will start to 
detect thermal radiation from the black hole. This retarded time corresponds to the 
moment at which the function U(u) enters the exponential regime. We can then 
rewrite the previous expression, valid for u > U[,as 

U=U H - A e" K( """ l)/c . (7.3.16) 



Plugging this relation into eq. (7.3.14 1 and integrating in u gives 



Puts = ^exp[-za/(i H ]exp[-^ln(a/A )] 

x exp(— icouj) exp(— C )T(ico)/%) . (7.3.17) 

7.3.5 Wave packet formulation 

In order to obtain physically sensible results, it is a good precautionary measure to 

replace the monochromatic rays used so far by wave packets (see e.g. Il62l or the 
discussion in BSTl ). Positive-energy wave packets can be defined as 

f 4= -lAco < co - to; < lAco /r , „ 

F W Ul {co) e Va, 2 12 (7318) 

I otherwise , 



86 



CHAPTER 7. HAWKING RADIATION IN A COLLAPSE SCENARIO 



where U\ = uq + 2nl/ Aco with uq an overall reference and I an integer phase pa- 
rameter. The central frequencies of the wave packets are C0j = ]Aco, with Aco their 
width. 

Then, from the expression fi W j,ui;u' = / PuHv'PtOjW (^) / and assuming that the 
wave packets are sufficiently narrow (Aa; <C coj), we obtain 

2 _ cAco sin 2 (z-zQ 1 
l/WI ~ 2 ncv'K (z-z,) 2 exp( ^)_ 1 ' (7 - 3 - 19) 

where we have defined 

z = C -^- In co'Aq ; z x = ^ («/ - Uj) . (7.3.20) 

Finally, integration in co' gives the number of particles with frequency coj detected 
at time u\ by an asymptotic observer: 

"+0O 



j Q dco' |/3 Wy , ii;;a/ | 2 (7-3.21) 
1 (7.3.22) 



exp(2ncco j / k) — 1 ' 

which reproduces Hawking's formula (in the absence of backscattering) and corre- 
sponds to a Planckian or black-body spectrum with temperature Tjj = k/ (2tzc). 

Note that the essence of the transplanckian problem is contained in this last step, 
since the integration must indeed be carried out up to infinite frequency to recover 



the standard Hawking result, an issue which we will further discuss in section 7.4.5 



7.4 Superluminally modified dispersion relations 

In this section we will indicate how the late-time radiation originating from the for- 
mation of a black hole through a collapse process can be calculated in the case of 
superluminal dispersion relations. 

We have discussed earlier in this thesis how superluminal dispersion relations 
can arise as a phenomenological description in quantum gravity scenarios with 
Lorentz violations at high energies, see chapter |6j and how they are also naturally 
found in BECs (chapter [2]). From a theoretical point of view, one can also introduce 
superluminally modified dispersion relations by adding a quartic term to the wave 
equation: 

(3, + d x v) (d t + vd x )<p = c 2 (d 2 x + -U 4 T ] <p , (7.4.1) 
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Figure 7.4: Behaviour of the effective phase C/ cp ^ and group speeds of 
light/ sound with respect to the wave number k. Due to the non-relativistic disper- 
sion relation, the effective velocities become 'superluminal' for k > ki (where k^ is 
the Lorentz violation scale). 



where ki (the Lorentz symmetry breaking scale) is the scale at which non-relativistic 
deviations in the associated dispersion relation 

(a> - vkf = c 2 k 2 (l + ^ j (7.4.2) 

become significant. The reasons for using this relation in particular were discussed 
in chapter |6j Let us just recall that in Bose-Einstein condensates, ki = 2/£, with £ 
the healing length of the condensate. We also repeat that in quantum gravity, this 
Lorentz violation scale is usually identified with the Planck scale, although analogue 
gravity models indicate that the two are not necessarily related. In any case, quali- 
tatively, our results will not depend on the specific form of the deviations from the 
relativistic dispersion relation but on their superluminal character. 

This dispersion relation leads to a modification in both the phase velocity 
and the group velocity Vg. For a right-moving wave, we have 

w do; ,„ A „. 

v P h = j = Ck, P h + v, v g = — = c k/g + v , (7.4.3) 

where we have introduced the effective fc-dependent phase and group speeds of 
light/ sound 

/ W 1+2 F 

Ck,ph = c \l 1 + J2 ' C ^ = C ^=^' ( 7A4 ) 

1 + l 



respectively (see Fig. 7.4). Both and c^g become larger than c (i.e., become 
'superluminal') as k increases above k^. 

At first sight, it seems obvious that the ray equation should be defined in terms 
of the group velocity Vg. Nevertheless, the question of whether the velocity relevant 
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for Hawking radiation is the group or the phase velocity seems to be tricky [108: 121 ; 
I122| . For example, (7.4.9) below suggests that the phase velocity might be relevant. 
We limit ourselves to remark that C/ c „ and show the same qualitative behaviour. 
Hence our results are independent of this issue and we will simply write (or 
cu w ') when wishing to emphasise the frequency-dependence). There will then be 
a horizon, which becomes frequency-dependent, when + v = 0, irrespectively 
of whether „ or „^ is used for c^. Moreover, since C/ c becomes arbitrarily high 
for increasing wave number k, given a certain |u| m ax at a particular instant of time, 
there will be a critical co' c such that waves with an initial frequency to' > co' c do not 
experience a horizon at all. The only exception to this rule occurs when the velocity 
profile ends in a singularity v — > — oo, which implies to' c — > oo. 

Our aim is to calculate the black hole radiation with superluminally modified 
dispersion relations. We will now repeat the main steps of section [73} and point out 
where and how these modifications must be taken into account. 



7.4.1 Generalisation of inner product 



The essential point with regard to the pseudo-scalar product ( |7.3.9 ) is the following. 



Its explicit form for t =constant is not changed by the presence of the d\ term in the 
wave equation (see also the discussion in [42]). Indeed, it is still a well-defined inner 
product, and in particular it is still a conserved quantity, since 

dt((p\,(pi) = J dx <pi d\ <p\ , (7.4.5) 

which can be seen, by repeated integration by parts, to vanish under the usual as- 
sumption that the fields die off asymptotically. Note that the modification of the dis- 
persion relation singles out a preferred time frame: the 'laboratory' time t. Chang- 
ing to another time t will in general lead to a mixing between t and x, and hence the 
simple relations given here would no longer be valid. 

Using the preferred time t, and making exactly the same change of coordinates 
as in the case of standard dispersion relations, we can again transform the inner 
product, evaluated at t — > +oo, into the expression ( 7.3.11| > in terms of u and w. As 



in the standard case, only the first term is relevant for the Hawking process: 



ic 
2 



r+oo 

/ du [<pid u <p* 2 - <P2d«<PlL=+co • ( 7A6 ) 
J — oo 



Note that we are now using u and w merely as a perfectly good set of (auxil- 
iary) coordinates, in order to cast the inner product into a useful form. However, 
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here they do not have the null character they had in the case of standard dispersion 
relations. 



7.4.2 Rainbow null coordinates 

Let us define some sets of spacetime functions that will prove to be useful in what 



follows. Given a fluid profile of the type described in section 7.2 we can integrate 
the ray equation 

dx 

— = c k(at) (t,x)+v(t,x) t (7A.7) 

starting from the past left infinity towards the right. The ray has an initial frequency 
to' and an associated initial wave number k', from which we can deduce the value of 
C/ c / . In the left region, where the velocity profile is dynamic but position-independent 



(v(x, t) = v(^(t)), with £(f) the auxiliary function introduced in section 7.2), k = k' 
can be considered as fixed while the frequency changes (this is what happens in a 
mode solution of equation ( 7.4.1] > in this region). Then, we can define the function 



U w <{t,x) = to{t)dt-^x, (7.4.8) 

where to{t) is the instantaneous frequency of the particle at each time f, defined 
through the dispersion relation and such that Q{t — > — oo) = to'. 

When the ray reaches the kink it passes into a stationary region in which the 
velocity profile only depends on the position (v(t, x) = v(x)). At the kink, the ray 
still has the initial wave number k' and a frequency oj. In its propagation towards the 
right this frequency now remains fixed while its wave number becomes a decreasing 
function k(x) of the position, such that fc(£(fo)) = k! , with to the moment at which 
the kink is crossed. The final frequency of the ray will then simply be to and its final 
wave number k = lim x ^ +00 k{x). In this region, then, the function U w /(t,x) can be 
expressed as 

U w >{t,x) =t- I ^-dx . (7.4.9) 

J CO 

Note here that to and k(x) both depend on the initial to'. 

The same can be done by integrating the ray equation starting from the future. 
In this way we can define u w functions. (The same procedure can be used to define 
W w ! and w w ). It is worth noting that U w i and are not null coordinates in the 
usual geometric sense, since they are frequency-dependent, but they nevertheless 
share many properties with null coordinates. 
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7 A3 Bogoliubov jS coefficient 



The calculation of the inner product (7.4.6) involves the limit w — > +oo (equiva 



lently, Wa> — > +oo). We can therefore change variables in the inner product from u to 
Ucv = Uio(u) = Uu){u,zv —¥ +oo). The combination of the derivative and the integral 



means that the form of (7.4.6) is preserved in the new integration variable u w . We 



then can write the Bogoliubov jS coefficients relevant for the Hawking process as 

ic /'+ 00 

/W = -y / dUallptfduulpa - ipujdu^U^+co . (7.4.10) 

Now, assuming that the profiles vary slowly (in scales much larger than the 
Planck distance), the right-moving positive-energy modes associated with past and 
future infinity can be approximated by the following simple expressions: 

ip'ut ~ 1 e-MUj f (7.4.11a) 

V27T C CO 1 

xp w n 1 e~ i(VUw , (7.4.11b) 

V27T C CO 

so that the Bogoliubov coefficients read 

Puicv' ~ ^\f^7 [dUu g-MUjMg-iuUu (7 4 12 ) 

2n V co 'J 

In analogy with the standard case, all information about the radiation is seen to 
be encoded in the relation U^^Ua,). In this expression co' is the initial frequency of 
a ray at the past left infinity and co = to (to') its final frequency when reaching the 
future right infinity. 

The approximation just discussed for the Bogoliubov coefficients amounts to 
considering profiles that vary sufficiently slowly both with x and t. This is equiva- 
lent to considering large black holes. In general the quartic term in the wave equa- 
tion, or equivalently the quartic modification of the dispersion relation, introduces 
a new source of backscattering, on top of the usual angular-momentum potential 
barrier which we have already neglected. In our approximation, this additional 
backscattering (beyond the standard grey-body factors) has been neglected. For 
large black holes this contribution will in any case be very small as has been ob- 
served in numerical simulations [108]. In addition we are also neglecting any re- 
flection caused by the kink. However, let us remark that given an approximate 
scheme for calculating W w /(m w ) for general profiles v(t,x), the same U w i(uco) ob- 
tained from a profile with a kink would be obtainable from one (or several) specific 
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v(t,x), this time perfectly smooth and thus causing no further backscattering. Our 
results, which rely only on the specific form of the relation U^iu-aj), are therefore 
valid beyond the specific configurations with a kink presented in this chapter. 



7.4.4 Relation 24/ (i/o,) 

Our next task is to calculate the uniparametric family of functions U w i and the rela- 



tion between lA^i and u w for different configurations. As explained in section 7.4.2 



the relation W^(m w ) is obtained by integrating the ray equation dx/dt = c\ + v us- 



ing the profiles discussed in section 7.2 These can be described by means of sta- 



tionary profiles v(x) and an auxiliary function £(£), see the definition of v(t,x) in 



eq. ( |7.2.3[ ). We will use a straightforward extension of the procedure established 
in §F7\ (see also [18j for a summary) for a relativistic dispersion relation. Care must 
be taken, however, with the quantities that depend on the frequency co'. In particu- 
lar: 

• We will denote by x^ /CO i and f# ^ the position and the time at which the hori- 
zon associated with a particular initial frequency co' is formed. 



• The surface gravity k w /, defined in ( 7.2.4[ >, allows to write, for all to' < co' c and 
for x close to x^y: 

1 

v(x) w -c k + -k uI {x - x H/(V> ) , (7.4.13) 

up to higher-order terms in x — Xh,o/- 

• In a similar way, £ (f ) can be linearised near f# by introducing an a/-dependent 
parameter X w r. 

£(t) ~ XHrf ~ - tH,cv') , (7-4.14) 

again up to higher-order terms. 

Let us consider laboratory times t > tfj,o = t Hr(V > =0 such that £(f) has already 
crossed the classical horizon at x = 0: £(£ > tH,o) < 0. As we explained earlier 
we can define a critical frequency co' c as the minimum initial frequency such that 
Cfc + v(£(t)) > for all co' > co' c , or, in other words, the minimum frequency which 
at that particular time has not experienced any horizon yet. 

We first integrate the ray equation in the dynamical part of the profile, where 
v = v(£(t)), for rays crossing the kink just before the formation of the horizon. The 
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integral in this region runs between the asymptotic past and the point (io,w'/ x o,o;') 
at which the kink separating the dynamical and the stationary regions is crossed for 
this particular frequency: 



■0,w' 



dt 



^(£0)) + c fc(a/) 



X 0,u>' 



dx . 



(7.4.15) 



Writing c k i = lim^-oo cu w >) for the speed of propagation in the asymptotic past 
(and by extension in this whole first region, since k = k' remains constant, see sec- 
tion 7.4.2), using the property U w i{t,x) — > t — x/cy for t,x — > — oo, and noting that 



x o,cv> = this leads to 

which can be written as 



+ 



0,w' 



C k > J-co 



c 0,a/ — l H,a>' 



(7.4.16) 



(7.4.17) 



for small values of \to /CO / — t HA ,r\, where 



U 



H,co' = l H,cv' 



+ 



i r'H, 



(7.4.18) 



c k' c k' J -oo 

is the ray constituting the horizon associated with the frequency a/. 

In the stationary part of the profile, where v = v{x), we can proceed in a similar 
way and write 



lim t 0a}/ 

Xt;tf—*+00 1 



+ 



dx 



Xf 

c k,f ^(W) v(x)+c k(c 



(7.4.19) 



where c k j = lim^ +00 Cu w i\ is the speed of propagation in the asymptotic future. 
Using the identity 

£(*<W') x f _ dx 



c Kf c Kf 



Ck,f 



(7.4.20) 



eq. (7.4.19 1 can be cast in the following form: 

"+00 



U L0 — t 0,LO l T 



c hf ^(W) 



dx 



+ C k (co>) C k/f 



(7.4.21) 



We now focus on the integral on the right-hand side and split the integration interval 
into three parts: 



+ 00 



Xj rx u r+oo 

+ + / 



(7.4.22) 
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with xj close to £(£o,a/) Xji sufficiently large so that, effectively, v(x) — > and 
c fc(a/) ~^ c fc,/ f° r * — The value of the third integral now approximately van- 
ishes. The value of the first integral evaluated at its upper integration boundary 
%\, together with the complete second integral, can be absorbed into a largely irrele- 
vant bulk integration constant C\. What remains is the behaviour of the first integral 
around its lower integration boundary £(fo /W ')- We are interested in its behaviour 
for small values of | f ,o/ 
obtain 



t H/lv i | . Using ( |7.4.13| > to approximate the integrand, we 



t 



0,u>' 



at 



0,cv' 



Ck,f 



1 

c Kf 



-Ax , 



(7.4.23) 



from which we conclude that 



x H,o 



-Kjllco/C 



(7.4.24) 



Both regions are connected by combining eqs. (7.4.17) and ( 7.4.24| >. Making use 



of eq. ( |7.4.14| >, this gives 



U„, = U 



H,io> 



A w , e - K "' u " /c 



(7.4.25) 



This relation is valid for frequencies for which a horizon is experienced, i.e. for 
to' < co' c and times u w > Ui fU ,r, where U\ iW i is the threshold time defined in sec- 
tion 7.3.4 which, unsurprisingly, has become frequency-dependent. Again, as in the 
case of standard dispersion relations — section 7.3.4 — we can write 

A e-^ /{u 



a'' 



u 



H,oj> 



l,a>' 



)/c 



(7.4.26) 



where Ui^i is essentially (with possible higher-order corrections). Assuming 
that the collapse takes place rapidly, it is a good approximation to replace by 
%,o/ • Actually, as we will see shortly, this is a conservative estimate, in the sense that 
it slightly underestimates the superluminal correction to the radiation spectrum. 



7.4.5 Wave packet formulation 

The relation W w /(m w ) we are considering, interpolates between a linear behaviour at 
early times and an exponential behaviour at late times. It has the same form as the 
relation U(u) for standard dispersion relations, and so we can continue following 
the steps of the standard case. In particular, the equivalent of eq. ( |7.3.17[ ) is obtained 
by integrating out u w . Note that we use the subscript co to emphasise that the u w 
are not the null coordinates of the standard case, but this should not be interpreted 
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as an explicit function of co and so does not complicate the integration steps in u w 
and co. However, when integrating \Pu>j,ui;cv' 1 2 i n ^ to obtain N w . /Ul , see section 7.3.5| 
we must carefully consider the frequency-dependence of the relevant terms, i.e., of 
Uh,cv'' Ui i0J i and k w i. The term carrying Uh /C0 / is moduloed away in eq. ( |7.3.19| ), and 
we have replaced U\ i(jJ i by u IfCJ i , so the only relevant frequency-dependent factor that 
we are left with is the surface gravity k w i . 

Moreover, because of the critical frequency co' c in the horizon formation process, 
a finite upper boundary will also be induced in the integral. Indeed, frequencies 
co' > oo' c do not contribute to the radiation at all, since they do not experience a 
horizon. This is a delicate but crucial point. It was already observed long ago by 
Jacobson ||65l that trying to solve the transplanckian problem naively by imposing 
a cut-off frequency would seemingly extinguish Hawking radiation on a relatively 
short time scale. In our case, however, this cut-off is not imposed ad hoc, but appears 
explicitly because of the superluminal character of the system at high frequencies. 
Moreover, the critical frequency, and hence the upper boundary induced in the in- 
tegral, depend directly on the physics inside the horizon. Indeed, given a certain 
velocity profile, and in particular its behaviour near the centre of the black hole, the 



critical frequency can be calculated by setting — \v\ in eq. ( |7.4.4[ ) and extracting 
the corresponding critical frequency from the dispersion relation ( 7.4.2| >. We will see 



this effect graphically in section 7.5 



In analogy with eqs. (7.3.19) and (7.3.21 1, we now obtain the number of particles 
detected for each frequency C0j as 



LOj,Uj 



J dco \fi Wj ,ui* 



(7A.27) 
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t<-uj'=o (which corresponds to the standard k 
Uj^i). Changing the integration variable from co' 



to z, we finally obtain the central expression in our analysis: 
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dz 



sin 
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K.j 
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(7.4.28) 



where the upper integration boundary z c = \nco' c Ao. Note that, as we indicated 
earlier, the use of fij /a / in the definition of z\ iW i is a conservative stance. Indeed, 
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strictly speaking, we should write Z/ /Ct / ex (u/ — Ui /Ct /). For a fixed z c , a smaller value 
of Z/ ^,/, and hence a larger value of means that a larger part of the central peak 
of the integrand will be integrated over. Since we are replacing Uj^r by the upper 
bound Ui jW i, we are overestimating the resulting radiation, i.e., underestimating the 
modification with respect to the standard Hawking radiation. 

The expression ( |7.4.28 > brings out the two crucial factors mentioned earlier, and 
a third, corollary one. 



First, it shows the dependence of the total radiation on the critical frequency cv' c 
(through the integration boundary z c induced by it), as discussed just before 
obtaining formula (7.4.27). 



• Second, it shows the importance of the frequency-dependent k w i (to be com- 
pared with the fixed k of the standard case). Given a concrete profile v(t,x), 
the frequency-dependence of k w i can be derived explicitly, as we will discuss 
next. 



• Finally, as a corollary of the first point, it shows that, as U\ (and hence Z/^y) 
increases, a smaller part of the central peak of the integrand will be integrated 
over, so the radiation will die off as u\ advances. 



It is also a simple exercise to consider what this result implies for the case of a con- 
stant surface gravity k w i = kq. From ( |7.4.28| >, one obtains immediately that the ex- 
ponential Planckian factor leaves the integral, whereas the remaining terms lead to 
a pref actor smaller than one. 



7.4.6 Surface gravity 



A careful analysis has shown that, so far, most formulas for the standard case could 
be adapted to superluminal dispersion relations by replacing the relevant magni- 
tudes with their frequency-dependent counterparts. For example, f5 wco i was ob- 
tained by replacing U and u by U w i and u w , and in particular Uh, u\ and k by U.n,u'> 
Ui r!j0 i and k w /, respectively. Given a concrete profile, we can explicitly deduce the 
relation between k w i and to' , with the Lorentz violation scale ki as a parameter, as 
follows. The horizon for a particular initial frequency a/ is formed when 
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i.e., when v p h = c^ + v = 0, where we have used the phase velocity for concreteness. 
But again, qualitatively, our results would be similar if we take the group velocity 
v g instead of v^. 

Taking into account that fc# = k', the dispersion relation in the asymptotic past 
can be written as 



co 



12 



1 2 7,2 



\V 2 (X H ,U>>)\ 



1 



(7.4.30) 



Given a concrete profile, x H/C0 > can then be obtained and k w i calculated explicitly. 



For the profiles of the first type, see Fig. 7.1 the result is trivial. Indeed, \v(x)\ 
increases linearly between the horizon corresponding to co' = and the one for co' c , 
so we obtain a constant k w i = kq for all co' < co' c . 



For a Schwarzschild profile as in Fig. 7.2 on the other hand, 



v(x) 



c 2 /2k 
x + c 2 /2k ' 



(7.4.31) 



so we obtain 



dv 



dx 
1 



*0 
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2V2 



1 + 4/1+4 



CO 
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3/2 



3/2 



c 2 k\ 



(7.4.32) 



Note that k w i remains nearly constant until frequencies of the order of magnitude 
of the Lorentz violation scale are reached, and then starts to increase rapidly, see 
Fig. |7.5[ As we will see graphically in the next section, this can have an important 
qualitative influence on the radiation spectrum. 

We now have all the tools necessary to compute and plot ( |7.4.28 > for different 



parameters of the profiles described in section 7.2 



7.5 Graphical results and discussion 



We have integrated eq. ( |7.4.28 > numerically using the Gauss-Chebyshev quadrature 
method. The results are plotted in Figs. 7^6]- 7.8 and perfectly illustrate the three 



important factors that we deduced theoretically in section 7.4.5 Note that in all the 
figures we have plotted E = co 3 x N against co to make visual comparison with the 
usual thermal spectra in 3+1 dimensions easier. 
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Figure 7.5: Surface gravity k w i with respect to the (logarithm of the) frequency for 



a Schwarzschild-type black hole as in Fig. 7.2 with the Lorentz violation scale ki 
10 39 t _1 and the critical frequency co' c 
time unit.) 



13 x 10 39 t 1 (where t represents an arbitrary 



Fig. 7.6 shows the influence of the critical frequency co' c for a profile with constant 



K w i, as in Fig. 7.1 and U\ = «j /Ct / (i.e., immediately after the horizons have formed 
for all oJ < u)' c ). The constant k w i guarantees that the form of the Planckian spec- 
trum is approximately preserved. In particular, the peak frequency is not shifted. 
However, the intensity of the radiation decreases with decreasing critical frequency. 
Actually, only for extremely high critical frequencies is the original Hawking spec- 
trum recovered. For a critical frequency still well above the Lorentz violation scale, 
the decrease can be significant. For example, for co' c = 10 61 cki, the peak intensity 
decreases by nearly 30%, while for co' c = cki, the decrease is approximately 40%. 
At the other end, note that extremely low critical frequencies still leave a significant 
amount of radiation. For example, for co' c = 10~ 139 cki, still 20% of the original peak 
radiation is obtained. So Hawking radiation receives significant contributions from 
an extremely wide range of frequencies. 



Fig. 7.7 shows the radiation spectrum for a Schwarzschild-like profile and hence 
increasing k w /, as in Fig. 7.2 in particular for critical frequencies co' c close to the 
Lorentz violation scale. On top of the general decrease of the standard Planckian 
part of the spectrum (approximately 40%, as in the previous case) due to the fi- 
nite integration boundary induced by the critical frequency, the fact that the surface 
gravity is now frequency-dependent leads to an important qualitative change of the 
spectrum. The high-frequency tail of the spectrum is totally transformed. Actually, 
if the critical frequency is sufficiently higher than the Lorentz violation scale, the 
dominant source of radiation lies in the high-frequency region. Note that this effect 
is truly a consequence of the modification of the physics for frequencies above the 
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CO 



Figure 7.6: Influence of the critical frequency on the radiation spectrum for a black 



hole with velocity profile such that the surface gravity is constant, as in Fig. 7.1 and 
different values of the critical frequency. We have chosen c = 1, u\ = Uj = 0, Uh = 
A = 1 t, kq = 12 t _1 and the Lorentz violation scale ki = 10 39 t _1 (where t represents 
an arbitrary time unit), which amounts to considering a solar-mass black hole. From 
top to bottom we have plotted co' c = 10 3000 , 10 100 , 10 39 , 1, 1(T 39 , 1(T 100 , 1(T 300 . Note 
that the standard Hawking spectrum coincides perfectly with the upper curve, which 
effectively corresponds to the absence of a critical frequency. 



Lorentz violation scale. This can be appreciated by noticing that, at co' c = O.lcki, 
the whole tail-modifying effect has disappeared and the usual Planckian form of 
the thermal radiation spectrum is approximately recovered, although with a lower 
intensity, as the quantitative decrease described above still occurs. 



Finally, Fig. 7.8 shows the influence of the measuring time u\ (measured with re- 



spect to Ui /(V i) for a profile of the second type (increasing k w /). It is clearly seen that 
the radiation originating from the collapse process dies off with time, and actually 
dies off rather fast. For a solar-mass black hole and an co' c of the order of the Planck 
scale, the radiation would last less than a second. Note that this effect is a corollary 
of the existence of a critical frequency co' c , since in its absence, the integration bound- 



aries in eq. ( 7.4.28 1 would be infinite, and so the integral would be insensitive to a 
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Figure 7.7: Influence of a frequency-dependent surface gravity k w < on the radiation 
spectrum for a black hole with Schwarzschild-type velocity profile (surface gravity 
increases with frequency), as in Fig. 7.2 and different values of the critical frequency 
around the Lorentz violation scale kf. co' c / ck-^ = 13, 10, 7, 4, 0.1 (from top to bottom). 
The standard Hawking spectrum is depicted in dashed-dotted line for comparison. 
Numerical values of Kq, k^ etc as in Fig. |7.6| 



change u — > u + Au. 

This decrease with time is not wholly unexpected [58: 65 J. Indeed, the redshift 
in the frequency suffered by a wavepacket moving away from the black hole is ex- 
ponential, in analogy with the relation between the past and future null coordinates 



(see eqs. (7.3.161 and (7.4.25)). In terms of original frequencies co 1 and final frequen- 
cies co, this redshift can be written as 

-c 3 w/4M 



co oc co 'e 



(7.5.1) 



in units such that G = 1, and where u is the time of arrival at the asymptotic future 
infinity. Therefore, when tracing modes with a fixed final frequency co back in time, 
the corresponding original frequency co' is exponentially blueshifted with increas- 
ing u. In the case of superluminal dispersion, this implies that, as time advances, 
one inevitably reaches a moment when original frequencies above the critical fre- 
quency would be needed in order to sustain the radiation. But since these ultra- 
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Figure 7.8: Influence of the measuring time U\ on the radiation spectrum, for a 

and Lo' c /cki 



Schwarzschild-type velocity profile, see Fig. 



7.2 



13, compare with 
0,22,35,50,500 t with 



Fig. 7.7 Different values of u\ (from top to bottom): u\ 
t denoting an arbitrary time unit. The standard Hawking spectrum is depicted in 



dashed-dotted line for comparison. Numerical values of Kq, ki etc as in Fig. 7.6 



high-frequencies do not see a horizon, they do not suffer such an exponential shift, 
hence they simply do not contribute to the Planckian radiation. 



The combined effect of Fig. 7.7 and Fig. 7.8 leads to the following qualitative pic- 
ture for the further collapse towards a singularity once the initial or classical horizon 
has formed. As the interior gradually uncovers a larger portion of the Schwarzschild 
geometry, two competing processes will take place. On the one hand, the spectrum 
acquires ever larger contributions associated with higher and higher temperatures. 
On the other hand, the overall magnitude of the spectrum is damped with time. 
The question of which process dominates could depend on the fine details of the 
dynamics of the collapse, and might be further complicated by backreaction effects, 
which we have not considered in our analysis. 

It is anyhow clear, particularly for the case of a Schwarzschild collapse, that the 
relation between the critical frequency and the Lorentz violation scale can play a 
crucial role. Since the critical frequency is associated with the level at which the col- 
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lapse process saturates, it is reasonable to expect that, for the case of a gravitational 
collapse, this critical frequency will ultimately be related to the Planck scale E P i. 
Indications this expected relationship can be found both in dimensional estimates 
for quantum gravity based on quantum mechanics, special relativity and general 
relativity 115511 , and in toy models, for example within the context of loop quantum 
gravity [|27|. Since current observations seem to indicate that Lorentz violations at 
the Planck scale are strongly suppressed (if present at all), and hence Ei E P j, 
it seems that at least the Planckian form of the black hole radiation is rather robust, 
even if there turns out to exist superluminal dispersion at high frequencies. The full 
intensity of a truly Planckian spectrum, however, is only recovered if the critical fre- 
quency co' c —> +00, in other words, if there is a singularity at the end of the collapse, 
and if moreover the surface gravity is constant for all frequencies (which, we again 
insist, in a superluminal scenario is not the case for a Schwarzschild collapse). 

However, this robustness of the Planckian form of the spectrum seems to be 
more 'accidental' than intrinsic (unless a fundamental explanation can show why 
Ei ^> Epi should necessarily be the case). To illustrate this, remember that for 
BECs, we argued in chapter [6] that the analogue of the Planck scale E P i is related to 
the interatomic distance, which provides an intrinsic short-distance cutoff, whereas 
the Lorentz violation scale Ei is determined by the healing length £, such that 
Ei <C E P j. It is easy to see, then, that the mere creation of a black hole horizon 
in a BEC already implies that the critical frequency must lie above the Lorentz vio- 
lation scale Q Hence we expect ultraviolet contributions to spontaneously show up 
in any Schwarzschild-type collapse experiment in a BEC in which the collapse suffi- 
ciently crosses the zero-frequency horizon. In practical terms, in the model that we 
are considering, this involves an increasingly accelerating background flow veloc- 
ity, see Fig. 7.2 up to at least a few times the speed of sound. On the other hand, if 
the background flow velocity has a linear slope, as in Fig. 7.1 then such ultraviolet 
contributions should be absent. 

In any case, care should be taken when proclaiming the robustness of Hawking 
radiation in collapsing configurations. Indeed, our results show that, while it might 
be possible to recover the standard Hawking spectrum in some particular limit, it is 
certainly not the general rule in the case of superluminal dispersion. 



4 A straightforward calculation taking for example v = 2c and using the phase velocity | |7.4.4) 

4v / 3cfcj j . From Fig. 7.7 we see that this might already 



and the dispersion relation 1 7.4.2 1 gives co' c 
be sufficient to detect the ultraviolet contributions 
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7.6 Summary and discussion 



We have discussed the Hawking radiation for a collapsing configuration with super- 
luminal dispersion relations. Our approach followed the lines of Hawking's original 
calculation through the relation between the ray trajectories in the asymptotic past 
and those in the asymptotic future. We did therefore not need to impose any con- 
dition or assumption at the horizon. Rather, our main assumptions were the usual 
ones with regard to a gravitational collapse, namely a Minkowski geometry in the 
asymptotic past, and asymptotic flatness also in the future. 

An additional implicit assumption in our calculations is the following. We have 
assumed that the number of degrees of freedom of the quantum theory associated 
with the collapsing configurations is equal to that of the field theory in a station- 
ary Minkowskian spacetime, which is the starting point of the dynamical evolu- 
tion of our profiles. In this context, note that for a constant subsonic flow, given a 
real frequency co, the modified scalar field equation has two independent normalis- 
able solutions (with the standard caveats associated with the normalisation of plane 
waves), whereas the other two solutions become divergent in some region and hence 
must be discarded. However, in a constant supersonic flow, depending on the fre- 
quency range, one can find either two or four independent normalisable solutions 
(see Fig. 7.9 1, with a transition point in which there are three solutions (of which one 
is degenerate). Upon quantisation, these theories have a different number of degrees 
of freedom. Since a stationary black hole configuration consists of a subsonic and a 
supersonic region, the number of normalisable solutions for a given frequency can 
then either be two or three, depending on the frequency under consideration. But 
there always exist frequencies with three associated normalisable solutions. So the 
field theory over a stationary black hole configuration has more degrees of freedom 
than that for a Minkowskian (i.e., purely subsonic) spacetime. Our assumption can 
then be rephrased as follows. We assume that the introduction of dynamics does not 
change the number of degrees of freedom of the initial theory, and so we consider 
only how the two initial Minkowskian degrees of freedom are distorted by the dy- 
namics of collapse, and how this influences the standard mechanism for Hawking 
radiation in such a collapse process. 

Our results can be summarised as follows. 

Modifications of the dispersion relation cause the horizon, and various associ- 
ated quantities such as the surface gravity k, to become frequency-dependent. In 
particular, a critical frequency co' c naturally appears such that frequencies higher 
than co' c do not experience a horizon. More generally, it also means that the stan- 
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Figure 7.9: Superluminal dispersion relation for various values of v (left) and of co 
(right), scaled at c = 1. The intersection points between the line co — vk (in blue) and 
the curves ±ck^l + k 2 lk\ (in black) mark the real (normalisable) mode solutions 
for a given co. For a given frequency coo (left), there are always 2 normalisable solu- 
tions in a subsonic regime (\v\ < c), but there can be up to 4 in a supersonic regime 
(\v\ > c). For a given supersonic flow (right), depending on the value of co, there can 
be 2 or 4 normalisable solutions. 



dard geometric concepts traditionally used to study black holes must be handled 
with care. Nevertheless, through a detailed analysis, we have seen that the equa- 
tions related to the late-time radiation can be adapted quite straightforwardly from 
standard (relativistic) to superluminal dispersion relations. 

We analytically derived an approximate equation for the particle production 
at late times ( 7.4.28| > with superluminal modifications of the dispersion relations. 
This equation clearly shows that important modifications in the late-time radiation 
should be expected, first, due to the existence of the critical frequency co' c and the fi- 
nite upper boundary it induces in the integration, and second, due to the frequency- 
dependence of k w /. We integrated eq. ( |7.4.28| > numerically, and plotted the resulting 
spectrum, thereby confirming these expectations. 

We have seen that the standard Hawking spectrum is recovered only in a very 
particular case: in the limit when the critical frequency goes to infinity (i.e., when the 
profile for the velocity \v\ goes to a singularity) and moreover the surface gravity k u / 
is constant (linear velocity profiles of the type of Fig. 7.1 ). For lower critical frequen- 
cies, as long as k remains constant, the Planckian form of the Hawking spectrum 
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is roughly maintained and the peak remains at the same frequency, but the peak 
intensity decreases rapidly with decreasing co' c . A non-negligible radiation persists, 
however, even for extremely low critical frequencies. This clearly establishes the 
first point, namely the importance of the critical frequency to' c , which illustrates the 
statement made in the introduction that with superluminal dispersion relations, the 
interior of the black hole is probed and significantly affects the Hawking process. 

For a Schwarzschild collapse, the velocity profile generally leads to a frequency- 
dependent surface gravity k w i, increasing with the frequency a/. This means that 
the standard Hawking result cannot be recovered for a Schwarzschild black hole 
with superluminal dispersion relations, even when the velocity profile has a singu- 
larity. Actually, when the velocity profile has a limiting value, the same quantitative 
decrease of the Hawking part of the spectrum as before shows up. Moreover, if 
the critical frequency co' c is above the Lorentz violation scale, a drastic qualitative 
change of the radiation spectrum takes place, and for sufficiently high co' c (roughly a 
few times the Lorentz violation scale) the high-frequency part of the spectrum even 
becomes dominant. This shows the importance of the surface gravity, which again 
illustrates the role played by the interior of the black hole. We have indicated that 
the critical frequency is likely to be related to the Planck scale, or its analogue (asso- 
ciated with the interatomic distance) in the case of BECs. Observational limits on Ei 
from cosmology (namely, Ei ^> E P j, see chapter [6]) then suggest that astrophysical 
black holes would not suffer such strong ultraviolet (UV) contributions. However, 
this does not mean that the Hawking radiation is intrinsically robust to them. Since 
for BECs, the relation between Ei and Ep/ is precisely the opposite (E P j ^> El), the 
phononic radiation from a collapsing experiment in a BEC could in principle show 
such important UV contributions. 

Finally, we have also seen that, as a corollary of the existence of a critical fre- 
quency and of the finite upper boundary induced by it in the thermal response func- 
tion, the radiation spectrum of the collapse dies off as time advances. This effect can 
be understood by the fact that frequencies are exponentially redshifted when mov- 
ing away from the black hole. Therefore, when tracing back the outgoing radiation 
at a time u to the original frequencies that created it, these original frequencies in- 
crease exponentially with u. Inevitably, as soon as this tracing-back process reaches 
the critical frequency, the radiation is effectively switched off. This cut-off effect has 
long been understood to be one of the crucial aspects of the transplanckian problem 
related to Hawking radiation. Here, however, the cut-off is not imposed by hand to 
solve the transplanckian problem, but appears naturally because of the superlumi- 
nal dispersion. 
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A few observations might be useful to connect our work with existing results on 
Hawking radiation and its sensitivity to modified dispersion relations. The most 
general observation is that the 'robustness' of Hawking radiation, which is often 
considered to be a well-established result, actually depends crucially on a series of 
assumptions, usually imposed at or near the horizon. Such assumptions might be 
reasonable in the case of subluminal dispersion relations. But for the case of su- 
perluminal dispersion relations in a collapse scenario, as we have shown explicitly, 
the assumptions needed to reproduce the standard Hawking result depend on the 
physics inside the (zero-frequency) horizon, and moreover in a way which is not 
compatible with the Schwarzschild geometry. 

In particular, in |4T|, it was shown for a stationary scenario that Hawking ra- 
diation is robust with respect to superluminal modifications of the dispersion rela- 
tions, provided that positive free fall frequency modes were in their ground state 
just before crossing the horizon. However, it was also admitted that it is not clear 
whether this is the physically correct quantum state condition. In 111161 , three explicit 
assumptions were given for the previous condition to hold: freely falling frame, 
ground state and adiabatic evolution. In slightly more detail: 



1. If the Lorentz symmetry is violated in such a way that there is a preferred 
reference frame, then this should be the freely falling frame (and not the rest 
frame of the black hole, for example). 

2. High-frequency excitations must start off in their ground state (with respect to 
the freely falling frame) before leaving the region infinitesimally close to the 
horizon on their way out. 

3. The dynamics of all internal degrees of freedom should evolve much faster 
than any external time dependence, such that no high-frequency modes are 
excited by these external variations. 

In the same paper |116J, it is similarly pointed out that these assumptions might fail 
for a superluminal dispersion relation, since (as we also mentioned in the introduc- 
tion) for high-frequency modes this implies that one needs to make an assumption 
about the physics at the singularity. 

Rather than imposing any conditions on the state near the horizon, and hence 
ultimately at the singularity, we have fixed the initial geometry to be Minkowskian 
and the quantum field to be in the natural associated vacuum state, and evolved 
this into a black hole configuration. Since our results are seemingly in contradiction 
with the ones obtained in [41] and 111161 , it is worth examining in which sense the 
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conditions stated by those authors are violated in our approach. If one considers 
a collapse scenario, for example of a BEC, in a laboratory setting, then it is quite 
natural to assume that the 'freely falling frame' condition will be violated. Indeed, 
the superluminal modification is in this case associated with the existence of a privi- 
leged external reference frame: the black hole rest frame or lab frame (as we noted in 
sec. 7.4.1 1, and not the freely falling frame, as assumed in II116I1 . Note that this viola- 



tion of the free-fall frame condition automatically implies a violation of the ground 
state condition in the sense in which this condition is formulated in ||116| , namely as 
the ground state "with respect to the freely falling frame". 

Our results then suggest the following picture. The low-energy modes experi- 
ence the classical geometry, are therefore dragged along in the free-falling frame 
and so at the horizon they occupy the vacuum state associated with this free-falling 
frame, namely the Unruh vacuum. Hence they contribute to the black hole radi- 
ation in the traditional Hawking way. However, the ultra-high-frequency modes 
(above the critical frequency) do not see the horizon. Hence they do not couple to 
the classical geometry of the collapse, but rather remain connected to the external or 
laboratory frame, and therefore pass through the black hole (nearly) undisturbed. 
The ground state of these high-frequency modes is then not the vacuum associated 
with the freely falling frame, but the Boulware vacuum associated with the orig- 
inal Minkowski geometry, or in other words, with the stationary reference frame 
of the lab. So these frequencies above the critical frequency do not contribute to 
the thermal output spectrum. The overall radiation is a convolution of the contri- 
butions from all the different initial frequencies, where the surface gravity for each 
frequency can be interpreted as leading to an effective weighting factor. Depending 
on the internal physics of the black hole, this leads to either a spectrum of the tradi- 
tional Hawking form but with a reduced intensity (in the case of a constant surface 
gravity), or a modified spectrum where the high-frequency contributions dominate 
(as for a Schwarzschild profile, provided that the critical frequency lies sufficiently 
above the Lorentz violation scale). 

In any case, what our analysis shows with respect to the situation treated in [ j4l| 
and 111161 is the following. The robustness of Hawking radiation discussed by those 
authors is the robustness with respect to a superluminal modification under the 
precise conditions occurring in the standard case, in which all modes occupy the 
ground state associated with the free-falling frame just before crossing the horizon. 
This assumption implies a stationary configuration in which all modes, regardless 
of their frequency, occupy the Unruh vacuum. Our analysis precisely shows that 
the superluminality can lead to a spontaneous breaking of these conditions. The 
natural evolution of a collapsing configuration is such that the frequencies above the 
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critical frequency co' c are not in the Unruh vacuum associated with the freely falling 
frame, but remain in the Boulware vacuum associated with the initial asymptotic 
Minkowski condition. It should then come as no surprise that this can indeed have 
an important impact on the resulting radiation spectrum. 

An additional remark, already mentioned at the beginning of this section, is that 
we have only considered the evolution of the two Minkowskian degrees of freedom, 
whereas |[4"T1|116I1 deal with stationary configurations, and hence with a higher num- 
ber of degrees of freedom. It would certainly be interesting to see if it is possible to 
connect both situations in a consistent way. This would require the description of 
a possible transition from two to three degrees of freedom, which is not within the 
direct scope of the formalism that we have used, since this is essentially based on a 
ray-tracing method, see section 7.4.4 Our approach would obviously not be valid in 
the case of a stationary black hole configuration, since there the rays would simply 
break down at the horizon. Here, however, we have traced the rays that escape just 
before the formation of the horizon, and therefore reach the asymptotic infinity with- 
out ever having encountered a turning point that would signal such a breakdown. 
We therefore emphasise that our calculation considers exactly the same mechanism 
that produces Hawking radiation in the standard (Lorentz-invariant) collapse case. 

Our approach is therefore conservative in the sense that we have studied how 
the standard mechanism is deformed in the case of superluminal dispersion. The 
results that we have obtained are analytical. Indeed, although we have illustrated 
our main points graphically through numerical integration, the qualitative content 
of our results was obtained by inspection of formula (7.4.28) (or (7.4.27)). The fact 
that we have limited ourselves to study the Minkowskian degrees of freedom, or 
in other words: that we have used a ray-tracing method, means that our method is 
also approximative, since we have not taken into account the possible transition to a 
higher number of degrees of freedom. 

Our result should therefore be understood in a dual sense. 

First, we have established that the standard mechanism of Hawking radiation 
in a collapse scenario is generally not robust under superluminal dispersion. More- 
over, the equivalence between collapsing and stationary black holes (at least, in the 
straightforward extension from the standard case) is no longer valid in general, since 
the Unruh state is not necessarily the natural end-state of the gravitational collapse 
process. 

Second, it would be extremely interesting to see to what extent it is possible to 
go beyond the approximation that we have used, and see to what extent all degrees 
of freedom can be treated in a consistent way. In this sense, our work will hopefully 
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stimulate further investigation of these issues. 

A final observation concerns the recent article 11351 , in which the authors nu- 
merically simulate the creation of an acoustic horizon in a BEC, and analyse the 
production of correlated pairs of phonons through the so-called truncated Wigner 
method. Our findings for the case of a constant surface gravity seem to be in qual- 
itative disagreement with the discussion presented in 11351 , since the authors of this 
paper assert having observed a stationary Hawking flux. The source of this ap- 
parent discrepancy might reside in the fact that the correlation function which they 
study is normalised, and therefore probes the form of the spectrum (which we also 
find stationary with time), but does not provide details about the net amount of par- 
ticle production. Alternatively, it could be that the numerical method used in 11351 
automatically includes the creation of an additional, 'non-Minkowskian' degree of 
freedom. In that case, it would be extremely interesting to see whether we can un- 
derstand the exact physical mechanism that compensates for the decrease in the 
radiation coming only from the Minkowskian ones, such as in our analysis. A last 
remark is the following. If an experiment is realised in which the analogue Hawk- 
ing radiation in a BEC can indeed be detected, for example along the lines suggested 
in |[35j, this would certainly be a major advance (as well as a tremendous source of 
publicity). Additionally, it might help us to gain further insight into the (so far 
theoretical) discussion about the importance of the additional modes, possible ul- 
traviolet contributions and so forth. 



Chapter 8 

Emergent gravity from condensed 
matter models 



8.1 Introduction 



The various case studies that we presented so far were either derived from or at least 
inspired by the gravitational analogy in condensed matter systems, and in BECs in 
particular. To interpret these case studies as quantum gravity phenomenology as 
we did, assumes that Lorentz invariance will turn out to be broken by quantum 
gravity in a way which can, in a first approximation, be modelled by a superlumi- 
nal modification of the relativistic dispersion relations at high energies. In this final 
chapter, we will partly address the question of how seriously the condensed mat- 
ter analogy can or should be taken, and indicate some conceptual issues related to 
the construction of a model for emergent gravity based on the condensed matter 
analogy. 



A naive but nevertheless crucial observation in this respect is the following |131| . 
The fundamental 'quantum gravity' theory is generally assumed to have the Planck 
level as its characteristic scale. Expressed as a temperature, this Planck level lies at 
approximately 10 32 K. On the other hand, almost all of the observable universe has 
temperatures that barely exceed the cosmic background radiation temperature of a 
few Kelvins. Even the interior of a star such as the sun is more than 20 orders of 
magnitude colder than the Planck temperature, while the highest energies that are 
planned to be produced at the Large Hadron Collider are still roughly 15 orders of 
magnitude lower than the Planck scale. So the degrees of freedom of quantum grav- 
ity, independently of their fundamental structure, are probably effectively frozen out 
in most of our universe, just like in a condensed matter system in a low-temperature 
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laboratory. The physics that we observe might then well be due to collective exci- 
tations that result from the — comparatively tiny — thermal (or other) perturbations 
of this vacuum. If this observation is indeed relevant for quantum gravity, then one 
crucial consequence would be the following. In a condensed matter system, there is 
no direct quantisation path from the perturbations or 'effective' degrees of freedom 
(the phonons and quasi-particles) to the 'fundamental' degrees of freedom: the con- 
densed matter atoms. Likewise, a basic tenet of emergent gravity is that there might 
not be any direct path such as a quantisation procedure from the effective degrees 
of freedom of classical gravity to the fundamental ones of quantum gravity. 

Moreover, it is important to realise that many of the low-temperature effects that 
occur in condensed matter systems are not specific to one or another type of atoms 
or molecules, but are generic to a wide range of similar systems, e.g., bosonic sys- 
tems in the weakly interacting regime, or fermionic systems within the same low- 
temperature topology class in momentum space 1112.811 (of which the Fermi-point 
topology seems to be the most relevant for our universe |I124| , and is highly ro- 
bust with respect to fluctuations in the microscopic structure), in combination with 
some general features such as the relation between the characteristic scales (the in- 
teratomic distance and the scattering length, for example). From this point of view, 
knowledge of the relations between the characteristic scales such as the Planck scale 
and the Lorentz symmetry violation scale might actually be more relevant to calcu- 
late the high-energy corrections to classical gravity than an exact description of the 
fundamental, 'quantum-gravitational' degrees of freedom. In this sense, classical 
gravity should perhaps be interpreted much as thermodynamics 119811 : a set of rules 
for the collective behaviour at low temperatures of some degrees of freedom whose 
precise microscopic constitution can vary quite widely within a certain universality 
class. It is in this sense that we use the term 'emergent gravity' here. 

In the following, we will concentrate on one particular issue of this general idea 
of emergent gravity from condensed matter models, namely diffeomorphism in- 
variance. Indeed, one problem that immediately comes to mind from a relativistic 
point of view with respect to the condensed matter models is that they are 'back- 
ground dependent' in a trivial way: the metric description that emerges for the mo- 
tion of the phonons depends on the collective behaviour of the condensed atoms. 
In other words, the background is precisely determined by the condensed phase of 
the condensate. However, in the context of loop quantum gravity, for example, it 
has been argued that diffeomorphism invariance, understood as a requirement for 
the non-existence of any prior geometry, is an essential feature of general relativity, 
and should therefore also be an essential feature of quantum gravity l|104t 11091 . It 
seems that the background dependence of condensed matter models obstructs the 
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possibility of extending them to a diffeomorphism invariant model for gravity. We 
will see how this issue can at least partly be resolved, and discuss its relation to the 
problem of recovering the Einstein equations in models of emergent gravity 

First, though, although quite remote from the nucleus of this thesis, it is probably 
worth briefly mentioning an important reason why condensed matter models have 
attracted interest in the domain of quantum gravity, namely the accelerated expan- 
sion of the universe, which seems to require the existence of some form of repulsive 
'dark energy'. 



8.2 Quantum gravity and dark energy 

Quantum gravity is in the first place the endeavour to reconcile quantum mechanics 
and the general theory of relativity. There are many arguments why such a unifica- 
tion is needed, see e.g. [32j. The bulk of these go along the line of consistency argu- 
ments, for example that intrinsically classical fields are incompatible with quantum 
mechanics and could lead to violations of Heisenberg's uncertainty principle. How- 
ever, it is far from clear whether any such argument could have real empirical conse- 
quences, even in principle [33]. Another common argument is that general relativity 
predicts singularities, not only in mathematically exotic cases but also in situations 
of physical relevance such as black holes and initial 'big bang' singularities in cos- 
mological models. It is then generally assumed that a quantum theory of gravity 
will relieve this problem and avoid the occurrence of singularities. Nevertheless, 
again, this argument is far from water-tight: first of all, it is not obvious that a quan- 
tum theory of gravity is really needed to cure the singularities of general relativity. 
For instance, purely algebraic extensions of general relativity could solve the singu- 
larity problem without any need for a quantisation procedure ||79B . Furthermore, it 
could well be that gravity is intrinsically a non-quantum phenomenon Il29l , which 
simply does not apply at the energy scales at which the singularities are predicted 
to occur. The approach discussed in this thesis is close in spirit to other suggestions 
that gravity should not be quantised, for example thermodynamic interpretations 
of the Einstein equations H67H100I1 or Sakharov's induced gravity 111061 (about which 
we will have more to say in section 8.4[ >. 



Recently, it has become clear that, apart from the quest for a unification of general 
relativity and quantum mechanics, there is at least one empirical issue of importance 
related to (quantum) gravity: the accelerated expansion of the universe. This seems 
to require the existence of some sort of repulsive energy, generically baptised 'dark 
energy' [97] . It has been argued that the problem of dark energy indicates the need 
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for a reassessment of the traditional paradigm of quantum gravity, for example in 
the following way [99] . First, although there is a myriad of alternative approaches 
to dark energy, the cosmological constant is in many ways the most sensible can- 
didate, if only because none of the other approaches explains why the cosmologi- 
cal constant should be (approximately) zero, i.e., why the zero-point energy of the 
quantum matter fields is not huge, as a naive calculation from quantum field theory 
seems to indicate 111321 . Second, the equations of motion for matter fields are invari- 
ant under a shift of the matter Lagrangian C m by a constant. But the gravitational 
sector is obviously not invariant under such a shift. In [99], it was concluded that 
we need to rethink gravity in a way which can be described in terms of an effective 
action which is explicitly invariant under a shift of the energy-momentum tensor 
T-uv ~^ T„v + Ag„ v . The condensed matter perspective offers a related point of view, 
which shows that the original intuition that dark energy or the cosmological con- 
stant is simply the zero-point energy of all the quantum fields, i.e., the energy of the 
quantum vacuum, might prove to be right after all. The essence of the argument is 
the following H25HI2gHI27 l. 

The value of the quantum vacuum energy in a condensed matter system in equi- 
librium is regulated by macroscopic thermodynamic principles. The vacuum en- 
ergy density e vac = E vac /V of a quantum many-body system is obtained from the 
expectation value of the corresponding Hamiltonian which describes the effective 
low-temperature quantum field theory in second quantisation: Hqft = H — jiN. 



Here, H is a Hamiltonian of the form (2.2.2 1, N is the atom number operator and the 



chemical potential ]i can be understood as a Lagrange multiplier which accounts 
for the conservation of the number of atoms Q So e vac = (H-piN) vac /V. The 
equation of state relating the energy density and the pressure of the vacuum of a 
quantum many-body system can then be obtained through the Gibbs-Duhem equa- 
tion E — jiN — TS = —pV, which for a pure vacuum state (T = 0) gives simply 
e V ac = — Pvac- This equation of state corresponds to a cosmological constant. 

It seems reasonable to require that the universe should be describable as an iso- 
lated system, i.e., without recurring to any external quantities such as an external 
pressure p e xt- This suggests the analogy with a liquid-like system, since liquids can 
be in a self-sustained equilibrium without external pressure. If p e xt = 0, then for 
a liquid-like pure vacuum state in equilibrium, also p vac = 0. The natural value 



1 If the relativistic quantum field theories of the standard model of particles can indeed be inter- 
preted as effective theories in a similar way, then the equivalent to this condition would of course 
not be a conservation of particle number, since the particles of the standard model correspond to the 
quasiparticles of the condensed matter analogy, but would rather be related to the number of micro- 
scopic degrees of freedom of spacetime. 
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for e vac at T = in such a system in equilibrium is then e vac = 0. In other words, 
according to this argument, a pure vacuum state in equilibrium in the absence of 
external pressure is non-gravitating. At T 7^ 0, thermal fluctuations lead to quasi- 
particle excitations, and hence to a matter pressure pm- In the absence of external 
pressure, this matter pressure is compensated by a non-zero vacuum pressure in 
order to restore the equilibrium: p vac + p^ = 0. The vacuum energy therefore natu- 
rally evolves towards the value e vac = Pm in equilibrium |[T9][, and thereby becomes 
gravitating. The system automatically adjusts itself to obey the macroscopic thermo- 
dynamic rules, and there is no need to know the precise microscopic constitution of 
the system to calculate the macroscopic quantities. 

From such a condensed-matter point of view, the cosmological constant mystery 
becomes a lot less unsurmountable: from having to explain why the cosmological 
constant is more than a hundred orders of magnitude smaller than the value based 
on a naive zero-point energy calculation B132L it is reduced to having to explain 
why it is slightly bigger than the equilibrium value which would exactly cancel the 
total (baryonic + dark) matter contribution: Q A « 0.7 versus Qm ~ 0.3 11971 in terms 
of density parameters O, = p\l p c , with p c the critical density for a flat universe 
(which seems to correspond very closely to the actual density of our universe). So, 
the condensed matter approach offers a qualitative framework to understand dark 
energy, at least in terms of the coincidence problem (why Oa ~ Om) an d the cos- 
mological constant problem (why the zero-point energy of the vacuum is not huge). 
Whether it also allows to extract a concrete quantitative prescription for the evolution 
of dark energy is an open problem, to which a possible first step has been described 
in G1IZ3. 

We now turn to the main subject of this chapter: diffeomorphism invariance in 
emergent-gravity models. 

8.3 Diffeomorphism invariance: kinematical versus 
dynamical aspects 

One of the important lessons from analogue gravity is to remind us of the fact that 
general relativity can be split up into kinematical and dynamical aspects. The kine- 
matical aspects of general relativity are encoded in the curved Lorentzian space- 
time, i.e., a manifold with a metric signature (—,+,+,+). The dynamical aspect of 
general relativity is essentially the interplay between this curved geometry and the 
matter distribution through the Einstein equations ("Geometry tells matter how to 
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move, and matter tells geometry how to curve" [90J). The crucial lesson from ana- 
logue gravity in this respect is that the emergence of a curved Lorentzian spacetime 
occurs naturally in a wide variety of systems, in particular in several types of con- 
densed matter systems, i.e., when a many-body system is sufficiently cooled. So the 
curved Lorentzian aspect of general relativity is in a certain sense 'natural'. How- 
ever, it seems that what really makes our universe very different from the effective 
spacetimes that are obtained in analogue models are the Einstein equations. No real 
laboratory system (condensed matter or other) reproduces the Einstein equations. 
One of the main objectives of the emergent-gravity programme (or 'paradigm', to 
use the Kuhnian terminology [77]), is precisely to elucidate the origin and the natu- 
ralness of the Einstein equations. 

When discussing background independence and diffeomorphism invariance, a 
similar distinction between kinematical and dynamical aspects can be made. Tech- 
nically speaking, diffeomorphism invariance is the requirement that the laws of 
physics should be formally invariant under transformations defined by a diffeo- 
morphism, i.e., an invertible bijective and smooth function / with smooth inverse 
f~ l defined on two differentiable manifolds. Diffeomorphism invariance (or general 
invariance, which should not be confused with general covariance)^is understood by 
most contemporary relativists as an essential feature of general relativity. This is 
why quantum gravity programmes such as loop quantum gravity and related spin 
foam approaches have developed a canonical quantisation scheme for general rela- 
tivity which is diffeomorphism invariant by construction Il6l 110511 . Diffeomorphism 
invariance encompasses the requirement that the laws of physics as experienced by 
an internal observer should not depend in any way on external or predefined struc- 
tures. In particular, the lesson from general relativity is understood to be that there 
should be no prior geometry, i.e., that the geometry is fully defined by the matter 
distribution. 

How the terms internal observer and external structure in the previous descrip- 
tion should be interpreted in the framework of the condensed matter model will 
become clear in the remainder of this chapter. But there is already an important 
point to be made here. If analogue gravity has reminded us that general relativity 



2 General covariance simply means that the laws of physics must transform properly under ar- 
bitrary coordinate transformations, which is essentially equivalent to the (physically irrelevant) re- 
quirement that they can be written in tensorial form. General or diffeomorphism invariance means 
that the laws of physics should be invariant, i.e.: not transform at all, under arbitrary coordinate 
transformations. See, e.g., J 134 1 for a discussion. Of course, terminology should not become sacro- 
sanct. However, at least historically the confusion between both principles has been more than just 
a linguistic issue Il94l . 
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encompasses two different aspects, kinematical and dynamical ones, then perhaps 
the same distinction should be made when discussing diffeomorphism invariance. 
Indeed, it is in principle perfectly possible to have a system in which part of the 
physical laws — in particular, the equations of motion governed by the geometry 
(taken as instantaneous or stationary) — are diffeomorphism invariant, whereas an- 
other part — the dynamics that governs the evolution of this geometry — is not. We 
will therefore discuss diffeomorphism invariance in condensed matter systems at 
these two levels: the kinematical level and the dynamical level. 



8.4 Diffeomorphism invariance at the kinematical level 

In general relativity, given a particular geometry, there are many metrics that cor- 
respond to this geometry. One could formalise this by saying that the geometry 
is defined by \gyv\, the equivalence class of all diffeomorphism equivalent metrics. 
The question with respect to condensed matter models is then whether the effective 
spacetimes that are obtained are really described by geometries in the sense just de- 
fined, rather than by individual metrics that are distinguishable due to the presence 
of the background structure. 

To reformulate the question in terms of background (in)dependence, remember 
that in condensed matter models, we are dealing with bi-metric systems. The atoms 
of the condensate can be thought of as obeying the Newtonian laws of physics, and 
in particular Galilean relativity, since v <C c^ght for the velocity v of the condensed 
matter fluid. The excitations inside the condensed matter systems, on the other 
hand, see a relativistic geometry. This relativistic geometry is in a certain sense 
embedded in the laboratory frame, which provides a fixed background, and hence 
a privileged coordinate system. Then, if background independence is apparently 
violated, how could such a setup lead to a diffeomorphism invariant theory? 

A naive answer to these questions is the following. An observer inside, for exam- 
ple, a submarine could perform an 'acoustic' experiment of the Michelson-Morley 
type [88J in the water surrounding him to find out whether the submarine is at rest 
or moving with respect to the water. An acoustic signal should here be interpreted 
as any signal in the form of a relativistic massless field perturbation, moving at the 
speed of sound c s (rather than the speed of light Ci ig ht)|^] Then, in the simplest form 
of the Michelson-Morley experiment^] the observer would take an interferometer 

3 In this section, we will explicitly differentiate the speed of sound c s and the speed of light Ci ig ^t 
in order to avoid any possible confusion. 

4 Modern versions such as [91 1 are of course technically much more refined than the original 
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with two perpendicular arms of length L and send acoustic signals along the arms. 
Mirrors at the end of the interferometer arms would reflect the signals back to a com- 
mon point, where the observer looks at the interference pattern. If the interference 
fringes move when the interferometer is rotated, then the velocity of the submarine 
is not equal along the directions of both arms, and hence the observer concludes 
that the submarine moves with respect to the 'absolute' reference frame provided 
by the water. So, the observer can determine his state of movement with respect 
to the water, and diffeomorphism invariance is broken by the privileged reference 
frame provided by the background fluid. 

However, the previous argument is not a satisfactory answer to the problem that 
we are discussing in the context of condensed matter models. The reason lies essen- 
tially in the characteristics of the observer, and of the building blocks available for 
him to carry out his experiment. To refine the answer, we need to go back to the con- 
cept of a bi-metric theory. The exterior world is the world of the laboratory, and of its 
Newtonian physics, where the basic building blocks are the condensate atoms and 
their microscopic physics. The interior world consists of everything that can be de- 
scribed in terms of collective excitations (phonons, quasi-particles) emerging inside 
the condensed matter system. An internal observer is then somebody who is limited 
to the observation of the interior world, and hence to the manipulation of these col- 
lective excitations emerging inside the condensed matter system. If the observer in 
the submarine is really an internal observer in a condensed matter world, the build- 
ing blocks available for him to build an interferometer would necessarily have to be 
made out of these same collective excitations emerging inside the condensed matter 
system. The fundamental signalling velocity of those building blocks is the speed of 
sound c s of the fluid system rather than the speed of light c^g^ of the laboratory. 

Let us work out this idea and rebuild the Michelson-Morley experiment from 
scratch for an internal observer inside a condensed matter system (see [21J for more 
details). We will take the simplest case of a homogeneous medium at rest in the 
laboratory frame. The associated internal or effective metric is then Minkowskian. 
The internal observer tries to build a 'quasi-interferometer' from quasi-particles and 
wants to obtain arms that are as rigid as possible. 

When Einstein introduced special relativity [48J, he postulated rigid bars (and 
perfect clocks) without further ado, although he has later repeatedly declared his 
dissatisfaction with the aprioriness of this assumption ||3T| . Let us see what kind of 
structures an internal observer such as the one we are contemplating would have to 



Michelson-Morley experiment, but the argument given here is unaffected by the concrete implemen- 
tation. 
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look for, when trying to construct rigid bars for his experiment. 

Phenomenologically speaking, we know that rigid structures get their rigidity 
from the fact that their constituents are arranged in a regular and stable way In par- 
ticular, the equilibrium distance uq between two constituents (atoms, for example) 
is determined by a minimum in the interaction energy E(a) between both atoms: 



Rigidity assumes a sharp minimum at flo i n the E(a) function, such that a small 
deviation 5a from the equilibrium distance uq would lead to a large energetic dis- 
advantage AE. Generally speaking, this corresponds to the requirement that the 
interatomic potential must be, for example, of the Lennard-Jones type. Such po- 
tentials are typically realised in nature by the dipolar interaction between globally 
neutral objects composed of charged substructures. At large distances a ^> uq, the 
forces between the objects will be approximately zero. When they are brought closer 
together, but still a > uq, opposite charges will rearrange themselves and the two 
objects will attract each other, with a force which becomes stronger as a decreases. 
However, if the objects are forced to a separation a < ao, then the charges will re- 
distribute and due to the repulsion between like charges, a strong repellent force 
between the objects will come into play. 

Note that we are using the term 'charge' here in a generalised way. Any charge 
with the characteristics just described serves our purpose. In particular, in 3 He- 
A, chiral fermionic quasiparticles arise at low temperature together with an effec- 
tive electromagnetic field A v ||124| . This effective theory is described by a U(l) 
charge symmetry identical to the one of quantum electrodynamics, but with mass- 
less fermions. This effective, 'quasi-electrical' charge is a perfect example of what 
we have in mind. 

Since the observer wants to keep the bars at rest within the system, we fur- 
thermore assume that he has massive charged quasi-particles at his disposal. This 
mass could for example originate in multi-component condensed matter systems 
(see [120J for the case of a two-component BEC), in quasiparticles with an intrinsic 
internal structure (e.g., a small pipe through which the background fluid is pumped Il89| ), 
or from a symmetry-breaking mechanism like the Brout-Englert-Higgs mechanism ||50l 
|63| which is generally assumed to provide mass to the particles of the standard 



5 It is perhaps useful to point out that the idea for the Brout-Englert-Higgs mechanism is an ex- 
cellent example of a condensed matter analogy applied in high-energy physics, since it was inspired 
precisely by the effective mass acquisition process in superconductors [4l l92l . 




dE 



(8.4.1) 
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So assume charged massive 'elementary quasi-particles' are available to the in- 
ternal observer, which he can use as building blocks to compose globally neutral 
constituents. Also, assume that these quasi-particles interact through the use of 
photon-like signals moving at a velocity c s , which can be described as relativistic 
collective excitations corresponding to an effective field A^, as in the case of 3 He- 
A mentioned before. Furthermore, we assume that if there are different types of 
elementary quasi-particles, they all interact with signals propagating at the same ve- 
locity c s . This condition is related to the principle of equivalence in general relativity, 
and can easily be fulfilled if one assumes (as in the BEC model discussed previously 
throughout this thesis) that the velocity of these wave signals is determined by the 
properties of the background condensate, and not by the intrinsic properties of the 
signals themselves. 

One might argue that we of course know that all these conditions are satisfied, 
since what we here call the quasi-particles created by the hypothetical microscopic 
condensed-matter-like system are precisely what we, as internal observers of our 
universe, perceive as the 'real' particles of the standard model, and hence their con- 
glomerates interact through electromagnetic signals. The point is that any type of 
interaction with any type of charge and any type of mass-acquisition process would 
do, as long as there is a constant signalling speed (or at least an invariant limit- 
ing speed) and a Lennard-Jones type of potential. Hence the construction that the 
internal observer in the condensed-matter system adopts to build rigid bars will 
inevitably lead to equations of exactly the same form as those that we know from 
classical electrodynamics, as we will now illustrate. 

Let us first consider a single massive quasi-particle with a charge q, acting as a 
source for the relativistic field A^. This field will then necessarily obey the following 
equations: 



DA H -d H (d v A v ) =j H ; (8.4.2) 
]]i = {-^ 3 [x-x(i)]^(v/c s )^ 3 [x-x(0] •} (8.4.3) 

For a source at rest, and taking the Lorenz gauge (d^A^ = 0), the solution is 



Now assume that the source moves at a velocity v, for example in the x direction. 
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Then, because of the finite speed of signal propagation c s , the solution becomes 



[( 7s (x - vt) - 7s x ) 2 + (y - y ) 2 + (z - zo) 2 ] 1 / 2 ' 

MX) = [(7s(x - vt) - ls xtf + (y - y ) 2 + (z - z ) 2 ]V2 ; (8A5) 
A y (x) = A z (x) = . 

From this solution, we see that the fields decay faster in the x direction than in the 
orthogonal y and z directions. The ratio between both decays is given by an acoustic 
or sonic Lorentz factor 

7s = (l-z; 2 /c 2 r 1/2 . (8.4.6) 

Since the internal observer is building an interferometer arm, he will need to add 
more quasi-particles. If a single quasi-particle acquires a 7 S factor in its direction of 
movement, then two particles in movement, aligned along the direction of move- 
ment, will experience a modified interaction energy potential given by E'(a) = 
£(7 s a), where the prime denotes quantities in the co-moving reference frame. As 
before, the equilibrium distance is given by the minimum of this potential: 



dE'(a) 
da 



E' Q ,a' Q . (8.4.7) 



We obtain 



dE'(a) dE(a') , 

= = 7s -jV 2 ; a = 7s« , (8.4.8) 

da da' 

which shows that the minimum now occurs when a' = ao, i.e., when a = 7<T lfl o- So 
the 'real' distance (in the sense of the distance measured in the absolute reference 
frame provided by the laboratory) between the two quasi-atoms has decreased by 
an acoustic Lorentz factor 7 S due to their velocity with respect to the medium. 

This factor 7 S is precisely the length contraction (the contraction of an interfer- 
ometer arm oriented in the direction of motion) needed for the interference pattern 
to remain unaffected by a uniform motion. The conclusion is then the following. A 
Michelson-Morley type of experiment using a quasi-interferometer does not allow 
internal observers to distinguish between rest and (uniform) movement with respect 
to the condensed matter medium (remarks to this effect were also made in ||811ll24| V 

We have explicitly worked out the case of a homogeneous condensate, in which 
the internal geometry is Minkowskian. But the essence of the argument is the same 
for more complicated situations: the internal observer has no way of connecting to 
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the microscopic external world, at least not using only internal geometrical tools. 
Therefore, from an internal point of view, there is no way of connecting to the back- 
ground structure and discriminating between the various metrics that describe the 
same given geometry. In other words: diffeomorphism invariance is effectively re- 
alised at the kinematical level. 

Note that we have assumed that the wave perturbations propagate at the con- 
stant speed of sound c s . As we have seen in the case of BECs, high-frequency per- 
turbations might acquire frequency-dependent speeds of propagation, which incor- 
porate sub- or supersonic corrections with respect to the (low-frequency) speed of 
sound c s . Therefore, diffeomorphism invariance at the kinematical level is respected 
to the same order at which Lorentz invariance is, but can be croken at high frequen- 
cies. In this sense, diffeomorphism invariance might be 'emergent', at least at the 
kinematical level, in the sense of being valid at the effective or internal level, but not 
necessarily at the fundamental level, i.e., in the external world. To put it differently: 
if local Lorentz invariance is a low-energy effective symmetry, violated at some high 
energy, then the same is true for diffeomorphism invariance at the kinematical level. 



8.5 Diffeomorphism invariance at the dynamical level 

In the previous section we discussed diffeomorphism invariance at the kinematical 
level. This requirement implies that the internal phenomena or observers cannot see 
the external degrees of freedom, but only the internal geometry. We argued that this 
requirement is fulfilled in condensed matter models for emergent gravity. Similarly, 
diffeomorphism invariance at the dynamical level can be understood as the require- 
ment that the dynamical equations for the geometry only depend on the internal 
degrees of freedom, and not on any external or pre-existing degrees of freedom. 

The simplest dynamical equations that are diffeomorphism invariant and de- 
scribe gravity in a way which agrees with observation are the Einstein equations [136J 
The dynamical equations obeyed by the geometries of condensed matter systems are 
not diffeomorphism invariant, they are hydrodynamical. For example, in the case 
of BECs, the Gross-Pitaevskii equation depends on external or microscopic degrees 
of freedom such as the interatomic potential and the scattering length. 

More in general, since it is a key requirement of quantum gravity to reproduce 
the Einstein equations in some adequate limit, it is essential within the condensed 
matter approach for emergent gravity to understand how the Einstein equations 
could be recovered, and if this is not possible, to at least understand why. This 
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is largely an open problem, although some promising steps have been taken, see 
e.g. |[129J. Here, we will limit ourselves to some basic remarks with respect to the 
possibility of obtaining the Einstein equations based on an induction scheme a la 
Sakharov. We emphasise that, from the point of view advocated here, the problem 
of dynamical diffeomorphism invariance and the problem of recovering the Einstein 
equations are essentially identical (up to higher-order curvature terms). 

The basic idea behind Sakharov's induced gravity proposal [106] can be under- 
stood as follows [|26l [119]. It is well known that the classical Einstein equations 
G^v = kTji V (with k = 8nG/c A ) can be obtained via an extremal procedure from 
an action 

S = Sa + Seh + S m , (8.5.1) 

where A is the cosmological constant, Seh the Einstein-Hilbert action and S m the 
action for the matter. In particular, in units such that h = c = 1: 



d*x. 



R 
16nG 



(8.5.2) 



with g the determinant of the metric, R the Ricci scalar and C m the matter La- 
grangian. 

By analogy, one can seek a quantity W, called the effective action, such that min- 
imising it with respect to the inverse metric g^ v reproduces the semiclassical Einstein 
equations 

G }iv = k{T }iv } , (8.5.3) 

in which the gravitational field can be treated classically even though the matter 
fields may be quantised. 

Assume that a Lorentzian manifold with metric g^ v is given, with quantum fields 
propagating on it, but without any prior dynamics. Then, such an effective action 
can indeed be constructed at one-loop order, based, e.g., on path-integral quanti- 
sation techniques 11261 . For a scalar field with mass m, for example, this one-loop 
effective action Wi can be written as I1119H : 

W 1 = -^\ndet{u + m 2 + xR) = -^Trln (u + m 2 + ^r) , (8.5.4) 

where x represents the coupling between the scalar and the gravitational field. 

One can isolate the potential high-frequency divergences in W\ by using a reg- 
ularisation technique, for example the introduction of an explicit ultraviolet (UV) 
cut-off a uv , and expanding this action in orders of the curvature. Comparing this 
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to the generalised Lagrangian for Einstein gravity (including all allowed invariant 
curvature-squared correction terms) 

R 



J d 4 x^g 



-A - + KiC^C^ + K 2 R 2 + £ t 



(8.5.5) 



where Cp lvpa is the Weyl curvature tensor, one sees that regularising the divergences 
of the one-loop effective action Wj automatically induces a term which can be inter- 
preted as the cosmological constant, one which provides the Newton constant, and 
curvature-squared terms in both C^ V p a O lv P a and R 2 . In particular, such a compari- 
son leads to 



A = A -J^(-l) F f A ((x UV/ m i/ g FV ) + (8.5.6) 

i 

E(- 1 ) F /i/G( a uv,m ! ,^ v ) + ... , (8.5.7) 



G G 



for the cosmological and gravitational 'constants' A and G, plus similar equations 
for the second-order curvature terms II119I1 . In these equations, the subindex refers 
to an arbitrary reference metric go (for example, Minkowski), the summation is over 
all particle types i, with the factor (— 1) F indicating that a relative minus sign should 
be added for fermions, and the ellipses represent non-UV-divergent terms. 

Several strategies are possible on what to do with this result. Sakharov's own 
proposal 111061 implicitly assumes dominance of the one-loop terms such that Ao and 
1 / Go can be neglected in the above equations. Furthermore, he explicitly suggests 
using the experimentally known value of G to estimate the ultraviolet cut-off a uv , 
and emphasises that the magnitude of the gravitational interaction depends on the 
masses of the particles and their equations of motion (through the original metric). 
An alternative strategy 111 1911 would be to aim for one-loop finitess of A and G (and, 
if possible, also of the second-order curvature terms), and use this requirement to 
extract conditions on the particle families i. However, this turns out to produce 
unrealistic constraints, certainly not satisfied by the standard model and its most 
common extensions. This is mainly due to the fact that the obtention of a finite 
A requires not only a compensation between the number of fermionic and bosonic 
particle families, but also an incredible fine-tuning between their respective masses. 

In any case, for our present purposes, the key point is that the Einstein equations 
(plus possible higher-order curvature correction terms) can be obtained automati- 
cally from quantum field theory via an integration over the matter degrees of free- 
dom. This idea seems very appealing in the present context, and agrees well with 
the basic idea of gravity as a macroscopic emergent phenomenon^] Moreover, given 



"Sakharov himself talks of a "metrical elasticity" of space 1 106 
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the argument discussed in section 8.2 that the cosmological constant is regulated by 
thermodynamical principles, the fine-tuning problem just mentioned with respect 
to the one-loop finiteness of A would also automatically be solved. 

However, let us repeat the assumptions stated above. We assumed a Lorentzian 
manifold with metric g^ v , irrespective of how this is obtained, but without any prior 
dynamics. An Einstein-like dynamics is then spontaneously induced. In condensed 
matter models, the condition of no-prior-dynamics is violated due to the presence of 
hydrodynamics. An apparent way out would be to define an 'abstract' condensed 
matter system without any prior dynamics. However, in the approach discussed 
in this thesis, the effective Lorentzian spacetime is precisely obtained from the hy- 
drodynamic equations (from the Gross-Pitaevskii equation in the case of BECs, or 
more generally from the Euler and continuity equations, see chapter [2b. Therefore, 
a direct extension of the current setup does not seem to resolve the problem. One 
would need to redefine a microscopic system from scratch, in such a way that the ef- 
fective spacetime is obtained without needing to appeal to hydrodynamics, thereby 
allowing the Einstein equations to be obtained through Sakharov's mechanism. 

For example, in [fTOl HI] , it was shown that an effective Lorentzian spacetime 
can be obtained quite generally from the linearisation of a first-order Lagrangian 
jC((p, 9^</>) in terms of a single scalar field (p, and even from a wide class of second- 
order hyperbolic partial differential equations, without recurring to any Euler or 
continuity equation. A formal condition was also derived for the induced dynam- 
ics of the system to be 'Einstein-like'. However, this formal condition seems to re- 
quire some fine tuning of which it is not clear how it could be implemented phys- 
ically, while there also arise serious technical complications when considering sev- 
eral fields ||T2f. Alternatively, the following reasoning might provide interesting 
clues. 

In the context of fermionic systems, it has been argued that an induction pro- 
cedure in the spirit of Sakharov's proposal can nevertheless be carried out, even in 
hydrodynamic systems ||124U129H l30l. Moreover, it is argued that the resulting dy- 
namics depends on the relation between the characteristic scales in the system. More 
precisely, in some fermionic systems such as 3 He-A, apart from the Lorentz violation 
scale Ei, a second characteristic energy scale E* appears, which determines the scale 
below which part of the fermions can form a bosonic condensate and the quasipar- 
ticles see an effective gravitational field. This energy scale E* can be interpreted as 
(the analogue of) the Planck scale (E* = E P {) [|23l [73| . As in the case of BECs (see 
the brief discussion in chapter [6]), these two scales Ei and E* are in principle un- 
related. According to the argument in ||124t |129t 11301 , the dynamical laws for this 
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gravitational field can be obtained by an induced-gravity procedure in the spirit of 
Sakharov's proposal. 

In general, this will generate the Einstein equations plus non-diffeomorphism- 
invariant higher-order terms in powers of E*/El- Since for all known real labo- 
ratory condensed matter systems, E* ^> El, these higher-order terms are actually 
dominant and the dynamical equations satisfied by the laboratory geometries are 
not diffeomorphism invariant, but hydrodynamical. If however E* <C El, then the 
induction procedure would produce the Einstein equations plus higher-order per- 
turbative corrections in E*/E^. 

Note that this argument is in line with the experimental constraints mentioned in 
chapter [6j namely that extrapolations from recent cosmic ray and other high-energy 
experiments show that there exist stringent bounds on the most commonly expected 
types of Lorentz violation at the Planck level W68ll83l , and hence, if Lorentz violations 
occur at all, they probably occur at energies much higher than the Planck scale: 
El E P j. This observation adds strength to the argument that we just sketched, 
and reinforces the idea of treating gravity in an emergent framework. However, it 
is probably fair to admit that we need a more detailed understanding of this line of 
argumentation before any strong conclusions can be drawn from it with respect to 
the naturalness of the emergence of the Einstein equations. 

As a final note, the previous reasoning was defended in the context of systems 
where the fundamental degrees of freedom are fermionic, leading to fermionic quasi- 
matter and collective bosonic fields. However, an induction procedure such as 
Sakharov's can in principle also be carried out in a bosonic system. Whether the 
Einstein equations can indeed be derived in a purely bosonic system is largely an 
academic question, since the degrees of freedom of the effective matter would then 
in principle also be purely bosonic, contrarily to the case of our universe. However, 
in terms of the naturalness or universality of the Einstein equations, the question 
might have some relevance. We do currently not have a clear view on this question. 
On the one hand, the fact that our universe contains effective fermionic matter, and 
that the Einstein equations — at least at the classical level — essentially depend on this 
effective matter content for the interplay with the geometry, seems to indicate that 
a system with only bosons at the fundamental level (and certainly a relatively sim- 
ple system such as a BEC, let alone in the weakly interacting regimes that we have 
discussed in this thesis) will not generate sufficient complexity to reproduce the Ein- 
stein equationsF] On the other hand, however, in the light of our above argument 



7 A similar conclusion is reached by the authors of 1 59 1, who study a toy model based on a BEC but 
with massive quasi-particles, and obtain that the analogue gravitational dynamics is encoded in the 



8.6. SUMMARY 



125 



on the relation between the Einstein equations and diffeomorphism invariance, it 
might be sufficient to define a sufficiently complex diffeomorphism invariant sys- 
tem, whether fermionic or bosonic, to obtain the Einstein equations. 

In any case, although it is not clear if Bose-Einstein condensates are a relevant 
model for the dynamical aspects of gravity, they certainly have a lot to teach us 
about its kinematical aspects. 

8.6 Summary 

In this chapter, we have discussed some aspects of a model for emergent gravity 
based on the condensed matter analogy. In such a model, classical gravity is thought 
of as an effective property emerging at low energies from the collective behaviour 
of a large number of microscopic constituents or 'spacetime atoms'. The quantum- 
gravitational degrees of freedom would then not be obtainable directly (for example, 
through a quantisation procedure) from the classical degrees of freedom of gravity. 

A first motivation for such an approach is the observation that most of our uni- 
verse is extremely cold compared to the Planck scale, hence the quantum-gravitational 
degrees of freedom are probably effectively frozen out in the physics that we ob- 
serve. A second motivation comes from the accelerating expansion of our universe, 
which seems to require some form of 'dark energy'. We discussed the argument 
that, from a condensed matter point of view, the quantum vacuum energy is really 
the natural candidate for such a dark energy, and provides a framework to under- 
stand both the coincidence problem (why Ci A ~ CIm) and the cosmological constant 
problem (why the zero-point energy of the vacuum is not huge). 

The main part of this chapter was dedicated to a discussion of diffeomorphism 
invariance in emergent gravity. We argued that this issue should be split up into 
a kinematical and a dynamical aspect. With respect to the kinematical aspect, we 
showed that diffeomorphism invariance can emerge as a low-energy effective prop- 
erty at the same level as Lorentz invariance. We discussed how the dynamical as- 

form of a modified Poisson equation. In [60], the same authors propose an alternative way to obtain 
effective gravitational dynamics in an abstract toy model related in spirit to 1 10] [11] Il2l mentioned 
above. Interestingly this model is based on a fundamentally Euclidean (and hence 'timeless') theory 
but nevertheless an effective Lorentzian signature is obtained for the perturbations. Although again 
it seems that there is not a sufficient amount of complexity being generated to obtain the Einstein 
equations, the effective dynamics is nevertheless explicitly diffeomorphism invariant. In particular 
the gravitational dynamics is obtained in the form of a scalar theory of gravity namely Nordstrom 
gravity. 
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pect of diffeomorphism invariance is related to the problem of recovering the Ein- 
stein equations in a model based on the condensed matter analogy. We made some 
observations with respect to the possibility of achieving this through an induction 
procedure such as Sakharov's proposal. The essential problem, as we pointed out, is 
that Sakharov's procedure assumes the absence of any prior dynamics, whereas in 
the condensed matter analogy, the hydrodynamical equations are essential to obtain 
the emergence of an effective spacetime. 

We mentioned statements in the literature, in the context of fermionic systems 
such as 3 He-A, which argue that a procedure similar to Sakharov's approach could 
be used to calculate the effective gravitational dynamics, even in a hydrodynamic 
system, and would lead to the Einstein equations plus higher-order correction terms, 
provided that the Lorentz violation energy scale lies far above the (analogue) Planck 
scale. There exists no known laboratory condensed matter system which fulfils this 
requirement. However, it is not clear whether this is due to some fundamental rea- 
son, or whether it is precisely an illustration of the fact that quantum gravity can 
be treated as a condensed-matter-like system, albeit — precisely for this relation be- 
tween the characteristic scales — an exceptional one in comparison with laboratory 
condensed matter systems. 

Finally, we noted that such an induction procedure could in principle also be 
applicable to bosonic systems, although it is not clear at present whether this would 
be sufficient to recover the Einstein equations. 



Chapter 9 

General summary and outlook 



We have discussed various aspects of the gravitational analogy in Bose-Einstein 
condensates (BECs), and of how such an analogy can serve as a source of inspiration 
in quantum gravity both in terms of quantum gravity phenomenology and in an 
attempt to understand some more conceptual issues related to quantum gravity. 

We started with a review that led us from mean-field theory in BECs to the emer- 
gence of an effective acoustic metric for linear perturbations of the background con- 
densate in the hydrodynamic limit. We studied the validity of this hydrodynamic 
limit, and emphasised that it is violated at high frequencies and in the presence of 
a horizon. In such cases, the effective geometric or relativistic description breaks 
down, giving way to the underlying microscopic theory. This behaviour is similar 
to what is expected in gravity, with the important difference that, in the case of BECs, 
the underlying microscopic theory is well understood. Moreover low-temperature 
experiments with, for example, dilute alkali gases are in accordance with the the- 
oretical treatment of BECs to a high degree of accuracy. We also stressed that the 
full dynamical aspect of general relativity (the Einstein equations) are not included 
in the analogy, at least not in its direct application in BECs, a problem to which we 
came back in the previous chapter. 

A first, general and well-known conclusion is that some phenomena related to 
black hole physics, such as Hawking radiation, can in principle be simulated in 
BECs, both theoretically and experimentally. Several of the issues that we discussed 
in this thesis directly affect Hawking radiation and its possible experimental realisa- 
tion. Generally speaking, though, the original contribution of this thesis lies beyond 
the basic observation of the possible simulation of black-hole related phenomena. 
In particular, we focused on how the gravitational analogy in BECs can be extended 
to study some aspects of quantum gravity, both in the sense of quantum gravity 
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phenomenology and in terms of constructing a toy model for emergent gravity. 

The crucial element in the theory of BECs that makes it relevant for the study of 
quantum gravity phenomenology is the following. Beyond the hydrodynamic ap- 
proximation, the first corrections to the relativistic description in a BEC are encoded 
in the form of higher-order, superluminal modifications of the dispersion relation. 
These superluminal corrections are associated with the existence of a privileged ref- 
erence frame, and thereby lead to a breaking of the effective Lorentz symmetry. They 
therefore form an excellent testbed for quantum gravity scenarios which consider a 
similar possibility of high-energy violations of local Lorentz invariance, and in par- 
ticular of those scenarios, for example based on effective field theory, where the first 
modifications of the dispersion relation that are allowed by observations are of ex- 
actly the same kind as those that occur in a BEC. 

Two possible strategies could then be developed. The first one consists in elab- 
orating the gravitational analogy in a BEC, proposing a concrete configuration, and 
studying some of its aspects inspired by relativistic phenomena of interest. 

We applied this first strategy in the following way. We presented a simple (1+1)- 
dimensional black hole configuration, and studied its dynamical stability based on 
the full Gross-Pitaevskii equation (i.e., incorporating the superluminal corrections 
to the dispersion relation). We found that such configurations are devoid of dynam- 
ical instabilities under rather generic boundary conditions inspired directly by those 
used in the study of general relativistic black holes. This opened the way to study 
their stable dynamical modes or quasinormal modes, and moreover confirmed that, 
possible experimental complications apart, it should indeed be possible to repro- 
duce the phononic analogue of Hawking radiation in a BEC. 

With respect to the quasinormal modes of black hole configurations, we found 
that in the particular (l+l)-dimensional setup that we studied, no quasinormal modes 
could exist in the hydrodynamic al limit. We remarked that in general relativity, 
similarly, no quasinormal modes can exist in 1+1 dimensions, because of the con- 
formally flat character of (l+l)-dimensional spacetimes in combination with the 
conformal invariance of the d'Alembertian wave equation. In our BEC model, a 
continuous spectrum of short-lived quasinormal modes was found to exist beyond 
the hydrodynamical limit, i.e., when using the full dispersion relation. This result 
is surprising in the sense that in general relativity, even in 3+1 dimensions, the 
quasinormal modes form a discrete spectrum. We pointed out that the origin of 
these quasinormal modes lies in the superluminal character of the dispersion re- 
lation, which leads to a permeability of the (zero-frequency) horizon. Ultimately, 
this means that these quasinormal modes are characteristic of a relaxation at the 
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microscopic level, which expresses itself in the form of lifetimes of the order of the 
Lorentz violation scale (determined by the healing length, and hence ultimately by 
the scattering length, in the case of BECs). We therefore speculated that a similar 
continuous spectrum of short-lived 'planckian' quasinormal modes could show up 
in the spectrum of gravitational black holes, if quantum gravity also turns out to 
incorporate superluminal high-energy corrections. 

A second strategy consists in directly modifying the physics that describes phe- 
nomena related to black holes in a way which incorporates superluminal corrections 
at high energies. We applied this strategy to the Hawking radiation of a dynamically 
collapsing configuration. 

We stressed the fact that the superluminal dispersion causes the effective speed 
of light to become frequency-dependent, and hence also the horizon to become a 
frequency-dependent concept. In particular, the horizon for higher frequencies lies 
inside the zero-frequency or classical horizon, and moves towards the singularity as 
the frequency increases. As a corollary, several quantities associated with the hori- 
zon, such as the surface gravity, can also become frequency-dependent. Moreover, 
the superluminal dispersion means that, at every moment of the collapse there is 
a critical frequency such that higher frequencies simply do not experience a hori- 
zon. If the collapse saturates at some level — in other words: unless a singularity is 
formed at the end of the collapse — then there also exists a global critical frequency. 
We stressed the fact that the violation of Lorentz invariance, and in particular the su- 
perluminal character of the dispersion relation, implies that the usual equivalence 
between a collapsing configuration and a stationary black-hole configuration (with 
adequate initial conditions on the horizon) is no longer valid, at least not in a way 
straightforwardly adapted from the Lorentz-invariant case. An easy way to see this 
is that, because of the frequency-dependence of the horizon just mentioned, such a 
construction would require imposing conditions arbitrarily close to the singularity. 

Three important consequences distinguished the outcome of our model from the 
standard (Lorentz-invariant) Hawking result. First, we saw that the existence of the 
critical frequency makes the overall radiation weaker than the standard Hawking 
spectrum. Second, the frequency-dependence of the surface gravity implies that 
the radiation spectrum can also undergo a modification in its form. In particular, 
if the Lorentz violation scale lies sufficiently above the critical frequency, important 
ultraviolet contributions can show up which are associated with frequencies whose 
effective horizon lies well inside the zero-frequency horizon (in other words: close 
to the singularity). Finally, the radiation coming from the collapse process is not 
stationary, but dies out with time. 
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It is probably worth stressing that, in contrast to the extensive body of literature 
on the transplanckian problem in stationary black holes, our work provides the first 
calculation of the sensitivity of Hawking radiation to the influence of superlumi- 
nal dispersion relations in a collapsing configuration. Moreover, the results that we 
discussed were obtained analytically, i.e., we have used numerical integration to il- 
lustrate them, but they can qualitatively be derived from the analytical formula that 
we obtained for the late-time radiation. This has necessarily involved some (phys- 
ically motivated) approximations, in particular we have only traced the evolution 
of the original (Minkowskian) degrees of freedom. It would certainly be interest- 
ing to see whether a consistent treatment can be given which takes into account the 
possibility that the number of degrees of freedom abruptly increases during the col- 
lapse process, and how this would affect our results. One result that seems firmly 
established, though, is that the equivalence between collapsing configurations and 
stationary ones is no longer valid in the case of superluminal dispersion, at least not 
in a form straightforwardly extended from the standard case. In any case, apart from 
the intrinsic interest of our results, this should also be considered as an invitation to 
further study the robustness of Hawking radiation in collapsing configurations. 

In the final part of this thesis, we illustrated how the basic ideas of the condensed- 
matter analogy can be applied to study some more conceptual problems related to 
quantum gravity. We gave several motivations for such an 'emergent gravity' ap- 
proach, in which classical gravity is treated as a low-energy effective description 
emerging from the collective behaviour of a large quantity of microscopic degrees 
of freedom ('spacetime atoms'). In particular, we mentioned two reasons to take 
such an approach seriously: the ultra-cold temperatures in most of the universe 
(compared to the Planck scale) and the accelerating expansion of the universe, i.e., 
dark energy. 

We then discussed the question of diffeomorphism invariance in such a model 
for emergent gravity based on the condensed matter analogy. We argued that, at the 
kinematical level, diffeomorphism invariance is realised as a low-energy, effective 
symmetry, in the same way as Lorentz invariance. The heart of our argument con- 
sisted in pointing out that an internal observer inside a condensed matter system, 
confined to the manipulation of effectively relativistic acoustic signals and quasi- 
particles, would have no way to detect his absolute state of movement with respect 
to the laboratory frame. In other words, at the kinematical level, the internal laws of 
physics do not depend on any prior or external geometry, and so diffeomorphism 
invariance is indeed realised at this level. 

At the dynamical level, the problem is somewhat more complicated. We argued 
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that the question of dynamical diffeomorphism invariance is essentially identical to 
the problem of recovering the Einstein equations, and briefly discussed a possible 
strategy, namely an induction procedure a la Sakharov. We pointed out that a crucial 
element in this respect seems to be the relation between (the analogue of) the Planck 
energy scale E P i and the Lorentz violation scale El- The fact that condensed matter 
systems do not reproduce the Einstein equations might simply be due to the fact 
that all known real laboratory condensed matter systems have E P i ^> El. It remains 
to be seen whether a fundamental reason can be given for this relation. 

The global conclusions of our thesis are then the following. 

First, we have confirmed that there is no theoretical obstacle to the simulation of 
Hawking radiation in a BEC. Indeed, dynamically stable black-hole configurations 
are in principle possible. 

Second, we have seen that the superluminal corrections to the dispersion rela- 
tion can cause strong qualitative modifications to relativistic phenomena. In the 
case of quasinormal modes, these modifications were quantitatively small. How- 
ever, in the case of Hawking radiation in a collapsing configuration, we found that 
the changes could be important both qualitatively and quantitatively. It is perhaps 
worth pointing out that we avoided at all times to make assumptions about the gen- 
eral relativistic singularity in a gravitational black hole. Therefore, our results can 
be interpreted as a consequence of the influence of microscopic corrections to the 
horizon only, not to the singularity. In this sense, the conventional wisdom which 
states that quantum gravitational effects near the horizon can in general (i.e., except 
for extremely small black holes) only lead to small modifications with respect to 
classical gravity — namely, modifications of the order of the Planck length divided 
by the radius of the black hole — should be taken with some caution. Indeed, the 
modifications that we have found are essentially applicable to black holes of any 
size, provided that the dispersion relations acquire superluminal modifications at 
high energies. There is no need to overthrow the common wisdom to interpret this 
: the key issue is that superluminal modifications imply that there will always be 
some high frequencies which experience the black hole as extremely small, thereby 
introducing important corrections, while other, even higher frequencies simply see 
no black hole at all. 

Finally, a third global conclusion is that it seems very worthwhile exploring how 
far the idea that classical gravity might emerge from the collective low-frequency 
behaviour of some microscopic degrees of freedom, in a sense related to thermody- 
namics, can be taken, and perhaps even should be taken if we really want to under- 
stand some aspects related to quantum gravity. 
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CHAPTER 9. GENERAL SUMMARY AND OUTLOOK 



An outline for future work is tentatively the following. 

First, several experimental condensed matter groups are showing an increas- 
ingly concrete interest in a possible experimental detection of Hawking radiation. 
Such an experiment would, at least in a first instance, be based on a stationary con- 
figuration, or on an externally imposed evolution such that the final state is enforced 
to resemble a stationary black hole. We can of course only be extremely enthusiastic 
about the possibility of such an experimental confirmation. At the same time, how- 
ever, we strongly feel that more effort is needed to understand the precise relation 
between such stationary or externally enforced configurations and the dynamical 
formation of a black hole through a collapse process, before extracting any general 
conclusions from such an experiment about collapsing black holes, or before pro- 
claiming the universal robustness of Hawking radiation. It might for example be 
interesting to go beyond the particular model that we have studied and examine the 
possible influence of superluminal dispersion on collapsing models that are closer 
to what is expected for astrophysical black holes and/ or that could be simulated in 
a laboratory. As we mentioned before, the work presented in chapter [7| apart from 
the intrinsic interest of the results obtained, should therefore be considered as an 
onset and an invitation to further study these aspects. 

Second, while important technical progress has been made in the traditional the- 
oretical or 'top-down' approaches to quantum gravity, none of these currently make 
any clear and definite predictions with respect to, for example, semiclassical correc- 
tions to general relativity. The approach to quantum gravity phenomenology that 
we have applied in this thesis therefore seems particularly relevant in bridging the 
gap between classical gravity and a possible 'full' theory of quantum gravity, since 
it is based on a concrete and real physical example of how effective low-energy sym- 
metries can be modified at high energies. It would, for example, be interesting to 
see if any quantitative predictions can be made with respect to the modification of 
the quasinormal mode spectrum due to possible superluminal corrections in an as- 
trophysically realistic situation involving black holes. 

Third and finally, and in our personal opinion perhaps most relevant, we believe 
that several of the problems encountered in quantum gravity are not (or at least: not 
only) of a technical kind, but are conceptual in origin. For example, none of the main 
approaches to quantum gravity offers a satisfactory explanation for the accelerating 
expansion of the universe. It therefore seems to us that there remains a lot to be 
learned from the condensed matter analogy, and in particular from how it can guide 
us in the construction of a model for emergent gravity. In this respect, the most 
urgent point is perhaps to acquire a more detailed understanding of the possible 
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emergence of the Einstein equations from the collective behaviour of microscopic 
degrees of freedom. 

As a final note, we have not explicitly defined what 'the BEC paradigm for emer- 
gent gravity' precisely consists in. Perhaps curiously, although the term 'paradigm' 
in the context of scientific practice was coined by Thomas Kuhn in The Structure of 
Scientific Revolutions f77\, he nowhere gives an explicit definition of it. Even more 
curiously, in the second (1970) edition of the same book, Kuhn added a postscript 
in which he sets out to clarify "the key difficulties (...) about the concept of a 
paradigm", but again fails to provide an explicit definition. The closest he arrives at 
a definition is the following: "The term 'paradigm' is used in two different senses. 
On the one hand, it stands for the entire constellation of beliefs, values, techniques, 
and so on shared by the members of a given community. On the other, it denotes 
one sort of element in that constellation, the concrete puzzle-solutions which, em- 
ployed as models or examples, can replace explicit rules as a basis for the solution 
of the remaining puzzles of normal sciences" (HZZL 1970 edition, p. 175). 

It should then come as no surprise that, just like the concept of a paradigm itself, 
the BEC paradigm for emergent gravity cannot be defined in a single, simple and 
precise sentence, but is a general framework of which we have set out to illustrate 
some of the most important aspects. However, the three points that we already men- 
tioned in the introduction (possibility of experimental simulation, complementing 
theoretical top-down approaches for quantum gravity with bottom-up approaches 
rooted in 'real' physics, and finally: take the idea of emergent gravity seriously) and 
repeated in our outline for future work perhaps give an appropriate summary for 
the BEC paradigm, and by extension for the condensed matter paradigm for emer- 
gent gravity, while at the same time indicating that it is very much a paradigm under 
construction. Then again, as the popular saying indicates, the road to success is always 
under construction, and perhaps this is also valid for the road to quantum gravity. 

As Karl Popper remarked, "Our knowledge can only be finite, while our igno- 
rance must necessarily be infinite" II103R . 
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